[yt-users] Analyzing DM-only data

kcai at pa.uky.edu kcai at pa.uky.edu
Tue Dec 22 10:58:49 PST 2009


Eric,
     No, 'Dark_Matter_Density' does not work, either...
-> similar error:

In [5]: pc.add_slice('Dark_Matter_Density', 0)
Exception AttributeError: "EnzoHierarchy instance has no attribute  
'_data_file'" in <bound method EnzoHierarchy.__del__ of  
<yt.lagos.HierarchyType.EnzoHierarchy instance at 0x2aa1cfbcf8>> ignored
---------------------------------------------------------------------------
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'

I think it is I did not set the "data_style" right...

Thanks,
Kai


Quoting Eric Hallman <Eric.Hallman at colorado.edu>:

> ,
>   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 at pa.uky.edu wrote:
>
>> 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.
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
> 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
> http://solo.colorado.edu/~hallman/



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




More information about the yt-users mailing list