Hi John & Nathan,
Yes, I'm having trouble using the time series analysis machinery with my FLASH simulations. I've checked that the plot files and particle files all have the same current_time variable, so there's no issue there. I'm using a development version of yt with a changeset ID of 5cf30893afdf. This changeset includes Nathan's changes supporting cut_region functionality for FLASH plot files. Here's what I see when I try to do a time series analysis:
###########################
In [9]: ts = yt.load("mysim_hdf5_plt_cnt_*",
particle_filename="mysim_hdf5_part_*")
In [10]: for ds in ts:
print ds.current_time
....:
---------------------------------------------------------------------------
IOError Traceback (most recent call last)
<ipython-input-10-fccd01bb5065> in <module>()
----> 1 for ds in ts:
2 print ds.current_time
3
/home/emrys/yt/yt/data_objects/time_series.pyc in __iter__(self)
179 for o in self._pre_outputs:
180 if isinstance(o, string_types):
--> 181 ds = self._load(o, **self.kwargs)
182 self._setup_function(ds)
183 yield ds
/home/emrys/yt/yt/data_objects/time_series.pyc in _load(self, output_fn, **kwargs)
396 elif self._mixed_dataset_types:
397 return load(output_fn, **kwargs)
--> 398 ds = load(output_fn, **kwargs)
399 self._dataset_cls = ds.__class__
400 return ds
/home/emrys/yt/yt/convenience.pyc in load(*args, **kwargs)
84 candidates = find_lowest_subclasses(candidates)
85 if len(candidates) == 1:
---> 86 return candidates[0](*args, **kwargs)
87 if len(candidates) == 0:
88 if ytcfg.get("yt", "enzo_db") != '' \
/home/emrys/yt/yt/frontends/flash/data_structures.pyc in __init__(self, filename, dataset_type, storage_filename, particle_filename, units_override, unit_system)
209 self._particle_handle = HDF5FileHandler(self.particle_filename)
210 except:
--> 211 raise IOError(self.particle_filename)
212 # Check if the particle file has the same time
213 if self._particle_handle != self._handle:
IOError: mysim_hdf5_part_*
###########################
It looks to me like HDF5FileHandler() just does not like the wildcard. I've also tried, unsuccessfully, to feed a list of particle files to yt.load(). I thought perhaps that I should be using the DatasetSeries framework directly, but I'm not clear on how to pass in the list of particle files, as DatasetSeries does not understand the particle_filename keyword argument.