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

[Matthew Turk]

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