
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
participants (1)
-
Christine Simpson