[yt-users] Images of extracted regions

[Elizabeth Tasker]

I've just seen that "data_source" is an option for add_projection, but 
it still fails with "AttributeError: 'AMRCoveringGrid' object has no 
attribute '_get_cut_mask'"


Elizabeth


Elizabeth Tasker wrote:
> Hi Britton,
>
> Sorry, could you give me an example of how the call should look?
>
> I tried:
>
> p = pc.add_slice("NegEscapeVelocity", 2, field_parameters = 
> {'disk_radius':profile["DiskRadius"], 
> 'potential':profile["PotentialField"]}, data_source=extractCube)
>
> which fails because data_source is unrecognised (and indeed, isn't in 
> the function description: 
> http://yt.enzotools.org/doc/reference/api/generated/yt.visualization.plot_collection.PlotCollection.add_slice.html?highlight=add_slice#yt.visualization.plot_collection.PlotCollection.add_slice) 
>
>
> And I tried:
>
> p = pc.add_slice(extractCube, "NegEscapeVelocity", 2)
>
>
> which also failed with:
>
>
> ValueError                                Traceback (most recent call 
> last)
>
> /1/home/taskere/yt/scripts/iyt in <module>()
> ----> 1                                         
> 2                                         
> 3                                         
> 4                                         
> 5                                   
> /1/home/taskere/yt/yt/visualization/plot_collection.py in 
> add_slice(self, field, axis, coord, center, use_colorbar, figure, 
> axes, fig_size, periodic, obj, 
> field_parameters)                                                         
>    402             if field_parameters == None: field_parameters = 
> {}                                               403             obj = 
> self.pf.hierarchy.slice(axis, coord, 
> field,                                             --> 
> 404                             center=center, 
> **field_parameters)                                               
> 405         p = 
> self._add_plot(SlicePlot(                                                                      
>
>    406                          obj, field, 
> use_colorbar=use_colorbar,                                          
> /1/home/taskere/yt/yt/data_objects/data_containers.py in 
> __init__(self, axis, coord, fields, center, pf, node_name, **kwargs)
>    935         self.coord = coord
>    936         if node_name is False:
> --> 937             self._refresh_data()
>    938         else:
>    939             if node_name is True: self._deserialize()
>
> /1/home/taskere/yt/yt/data_objects/data_containers.py in 
> _refresh_data(self)
>    273         """
>    274         self.clear_data()
> --> 275         self.get_data()
>    276
>    277     def keys(self):
>
> /1/home/taskere/yt/yt/data_objects/data_containers.py in 
> get_data(self, fields)
>    780         # We take a 3-tuple of the coordinate we want to slice 
> through, as well
>
>    781         # as the axis we're slicing along
>
> --> 782         self._get_list_of_grids()
>    783         if not self.has_key('pdx'):
>    784             self._generate_coords()
>
> /1/home/taskere/yt/yt/data_objects/data_containers.py in 
> _get_list_of_grids(self)
>    994
>    995     def _get_list_of_grids(self):
> --> 996         goodI = ((self.pf.h.grid_right_edge[:,self.axis] > 
> self.coord)
>    997               &  (self.pf.h.grid_left_edge[:,self.axis] <= 
> self.coord ))
>    998         self._grids = self.pf.h.grids[goodI] # Using sources 
> not hierarchy
>
> ValueError: setting an array element with a sequence.
>
>
>
> Elizabeth
>
>
>
>
>
> Britton Smith wrote:
>> The first is definitely possible.  Just add data_source=extractCube 
>> to your call to add_projection or add_slice or whatever.  The second 
>> probably won't work because the projection or slice routine will look 
>> for other things like cell dxs that won't be present from just a 
>> similarly sized array.  If you want to do something like that, you're 
>> probably better off making a new field and then projecting that with 
>> the above data_source keyword argument in there.
>>
>> Britton
>>
>> On Thu, Aug 18, 2011 at 10:59 AM, Elizabeth Tasker 
>> <taskere at mcmaster.ca <mailto:taskere at mcmaster.ca>> wrote:
>>
>>     Hi,
>>
>>     Is it possible to make a slice or projection of a extracted, fixed
>>     resolution region? i.e. the resultant object of:
>>
>>     extractCube = pf.h.covering_grid(extract_level,
>>                                            left_edge=extractLE,
>>                                            right_edge=extractRE,
>>                                            # How many dimensions along
>>     each axis
>>                                            dims=extractDims,
>>                                            # And any fields to preload
>>     (this is optional!)
>>                                            fields=["PotentialField"],
>>                                            num_ghost_zones = 3)
>>
>>
>>     Or of an array later made to the same size? e.g.
>>
>>     cloud = na.zeros(extractCube["x-velocity"].shape, dtype='float64')
>>
>>
>>     Elizabeth
>>     _______________________________________________
>>     yt-users mailing list
>>     yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
>>     http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>   
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org