[yt-users] Analyzing DM-only data

[kcai at pa.uky.edu]

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

Quoting kcai at 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:
>
> data_style="enzo_packed_3d"
>
> to the StaticOutput initialization.  Britton Smith has done some work
> recently on improving this, so he might be able to chime in with more
> info.
>
> -Matt
>
> 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
> with
>
> 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
>       2
>       3
>       4
>       5
>
> /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)
>     183
>     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,
> use_colorbar=use_colorbar,
>     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
> StaticOutput.")
> --> 146             self.__hierarchy = self._hierarchy_class(self,
> data_style=self.data_style)
>     147         return self.__hierarchy
>     148
>
> /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,
> testGridID)
> --> 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
> _setup_data_queue(self)
>     137
>     138     def _setup_data_queue(self):
> --> 139         self.queue = _data_style_funcs[self.data_style][5]()
>     140
>     141     def _setup_grid_corners(self):
>
> KeyError: 'enzo_packed_3d'
>
> Thanks a lot!
> Kai
>

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



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