[yt-dev] Issue #610: Trouble plotting a simple ramses dataset (yt_analysis/yt)

[Nathan Goldbaum]

I've got a number of 2D datasets if you'd like to take a look.  See also
the comment I made on the issue about a 3D sedov-taylor problem that fails
with a different traceback.

I'm going to try debugging, but I've never touched this code before so I
may be useless ;)


On Mon, Jul 22, 2013 at 2:10 PM, Matthew Turk <matthewturk at gmail.com> wrote:

> Hi Nathan,
>
> I haven't tested 2D, but I have not ever worked with it, either. This
> error has also in the past shown up when the field counts did not match
> what was in the file.
>
> Matt
> On Jul 22, 2013 5:02 PM, "Nathan Goldbaum" <issues-reply at bitbucket.org>
> wrote:
>
>> New issue 610: Trouble plotting a simple ramses dataset
>>
>> https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ramses-dataset
>>
>> Nathan Goldbaum:
>>
>> Using the following script:
>>
>>
>> ```
>> #!python
>>
>> from yt.mods import *
>> pf = load('/pfs/hennebel/KH/output_00037/info_00037.txt')
>> slc = SlicePlot(pf, 'z', 'Density')
>> ```
>>
>> I get the following traceback:
>>
>>
>> ```
>> #!python traceback
>>
>> Traceback (most recent call last):
>>   File "test.py", line 3, in <module>
>>     slc = SlicePlot(pf, 'z', 'Density')
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/visualization/plot_window.py", line
>> 1254, in __init__
>>     slc = pf.h.slice(axis, center[axis],
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/data_objects/static_output.py",
>> line 218, in hierarchy
>>     self, data_style=self.data_style)
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py",
>> line 303, in __init__
>>     super(RAMSESGeometryHandler, self).__init__(pf, data_style)
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/geometry/geometry_handler.py", line
>> 68, in __init__
>>     self._setup_geometry()
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/geometry/oct_geometry_handler.py",
>> line 50, in _setup_geometry
>>     self._initialize_oct_handler()
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py",
>> line 308, in _initialize_oct_handler
>>     for i in range(self.parameter_file['ncpu'])]
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py",
>> line 72, in __init__
>>     self._read_amr()
>>   File
>> "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py",
>> line 250, in _read_amr
>>     assert(n == ng)
>> AssertionError
>> ```
>>
>> The dataset is a simple 3D kelvin helmholtz test.  I've uploaded a copy
>> to my website (warning, 5GB): http://ucolick.org/~goldbaum/files/kh.tar
>>
>> I don't know much about this dataset as I got it precanned from Patrick
>> Hennebelle so sorry if I can't provide much detail.  FWIW, this is on the
>> current tip of 3.0 (250ee66bf017).
>>
>>
>> _______________________________________________
>> yt-dev mailing list
>> yt-dev at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>>
>
> _______________________________________________
> yt-dev mailing list
> yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20130722/ed8a80c7/attachment.html>