[yt-users] 2d radial profiles

[gso at physics.ucsd.edu]

I use this to plot a profile from domainleft to domainright, hope this may
help
this will plot on the x axis radius from the origin (domainleft) and y
axis the velocity

fn = "DD0000/DD0000" # parameter file to load
pf = load(fn) # load data
zmax = 1.67e-18
zmin = 9.2e-20
pc = PlotCollection(pf, center=[0.0, 0.0, 0.0])

import pylab

ray = pf.h.ray(pf["DomainLeftEdge"], pf["DomainRightEdge"])

xval = ray["RadiusCode"]
#xval = (ray["x"]**2+ray["y"]**2+ray["z"]**2)**0.5
yval = ray["Density"]

pylab.plot(na.log10(xval), na.log10(yval))

pylab.savefig("test.png")


> Hi, all--
>
> Is it possible to make a 2d radial profile from a projection with the
> existing profile tool?
>
> I think a better question is how do I make an extraction from a
> projection that I can then hand to BinnedProfile1d?
>
> Thanks,
> d.
>
> --
> Sent from my Stone Tablet and carried by my Pterodactyl.
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>