[yt-users] cut_region index error with Tipsy

[Samantha Benincasa]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140716/bab30fa1/attachment.htm>