1 Nov
2014
1 Nov
'14
12:10 p.m.
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