[yt-users] cut_region index error with Tipsy

Cameron Hummels chummels at gmail.com
Wed Jul 16 13:04:16 PDT 2014


Hi Samantha,

I think that data_source expects geometrical objects as opposed to cut
regions.  So you can use ds.sphere or ds.region to pass into the
ProjectionPlot, but I don't think a cut_region will work, regardless of the
frontend (gadget, enzo, etc.).  But I may be wrong here.

Cameron


On Wed, Jul 16, 2014 at 12:57 PM, Samantha Benincasa <benincsm at mcmaster.ca>
wrote:

> Hi everyone,
>
> I'm working with a Gasoline dataset.  I'm trying to create a cut region
> and then image/analyze just the particles in that set.  I'm able to
> successfully create the region using:
>
> cr = ad.cut_region(["obj['Gas','Temperature'] < 1000"])
>
> However, when I try to image that region:
>
> proj = ProjectionPlot(pf, "z",
> ('Gas','Density'),weight_field=('Gas','Density'),data_source=cr)
>
> I get the following error:
>
> ---------------------------------------------------------------------------
> IndexError                                Traceback (most recent call last)
> <ipython-input-14-4d67ba49bdf1> in <module>()
> ----> 1 proj = ProjectionPlot(pf, "z",
> ('Gas','Density'),weight_field=('Gas','Density'),data_source=cr)
>
> /1/home/benincsm/yt/yt-3.0-test/yt/yt/visualization/plot_window.pyc in
> __init__(self, pf, axis, fields, center, width, axes_unit, weight_field,
> max_level, origin, fontsize, field_parameters, data_source, proj_style,
> window_size, aspect)
>    1217         proj = pf.proj(fields, axis, weight_field=weight_field,
>    1218                        center=center, data_source=data_source,
> -> 1219                        field_parameters = field_parameters, style
> = proj_style)
>    1220         PWViewerMPL.__init__(self, proj, bounds, fields=fields,
> origin=origin,
>    1221                              fontsize=fontsize,
> window_size=window_size, aspect=aspect)
>
> /1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/construction_data_containers.pyc
> in __init__(self, field, axis, weight_field, center, pf, data_source,
> style, field_parameters)
>     226         self.data_source = data_source
>     227         self.weight_field = weight_field
> --> 228         self.get_data(field)
>     229
>     230     @property
>
> /1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/construction_data_containers.pyc
> in get_data(self, fields)
>     266         # We do this once
>     267         for chunk in self.data_source.chunks([], "io", local_only
> = False):
> --> 268             self._initialize_chunk(chunk, tree)
>     269         # This needs to be parallel_objects-ified
>     270         with
> self.data_source._field_parameter_state(self.field_parameters):
>
> /1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/construction_data_containers.pyc
> in _initialize_chunk(self, chunk, tree)
>     342
>     343     def _initialize_chunk(self, chunk, tree):
> --> 344         icoords = chunk.icoords
>     345         xax = self.pf.coordinates.x_axis[self.axis]
>     346         yax = self.pf.coordinates.y_axis[self.axis]
>
> /1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/selection_data_containers.pyc
> in icoords(self)
>     733     @property
>     734     def icoords(self):
> --> 735         return self.base_object.icoords[self._cond_ind,:]
>     736
>     737     @property
>
> IndexError: index 9025361 is out of bounds for axis 0 with size 9025360
>
> Is the imaging step not yet supported in yt-3.0?
>
> Thanks,
> Sam
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>


-- 
Cameron Hummels
Postdoctoral Researcher
Steward Observatory
University of Arizona
http://chummels.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140716/14485841/attachment.html>


More information about the yt-users mailing list