[yt-dev] Issue #901: Derived quantity total_mass fails on DM-only datasets (yt_analysis/yt)

[Cameron Hummels]

New issue 901: Derived quantity total_mass fails on DM-only datasets
https://bitbucket.org/yt_analysis/yt/issue/901/derived-quantity-total_mass-fails-on-dm

Cameron Hummels:

In trying to calculate the mass of objects in an enzo dark matter only dataset, yt fails. It expects a ('gas', 'particle_mass') field to exist for this dataset, because that is what the field_info object thinks exists for this dataset.  I think this may be indicative of more significant problems for particle-only datasets (at least in the enzo frontend), and i was unable to fix it with some cursory changes.

My demonstration script is:

import yt
ds = yt.load('DD0046/DD0046')
ad = ds.all_data()
print ad.quantities.total_mass()

and a sample 2MB DM-only dataset for this is here: http://www.astro.columbia.edu/~chummels/DD0046.tar