Hi John,<br>I think that if you want to compute "total" instead of mean you need to set weight_field to None in create_profile.<div>Cheers,</div><div>Kacper</div><br><div class="gmail_quote">On Sat Nov 01 2014 at 6:10:50 AM John Regan <<a href="mailto:johnanthonyregan@gmail.com">johnanthonyregan@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi All, <div><br></div><div>Is the accumulation flag working in 3.0? </div><div>I tried to plot the enclosed mass in a sphere and I got some funny results. </div><div><br></div><div><div>rpm = yt.create_profile(sp, 'radius',  'cell_mass',</div><div>                        units = {'radius': 'pc', 'cell_mass' : 'Msun'}, </div><div>                        weight_field='density', accumulation=True, fractional=False)</div></div><div><br></div><div>print "Mass = ", rpm["gas", "cell_mass"]<br></div><div><br></div><div>In this case bin[n-1] gives a mass of something like 0.5 Msun but when I print the totals quantity I get a value of several orders of magnitude higher and closer to what I would expect.</div><div><br></div><div>sp.quantities.total_quantity(["cell_mass", "particle_mass"])<br></div><div><br></div><div>Cheers,</div><div>John</div></div>
______________________________<u></u>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/<u></u>listinfo.cgi/yt-users-<u></u>spacepope.org</a><br>
</blockquote></div>