[yt-users] Volume rendering - filenames

[Mike Butler]

Hello,

This is probably a simple thing to fix, but I've been trying to use the
volume rendering package on my data and am running into an error.  If I
follow the steps at http://paste.enzotools.org/show/404/
when I get to mi,ma = dd.quantities["Extrema"]("Density")[0]
I get the error message at the bottom of this email.  It seems that it is
looking for files ending in .cpu000#  rather than .grid.cpu000#  like the
default enzo output.  I've tried doing a rename but then I get the same
error but the other way around.  Making a copy of the file with each
filename works if I'm just doing a single image, but I'm trying to make a
movie so doing that for 100 files for 8 cpus each would be a bit tedious.

If anyone knows a way around this or I'm doing something wrong, I'd
appreciate any suggestions.  Thanks

-Mike

In [8]: mi,ma=dd.quantities["Extrema"]("Density")[0]
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', (105, 0))

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', (72, 0))

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', (119, 0))

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

/home/butler85/yt/src/yt-trunk-svn/scripts/iyt in <module>()
----> 1
      2
      3
      4
      5

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in
__call__(self, *args, **kwargs)
     70                           self._data_source.pf.h.io)
     71         if lazy_reader and not self.force_unlazy:
---> 72             return self._call_func_lazy(args, kwargs)
     73         else:
     74             return self._call_func_unlazy(args, kwargs)

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in
_call_func_lazy(self, args, kwargs)
     77         self.retvals = [ [] for i in range(self.n_ret)]
     78         for gi,g in enumerate(self._get_grids()):
---> 79             rv = self.func(GridChildMaskWrapper(g,
self._data_source), *args, **kwargs)
     80             for i in range(self.n_ret): self.retvals[i].append(rv[i])
     81             g.clear_data()

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in
_Extrema(data, fields, filter)
    385     mins, maxs = [], []
    386     for field in fields:
--> 387         if data[field].size < 1:
    388             mins.append(1e90)
    389             maxs.append(-1e90)

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in
__getitem__(self, item)
     41         return getattr(self.grid, attr)
     42     def __getitem__(self, item):
---> 43         return self.data_source._get_data_from_grid(self.grid, item)
     44
     45 class DerivedQuantity(ParallelAnalysisInterface):

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseDataTypes.pyc in
save_state(self, grid, field)
     41         old_keys = grid.data.keys()
     42         grid.field_parameters = self.field_parameters
---> 43         tr = func(self, grid, field)
     44         grid.field_parameters = old_params
     45         grid.data = dict( [(k, grid.data[k]) for k in old_keys] )

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseDataTypes.pyc in
_get_data_from_grid(self, grid, field)
   1660         else:
   1661             pointI = self._get_point_indices(grid)
-> 1662             if grid[field].size == 1: # dx, dy, dz, cellvolume
   1663                 t = grid[field] * na.ones(grid.ActiveDimensions,
dtype='float64')
   1664                 return t[pointI].ravel()

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseGridType.pyc in
__getitem__(self, key)
    133         """
    134         if not self.data.has_key(key):
--> 135             self.get_data(key)
    136         return self.data[key]
    137

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseGridType.pyc in
get_data(self, field)
    166                     return self.data[field]
    167                 try:
--> 168                     temp = self.hierarchy.io.pop(self, field)
    169                     self[field] = na.multiply(temp, conv_factor,
temp)
    170                 except self.hierarchy.io._read_exception, exc:

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DataReadingFuncs.pyc in
pop(self, grid, field)
     55         else:
     56             # We only read the one set and do not store it if it
isn't pre-loaded
---> 57             return self._read_data_set(grid, field)
     58
     59     def peek(self, grid, field):

/home/butler85/yt/src/yt-trunk-svn/yt/lagos/DataReadingFuncs.pyc in
_read_data_set(self, grid, field)
    232     def _read_data_set(self, grid, field):
    233         return HDF5LightReader.ReadData(grid.filename,
--> 234                 "/Grid%08i/%s" % (grid.id, field)).swapaxes(0,2)
    235
    236     def _read_data_slice(self, grid, field, axis, coord):

ReadingError: ReadHDF5DataSet:
/home/butler85/sim/starsfuvamrrec/Cloud_128AMR4StarsFUV_0101.cpu0000 does
not exist, or no read permissions