[yt-users] Reading Array Data

[John Regan]

Hi All,

I'm trying to get YT to read a vector which is stored as part of particle
structure in Enzo (a vector in the active particle structure). The data
looks OK when I view the raw HDF5 data using hdfview. I find that in some
snapshots YT can read the data but not in others. If there is more than one
particle things get worse for sure.

It could be that when writing to the HDF5 file I haven't put the data in a
format that is compatible with YT.

I've attached a simple script that will fall over when trying to read the
data and a link to the enzo output (size = 12 GBs). The active particles
are in files DD0120.cpu0022 (Grid 3072), DD0120.cpu0023 (Grid 3076) and
DD0014 (Grid 3080). YT falls over in the io_handler but it's not obvious to
me whether I need to modify that code or the format that I write the data
to the HDF5 file in.

Any help/pointers would be much appreciated.

Cheers,
John

Link to dataset:
https://www.dropbox.com/sh/o6acuvia1sqyp7y/AAC0lEI5lZ6bGa0tP2-G1jLea?dl=0
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20180111/63a1132e/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SmartSinkTest.py
Type: text/x-python
Size: 1048 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20180111/63a1132e/attachment.py>