I think the solution is that dark matter only runs do not have the field "Density", but only dark matter density is defined as a field. Density is a baryon field.

On Dec 22, 2009, at 9:18 AM, kcai@pa.uky.edu wrote:

Sorry - forgot to change the subject of last email...

Quoting kcai@pa.uky.edu:

Hi Matt,
  Last night I was browsing the archives of the yt-users list, and
came across this thread of emails:

[yt-users] corrupt enzo restart file?
Matthew Turk matthewturk at gmail.com
Thu Nov 12 09:30:36 PST 2009

Hi Stephen,

A problem with the way yt opens up dark matter only simulations is
that the data_style has to be specified.  This would be, for running
on trunk, done with the parameter:


to the StaticOutput initialization.  Britton Smith has done some work
recently on improving this, so he might be able to chime in with more


On Thu, Nov 12, 2009 at 9:27 AM, Stephen Skory <stephenskory at
yahoo.com> wrote:
Hi all,

I am getting this error <http://paste.enzotools.org/show/257/> when  I attempt to instantiate the grid hierarchy.

I am encountering the same "TypeError" when I tried to analyze a DM-only run.
But I use iyt (the interactive yt) - how do I add the data_style  then? I tried

pf = load("RD0000/RD0000", data_style="enzo_packed_3d")
and then
In [3]: pc = PlotCollection(pf, center=[0.5,0.5,0.5])
yt         INFO       2009-12-22 00:01:13,536 Created plot collection
with default plot-center = [0.5, 0.5, 0.5]

But when I tried to plot a slice, I got:
In [4]: pc.add_slice("Density", 0)
KeyError                                  Traceback (most recent call last)

/share/home/kca223/src/yt-1.5/scripts/iyt in <module>()
----> 1

/share/home/kca223/src/yt-1.5/yt/raven/PlotCollection.pyc in
add_slice(self, *args, **kwargs)
   180         *fig_size* in (height_inches, width_inches)
   181         """
--> 182         return self.__add_slice(PlotTypes.SlicePlot, *args, **kwargs)
   184     def add_slice_interpolated(self, *args, **kwargs):

/share/home/kca223/src/yt-1.5/yt/raven/PlotCollection.pyc in
__add_slice(self, ptype, field, axis, coord, center, use_colorbar,
figure, axes, fig_size, periodic, data_source, **kwargs)
   204             coord = center[axis]
   205         if data_source is None:
--> 206             data_source = self.pf.hierarchy.slice(axis, coord,
field, center, **kwargs)
   207         p = self._add_plot(ptype(data_source, field,
   208                          axes=axes, figure=figure,

/share/home/kca223/src/yt-1.5/yt/lagos/OutputTypes.pyc in  _get_hierarchy(self)
   144             if self._hierarchy_class == None:
   145                 raise RuntimeError("You should not instantiate
--> 146             self.__hierarchy = self._hierarchy_class(self,
   147         return self.__hierarchy

/share/home/kca223/src/yt-1.5/yt/lagos/HierarchyType.pyc in
__init__(self, pf, data_style)
   689                 break
   690         self.__guess_data_style(pf["TopGridRank"], testGrid,
--> 691         self._setup_data_queue()
   692         # For some reason, r8 seems to want Float64
   693         if pf.has_key("CompilerPrecision") \

/share/home/kca223/src/yt-1.5/yt/lagos/HierarchyType.pyc in
   138     def _setup_data_queue(self):
--> 139         self.queue = _data_style_funcs[self.data_style][5]()
   141     def _setup_grid_corners(self):

KeyError: 'enzo_packed_3d'

Thanks a lot!

This message was sent using IMP, the Internet Messaging Program.

This message was sent using IMP, the Internet Messaging Program.

yt-users mailing list

Dr. Eric J. Hallman
NSF Astronomy and Astrophysics Postdoctoral Fellow                    
Center for Astrophysics and Space Astronomy 
University of Colorado at Boulder  
hallman (at) casa.colorado.edu 
Phone: (312) 725-4626