[yt-users] disk Profiles

[Elizabeth Tasker]

In a fit of impatience born of being in the wrong time zone, can I wave this at people again? :)


On 2011-07-09, at 6:23 PM, Elizabeth Tasker wrote:

> Hi all,
> 
> I'm trying to create a profile object of the density in my disk as a function of radius. I don't actually want to plot it; I just need to grab an average to use as a background value elsewhere.
> 
> I think I've almost got it:
> 
> source = pf.h. disk((16,16,16), (0,0,1), 12, 15.6e-3/pf['kpc'])
> profile = BinnedProfile1D(source, 128, "RadiusCode", 0.5, 12.0)
> profile.add_fields("CellMass", weight="CellVolume")
> print profile["CellMass"]
> 
> I just have a couple of questions:
> 
> 1. Is there a way to print out the radial bins being used? i.e. the value at the centre of each bin? 
> 
> 2. RadiusCode is going to be the 3d radius, is that true? I can create my own variable for the in-disk radius, but is there a quantity already in yt for that to use with the disk profile?
> 
> Thank you~
> 
> Elizabeth
> 
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org