[yt-users] Index error while using covering_grid

[Yang Haifeng]

Sorry, forgot the output:

level    # grids           # cells         # cells^3
----------------------------------------------
  0         1            262144                64
  1         1             24576                30
  2         1             38720                34
  3         1             98304                47
  4         1            147456                53
  5         1            327680                69
  6         1            269120                65
  7         1            403200                74
  8         1            645504                87
----------------------------------------------
            9           2216704


t = 1.06267789e-01 = 1.14972625e+12 s = 3.64326264e+04 years

Smallest Cell:
    Width: 9.890e-12 Mpc
    Width: 9.890e-06 pc
    Width: 2.040e+00 AU
    Width: 3.052e+13 cm


2015-01-17 15:16 GMT-05:00 Yang Haifeng <hfyangpku at gmail.com>:

> Hi Nathan,
>
> Here is the output of ds.print_stats(). I think it hasn't run off the edge
> yet. The refine ratio is 2 and top grid is 64*3. ad_c["Density"] will not
> produce any error.
>
> 2015-01-17 12:48 GMT-05:00 Nathan Goldbaum <nathan12343 at gmail.com>:
>
> Hi Haifeng,
>>
>> On Sat, Jan 17, 2015 at 8:55 AM, Yang Haifeng <hy4px at virginia.edu> wrote:
>>
>>> Dear yt users,
>>>
>>> I am trying to make covering_grid from an enzo simulation. The code is
>>> very simple: (somewhat sloppy since it's running with ipython. )
>>>
>>> ds=yt.load("data0201")
>>>
>>
>> What is the output of "ds.print_stats()" for this dataset?  I'm wondering
>> if your covering grid is running off the edge of the simulation domain.
>>
>>
>>> ad_c=ds.covering_grid(level=3, left_edge=0.4375, dims=[64,64,64])
>>> ad_c["x"]
>>>
>>> where the ad_c["x"] will give an Index Error on my desktop. However, the
>>> same piece of code runs great in a different machine. The version of yt are
>>> the same in both machines, which is 3.0.2. Both of them installed yt
>>> through anaconda. The error message I got was following. Thanks for your
>>> help!
>>>
>>>
>>> ---------------------------------------------------------------------------
>>> IndexError                                Traceback (most recent call
>>> last)
>>> <ipython-input-6-ff78a8bccc54> in <module>()
>>> ----> 1 ad_c["x"]
>>>
>>> /home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/data_objects/data_containers.pyc
>>> in __getitem__(self, key)
>>>     243             if f in self._container_fields:
>>>     244                 self.field_data[f] = \
>>> --> 245
>>> self.ds.arr(self._generate_container_field(f))
>>>     246                 return self.field_data[f]
>>>     247             else:
>>>
>>> /home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/data_objects/construction_data_containers.pyc
>>> in _generate_container_field(self, field)
>>>     587             np.multiply(rv, self.dds[2], rv)
>>>     588         elif field == ("index", "x"):
>>> --> 589             x = np.mgrid[self.left_edge[0] + 0.5*self.dds[0]:
>>>     590                          self.right_edge[0] - 0.5*self.dds[0]:
>>>     591                          self.ActiveDimensions[0] * 1j]
>>>
>>> /home/fxyang/software/anaconda/lib/python2.7/site-packages/yt/units/yt_array.pyc
>>> in __getitem__(self, item)
>>>     981
>>>     982     def __getitem__(self, item):
>>> --> 983         ret = super(YTArray, self).__getitem__(item)
>>>     984         if ret.shape == ():
>>>     985             return YTQuantity(ret, self.units)
>>>
>>> IndexError: too many indices for array
>>>
>>>
>>> --
>>> Best,
>>> Haifeng
>>>
>>>
>>> _______________________________________________
>>> 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
>>
>>
>
>
> --
> Best,
> Haifeng
>



-- 
Best,
Haifeng
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150117/4c909d01/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org