[yt-users] error in creating plot collection

[Christine Simpson]

Hi all,

I'm getting an error in trying to do some very simple analysis on one
output of an enzo simulation I have.  I don't have this issue with any
other output of the simulation.  Basically I can load the data and
that's it; every other routine I try on it crashes.  I have a feeling
this may be an issue with my data, not yt, but I was wondering if you yt
folks could give some insight as to what is crashing yt.  Any help would
be appreciated.  Here is a sample error where I try to create a plot
collection:

In [12]: pf = load("DD0008/output_0008")
yt         INFO       2011-02-25 17:03:41,236 Parameters: current_time
= 40.8164738158
yt         INFO       2011-02-25 17:03:41,236 Parameters:
domain_dimensions         = [128, 128, 128]
yt         INFO       2011-02-25 17:03:41,236 Parameters:
domain_left_edge          = [ 0.  0.  0.]
yt         INFO       2011-02-25 17:03:41,237 Parameters:
domain_right_edge         = [ 1.  1.  1.]
yt         INFO       2011-02-25 17:03:41,238 Parameters:
cosmological_simulation   = 1
yt         INFO       2011-02-25 17:03:41,238 Parameters:
current_redshift          = 6.36375016055
yt         INFO       2011-02-25 17:03:41,238 Parameters: omega_lambda
= 0.726
yt         INFO       2011-02-25 17:03:41,238 Parameters: omega_matter
= 0.274
yt         INFO       2011-02-25 17:03:41,238 Parameters:
hubble_constant           = 0.705

In [13]: pc = PlotCollectionInteractive(pf)
yt         INFO       2011-02-25 17:03:44,967 Getting the binary
hierarchy
ERROR: An unexpected error occurred while tokenizing input
The following traceback may be corrupted or invalid
The error message is: ('EOF in multi-line statement', (92, 0))

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call
last)

/share/home/01112/tg803911/yt_new/yt/scripts/iyt in <module>()
----> 1 
      2 
      3 
      4 
      5 

/share/home/01112/tg803911/yt_new/yt/yt/visualization/plot_collection.pyc in __init__(self, *args, **kwargs)
   1602         import pylab
   1603         self.pylab = pylab
-> 1604         super(PlotCollectionInteractive, self).__init__(*args,
**kwargs)
   1605 
   1606     def redraw(self):

/share/home/01112/tg803911/yt_new/yt/yt/visualization/plot_collection.pyc in __init__(self, pf, center)
    101         self.pf = pf
    102         if center == None:
--> 103             v,self.c = pf.h.find_max("Density") # @todo: ensure
no caching
    104         elif center == "center" or center == "c":
    105             self.c = (pf.domain_right_edge +
pf.domain_left_edge)/2.0

/share/home/01112/tg803911/yt_new/yt/yt/data_objects/static_output.pyc
in _get_hierarchy(self)
    157             if self._hierarchy_class == None:
    158                 raise RuntimeError("You should not instantiate
StaticOutput.")
--> 159             self.__hierarchy = self._hierarchy_class(self,
data_style=self.data_style)
    160         return self.__hierarchy
    161 

/share/home/01112/tg803911/yt_new/yt/yt/frontends/enzo/data_structures.pyc in __init__(self, pf, data_style)
    155             self.float_type = 'float64'
    156 
--> 157         AMRHierarchy.__init__(self, pf, data_style)
    158         # sync it back

    159         self.parameter_file.data_style = self.data_style

/share/home/01112/tg803911/yt_new/yt/yt/data_objects/hierarchy.pyc in
__init__(self, pf, data_style)
     68 
     69         mylog.debug("Parsing hierarchy.")
---> 70         self._parse_hierarchy()
     71 
     72         mylog.debug("Constructing grid objects.")

/share/home/01112/tg803911/yt_new/yt/yt/frontends/enzo/data_structures.pyc in _parse_hierarchy(self)
    217             return line.split()[2:]
    218         if os.path.exists(self.hierarchy_filename[:-9] +
"harrays"):
--> 219             if self._parse_binary_hierarchy(): return
    220         t1 = time.time()
    221         pattern = r"Pointer: Grid\[(\d*)\]->NextGrid(Next|
This)Level = (\d*)\s+$"

/share/home/01112/tg803911/yt_new/yt/yt/frontends/enzo/data_structures.pyc in _parse_binary_hierarchy(self)
    288         except h5py.h5.H5Error:
    289             return False
--> 290         self.grid_dimensions[:] = f["/ActiveDimensions"][:]
    291         self.grid_left_edge[:] = f["/LeftEdges"][:]
    292         self.grid_right_edge[:] = f["/RightEdges"][:]

ValueError: shape mismatch: objects cannot be broadcast to a single
shape