[yt-users] cut_region index error with Tipsy

[Cameron Hummels]

John is correct.  I just tested this with an enzo dataset, and it works
fine.  In fact, I just tested it with a sample Gasoline dataset, and I got
cut_regions working OK in ProjectionPlots.


On Wed, Jul 16, 2014 at 1:06 PM, John ZuHone <jzuhone at gmail.com> wrote:

> I've successfully used cut regions in projections before, so this isn't
> it, I don't think.
>
> On Jul 16, 2014, at 4:04 PM, Cameron Hummels <chummels at gmail.com> wrote:
>
> 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
>  _______________________________________________
> yt-users mailing list
> 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
>
>


-- 
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/6cc2e42c/attachment.html>