[yt-users] Active Particle Read Failure

[John Regan]

Hi All,

I'm having trouble getting YT3.1 to read an active particle from an Enzo
dataset. I can see the particle just fine with hdf5view so it appears to be
there but YT is falling over when it tries to read the particle.


The error I get is:


Traceback (most recent call last):
  File "CompareFieldsInCore.py", line 138, in <module>
    sphere = ds.h.sphere(centre, (10.0, 'kpc'))
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py",
line 378, in h
    self._index_proxy = IndexProxy(self)
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py",
line 129, in __init__
    ds.index
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/data_objects/static_output.py",
line 365, in index
    self, dataset_type=self.dataset_type)
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 211, in __init__
    GridIndex.__init__(self, ds, dataset_type)
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/geometry/geometry_handler.py",
line 58, in __init__
    self._detect_output_fields()
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 447, in _detect_output_fields
    ap_fields = self._detect_active_particle_fields()
  File
"/cosma/home/dp004/dc-rega4/YT/yt-x86_64/src/yt-hg/yt/frontends/enzo/data_structures.py",
line 407, in _detect_active_particle_fields
    g = gs[0]
IndexError: index 0 is out of bounds for axis 0 with size 0


As far as I can see YT isn't able to get the gridID for the particle (its
in grid 3). Any ideas what YT is missing here?

Cheers,
John

The dataset is here (~500 MBs)
https://dl.dropboxusercontent.com/u/18978426/Enzo/RD0008.tgz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160503/8832cb0b/attachment.htm>