<p>Hi Dave,</p>
<p>Looks like a leftover .yt. (and why newgridio-no?)</p>
<p>Adding a hash checking and invalidation system to the .yt files would be a great intro project for new developers and would be super useful...</p>
<p>-Matt</p>
<div class="gmail_quote">On Aug 19, 2011 5:05 PM, "david collins" <<a href="mailto:antpuncher@gmail.com">antpuncher@gmail.com</a>> wrote:<br type="attribution">> Hi, all--<br>> <br>> I have this simulation that has a single subgrid.  64^3 root, 24 zone<br>
> root grid (from 0.15625 to 0.34375 along each axis)<br>> <br>> With the old output style (enzo1.0.1) this doesn't give me any<br>> problems.  But with the new output version (packed IO, newgridio-no,<br>
> in mhdctmerge) I get the following message3<br>> <br>> <trace><br>>>>> pf = MyOutputs.FileStaticOutputFourPi('/Users/dccollins/ThisScratch/param//Turbulence/t64_1lev_m/./DD0001/data0001')<br>
>>>> pf.h<br>> Traceback (most recent call last):<br>>   File "<stdin>", line 1, in <module><br>>   File "/Users/dccollins/local/src/yt-hg/yt/data_objects/static_output.py",<br>
> line 159, in _get_hierarchy<br>>     self.__hierarchy = self._hierarchy_class(self, data_style=self.data_style)<br>>   File "/Users/dccollins/local/src/yt-hg/yt/frontends/enzo/data_structures.py",<br>
> line 158, in __init__<br>>     AMRHierarchy.__init__(self, pf, data_style)<br>>   File "/Users/dccollins/local/src/yt-hg/yt/data_objects/hierarchy.py",<br>> line 70, in __init__<br>>     self._parse_hierarchy()<br>
>   File "/Users/dccollins/local/src/yt-hg/yt/frontends/enzo/data_structures.py",<br>> line 220, in _parse_hierarchy<br>>     if self._parse_binary_hierarchy(): return<br>>   File "/Users/dccollins/local/src/yt-hg/yt/frontends/enzo/data_structures.py",<br>
> line 295, in _parse_binary_hierarchy<br>>     self.grid_dimensions[:] = f["/ActiveDimensions"][:]<br>> ValueError: shape mismatch: objects cannot be broadcast to a single shape<br>>>>> <a href="http://pdb.pm">pdb.pm</a>()<br>
>> /Users/dccollins/local/src/yt-hg/yt/frontends/enzo/data_structures.py(295)_parse_binary_hierarchy()<br>> -> self.grid_dimensions[:] = f["/ActiveDimensions"][:]<br>> (Pdb) self.grid_dimensions<br>
> array([[1, 1, 1],<br>>        [1, 1, 1]], dtype=int32)<br>> (Pdb) f['/ActiveDimensions'][:]<br>> array([[64, 64, 64],<br>>        [32, 32, 16],<br>>        [32, 32, 16]], dtype=int32)<br>> (Pdb)<br>
> </trace><br>> <br>> Why is self.grid_dimensions a 2d array?  Why don't these things match?<br>> <br>> Thanks!<br>> d.<br>> <br>> -- <br>> Sent from my computer.<br>> _______________________________________________<br>
> yt-users mailing list<br>> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</div>