<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Kai,<div> 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.</div><div><br></div><div><br><div><div>On Dec 22, 2009, at 9:18 AM, <a href="mailto:kcai@pa.uky.edu">kcai@pa.uky.edu</a> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Sorry - forgot to change the subject of last email...<br><br>Quoting <a href="mailto:kcai@pa.uky.edu">kcai@pa.uky.edu</a>:<br><br><blockquote type="cite">Hi Matt,<br></blockquote><blockquote type="cite"> Last night I was browsing the archives of the yt-users list, and<br></blockquote><blockquote type="cite">came across this thread of emails:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">[yt-users] corrupt enzo restart file?<br></blockquote><blockquote type="cite">Matthew Turk matthewturk at gmail.com<br></blockquote><blockquote type="cite">Thu Nov 12 09:30:36 PST 2009<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Hi Stephen,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">A problem with the way yt opens up dark matter only simulations is<br></blockquote><blockquote type="cite">that the data_style has to be specified. This would be, for running<br></blockquote><blockquote type="cite">on trunk, done with the parameter:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">data_style="enzo_packed_3d"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">to the StaticOutput initialization. Britton Smith has done some work<br></blockquote><blockquote type="cite">recently on improving this, so he might be able to chime in with more<br></blockquote><blockquote type="cite">info.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-Matt<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Thu, Nov 12, 2009 at 9:27 AM, Stephen Skory <stephenskory at<br></blockquote><blockquote type="cite">yahoo.com> wrote:<br></blockquote><blockquote type="cite"><blockquote type="cite">Hi all,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">I am getting this error <<a href="http://paste.enzotools.org/show/257/">http://paste.enzotools.org/show/257/</a>> when I attempt to instantiate the grid hierarchy.<br></blockquote></blockquote><blockquote type="cite">......<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I am encountering the same "TypeError" when I tried to analyze a DM-only run.<br></blockquote><blockquote type="cite">But I use iyt (the interactive yt) - how do I add the data_style then? I tried<br></blockquote><blockquote type="cite">with<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">pf = load("RD0000/RD0000", data_style="enzo_packed_3d")<br></blockquote><blockquote type="cite">and then<br></blockquote><blockquote type="cite">In [3]: pc = PlotCollection(pf, center=[0.5,0.5,0.5])<br></blockquote><blockquote type="cite">yt INFO 2009-12-22 00:01:13,536 Created plot collection<br></blockquote><blockquote type="cite">with default plot-center = [0.5, 0.5, 0.5]<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">But when I tried to plot a slice, I got:<br></blockquote><blockquote type="cite">In [4]: pc.add_slice("Density", 0)<br></blockquote><blockquote type="cite">---------------------------------------------------------------------------<br></blockquote><blockquote type="cite">KeyError Traceback (most recent call last)<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/scripts/iyt in <module>()<br></blockquote><blockquote type="cite">----> 1<br></blockquote><blockquote type="cite"> 2<br></blockquote><blockquote type="cite"> 3<br></blockquote><blockquote type="cite"> 4<br></blockquote><blockquote type="cite"> 5<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/yt/raven/PlotCollection.pyc in<br></blockquote><blockquote type="cite">add_slice(self, *args, **kwargs)<br></blockquote><blockquote type="cite"> 180 *fig_size* in (height_inches, width_inches)<br></blockquote><blockquote type="cite"> 181 """<br></blockquote><blockquote type="cite">--> 182 return self.__add_slice(PlotTypes.SlicePlot, *args, **kwargs)<br></blockquote><blockquote type="cite"> 183<br></blockquote><blockquote type="cite"> 184 def add_slice_interpolated(self, *args, **kwargs):<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/yt/raven/PlotCollection.pyc in<br></blockquote><blockquote type="cite">__add_slice(self, ptype, field, axis, coord, center, use_colorbar,<br></blockquote><blockquote type="cite">figure, axes, fig_size, periodic, data_source, **kwargs)<br></blockquote><blockquote type="cite"> 204 coord = center[axis]<br></blockquote><blockquote type="cite"> 205 if data_source is None:<br></blockquote><blockquote type="cite">--> 206 data_source = self.pf.hierarchy.slice(axis, coord,<br></blockquote><blockquote type="cite">field, center, **kwargs)<br></blockquote><blockquote type="cite"> 207 p = self._add_plot(ptype(data_source, field,<br></blockquote><blockquote type="cite">use_colorbar=use_colorbar,<br></blockquote><blockquote type="cite"> 208 axes=axes, figure=figure,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/yt/lagos/OutputTypes.pyc in _get_hierarchy(self)<br></blockquote><blockquote type="cite"> 144 if self._hierarchy_class == None:<br></blockquote><blockquote type="cite"> 145 raise RuntimeError("You should not instantiate<br></blockquote><blockquote type="cite">StaticOutput.")<br></blockquote><blockquote type="cite">--> 146 self.__hierarchy = self._hierarchy_class(self,<br></blockquote><blockquote type="cite">data_style=self.data_style)<br></blockquote><blockquote type="cite"> 147 return self.__hierarchy<br></blockquote><blockquote type="cite"> 148<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/yt/lagos/HierarchyType.pyc in<br></blockquote><blockquote type="cite">__init__(self, pf, data_style)<br></blockquote><blockquote type="cite"> 689 break<br></blockquote><blockquote type="cite"> 690 self.__guess_data_style(pf["TopGridRank"], testGrid,<br></blockquote><blockquote type="cite">testGridID)<br></blockquote><blockquote type="cite">--> 691 self._setup_data_queue()<br></blockquote><blockquote type="cite"> 692 # For some reason, r8 seems to want Float64<br></blockquote><blockquote type="cite"> 693 if pf.has_key("CompilerPrecision") \<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/share/home/kca223/src/yt-1.5/yt/lagos/HierarchyType.pyc in<br></blockquote><blockquote type="cite">_setup_data_queue(self)<br></blockquote><blockquote type="cite"> 137<br></blockquote><blockquote type="cite"> 138 def _setup_data_queue(self):<br></blockquote><blockquote type="cite">--> 139 self.queue = _data_style_funcs[self.data_style][5]()<br></blockquote><blockquote type="cite"> 140<br></blockquote><blockquote type="cite"> 141 def _setup_grid_corners(self):<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">KeyError: 'enzo_packed_3d'<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks a lot!<br></blockquote><blockquote type="cite">Kai<br></blockquote><blockquote type="cite"><br></blockquote><br><blockquote type="cite">This message was sent using IMP, the Internet Messaging Program.<br></blockquote><br><br><br>----------------------------------------------------------------<br>This message was sent using IMP, the Internet Messaging Program.<br><br>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org<br></div></blockquote></div><br><div> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>Dr. Eric J. Hallman</div><div>NSF Astronomy and Astrophysics Postdoctoral Fellow </div><div>Center for Astrophysics and Space Astronomy </div><div>University of Colorado at Boulder </div><div>hallman (at) casa.colorado.edu </div><div>Phone: (312) 725-4626 </div><div><a href="http://solo.colorado.edu/~hallman/">http://solo.colorado.edu/~hallman/</a></div></div></div></span></div></span></div></span> </div><br></div></body></html>