<div dir="ltr"><div><div><div><div>Hi everyone,<br><br>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:<br>
<br>cr = ad.cut_region(["obj['Gas','Temperature'] < 1000"])<br><br></div>However, when I try to image that region:<br><br>proj = ProjectionPlot(pf, "z", ('Gas','Density'),weight_field=('Gas','Density'),data_source=cr)<br>

<br></div>I get the following error:<br><br>---------------------------------------------------------------------------<br>IndexError                                Traceback (most recent call last)<br><ipython-input-14-4d67ba49bdf1> in <module>()<br>

----> 1 proj = ProjectionPlot(pf, "z", ('Gas','Density'),weight_field=('Gas','Density'),data_source=cr)<br><br>/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)<br>
   1217         proj = pf.proj(fields, axis, weight_field=weight_field,<br>   1218                        center=center, data_source=data_source,<br>-> 1219                        field_parameters = field_parameters, style = proj_style)<br>

   1220         PWViewerMPL.__init__(self, proj, bounds, fields=fields, origin=origin,<br>   1221                              fontsize=fontsize, window_size=window_size, aspect=aspect)<br><br>/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)<br>

    226         self.data_source = data_source<br>    227         self.weight_field = weight_field<br>--> 228         self.get_data(field)<br>    229 <br>    230     @property<br><br>/1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/construction_data_containers.pyc in get_data(self, fields)<br>

    266         # We do this once<br>    267         for chunk in self.data_source.chunks([], "io", local_only = False):<br>--> 268             self._initialize_chunk(chunk, tree)<br>    269         # This needs to be parallel_objects-ified<br>

    270         with self.data_source._field_parameter_state(self.field_parameters):<br><br>/1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/construction_data_containers.pyc in _initialize_chunk(self, chunk, tree)<br>    342 <br>

    343     def _initialize_chunk(self, chunk, tree):<br>--> 344         icoords = chunk.icoords<br>    345         xax = self.pf.coordinates.x_axis[self.axis]<br>    346         yax = self.pf.coordinates.y_axis[self.axis]<br>

<br>/1/home/benincsm/yt/yt-3.0-test/yt/yt/data_objects/selection_data_containers.pyc in icoords(self)<br>    733     @property<br>    734     def icoords(self):<br>--> 735         return self.base_object.icoords[self._cond_ind,:]<br>

    736 <br>    737     @property<br><br>IndexError: index 9025361 is out of bounds for axis 0 with size 9025360<br><br></div>Is the imaging step not yet supported in yt-3.0?<br><br></div>Thanks,<br>Sam</div>