I've just recently compiled enzo with grackle and run the enzo test problem: amr_cosmology.
I've loaded the results using:
ts = yt.load("/home/ian/enzo-grackle/run/CosmologySimulation/amr_cosmology/DD*/DD*")
Inspecting ts, I see that it is a DatasetSeries object, with outputs:
[DD0000, DD0000.boundary, DD0000.boundary.hdf, DD0000.configure, DD0000.cpu0000, DD0000.hierarchy, DD0000.memorymap,
DD0001, DD0001.boundary, DD0001.boundary.hdf, DD0001.configure, DD0001.cpu0000, DD0001.hierarchy, DD0001.memorymap,
DD0002, DD0002.boundary, DD0002.boundary.hdf, DD0002.configure, DD0002.cpu0000, DD0002.hierarchy, DD0002.memorymap,
As a first step in learning how to work with DatasetSeries, I'd like to generate a sequence of volume renders of density for each DD step.
I presume that this will require me to access DD0000, DD0001, DD0002...etc sequentially inside a loop or iteration, while skipping the intervening DD.boundary, DD.boundary.hdf, DD.configure, DD.cpu0000, DD.hierarchy and DD.memorymap outputs at each step.
I'm struggling to work out how to do this. Is anyone able to help me with a line or so of code to show me the way?