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.


On Jul 22, 2013 5:02 PM, "Nathan Goldbaum" <issues-reply@bitbucket.org> wrote:
New issue 610: Trouble plotting a simple ramses dataset

Nathan Goldbaum:

Using the following script:


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__
  File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/geometry/oct_geometry_handler.py", line 50, in _setup_geometry
  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__
  File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py", line 250, in _read_amr
    assert(n == ng)

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 mailing list