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_co...) 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@mcmaster.ca <mailto:taskere@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@lists.spacepope.org <mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
------------------------------------------------------------------------
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org