[yt-users] Images of extracted regions

[Elizabeth Tasker]

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
>