<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Sam,<div><br></div><div>Hmm, now you lay it out it seems kind of obvious. Thank you!!</div><div><br></div><div>Elizabeth</div><div><br></div><div><br><div><div>On 2011-07-11, at 11:20 PM, Sam Skillman wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hi Elizabeth,<br><br><div class="gmail_quote">On Sat, Jul 9, 2011 at 3:23 AM, Elizabeth Tasker <span dir="ltr"><<a href="mailto:taskere@mcmaster.ca">taskere@mcmaster.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi all,<br>
<br>
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.<br>
<br>
I think I've almost got it:<br>
<br>
source = pf.h. disk((16,16,16), (0,0,1), 12, 15.6e-3/pf['kpc'])<br>
profile = BinnedProfile1D(source, 128, "RadiusCode", 0.5, 12.0)<br>
profile.add_fields("CellMass", weight="CellVolume")<br>
print profile["CellMass"]<br>
<br>
I just have a couple of questions:<br>
<br>
1. Is there a way to print out the radial bins being used? i.e. the value at the centre of each bin?<br></blockquote><div><br></div><div>You can access this through profile["RadiusCode"] or profile.data._bins. Note that this will be the bin edges, and thus there will be 129 of them.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
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?<br></blockquote><div><br></div><div>
Not that I know of, but someone else may know better. </div><div><br></div><div>Best,</div><div>Sam</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<br>
Thank you~<br>
<br>
Elizabeth<br>
<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br>
_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></blockquote></div><br></div></body></html>