<div dir="ltr">Hi all,<div><br></div><div>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:</div><div>ds = yt.load("enzo_cosmology_plus/DD0046/DD0046")</div><div>ad = ds.all_data()</div><div>ad["H_nuclei_density"]</div><div><br></div><div>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?</div><div><br></div><div>Thanks,</div><div>Britton</div></div>