[yt-users] Active Particle Read Failure

[John Regan]

Hi Nathan,

Right you are! That was a bit lazy on my part not updating first. You guys
just move too quickly for me to keep up!

On Tue, May 3, 2016 at 2:03 PM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> You should try updating to the latest development version of yt --- yt 3.1
> is a couple years old at this point and there have been a number of
> bugfixes since then. Let us know if updating doesn't fix your problem.
>
>
> On Tuesday, May 3, 2016, John Regan <johnanthonyregan at gmail.com> wrote:
>
>> 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
>>
>>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160503/45579bba/attachment.html>