Hi John,
I think that if you want to compute "total" instead of mean you need to set
weight_field to None in create_profile.
Cheers,
Kacper
On Sat Nov 01 2014 at 6:10:50 AM John Regan
Hi All,
Is the accumulation flag working in 3.0? I tried to plot the enclosed mass in a sphere and I got some funny results.
rpm = yt.create_profile(sp, 'radius', 'cell_mass', units = {'radius': 'pc', 'cell_mass' : 'Msun'}, weight_field='density', accumulation=True, fractional=False)
print "Mass = ", rpm["gas", "cell_mass"]
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.
sp.quantities.total_quantity(["cell_mass", "particle_mass"])
Cheers, John _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org