Issue #610: Trouble plotting a simple ramses dataset (yt_analysis/yt)
New issue 610: Trouble plotting a simple ramses dataset https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ram... 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).
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@bitbucket.org> wrote:
New issue 610: Trouble plotting a simple ramses dataset
https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ram...
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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
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@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@bitbucket.org> wrote:
New issue 610: Trouble plotting a simple ramses dataset
https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ram...
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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
_______________________________________________ yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
participants (3)
-
Matthew Turk -
Nathan Goldbaum -
Nathan Goldbaum