[yt-dev] advice on derived field debugging

[Britton Smith]

Hi all,

There seems to be an issue accessing the nuclei density fields.  These are
derived fields that give the total number density of an element with names
like "H_nuclei_density."  At some point, these have stopped working and I
am trying to debug this.  To reproduce this, one can do the following:
ds = yt.load("enzo_cosmology_plus/DD0046/DD0046")
ad = ds.all_data()
ad["H_nuclei_density"]

I have traced this as far as the following.  In the call to
create_field_info, a YTFieldNotFound exception is raised on line 137 of
yt/fields/species_fields.py.  Here, the nuclei density field function is
trying to access a number density field that should exist, but doesn't seem
to at the moment.  Later on it does, but not at this moment.  Does anyone
have any advice for where to look next?

Thanks,
Britton
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20141002/5b455ccb/attachment.htm>