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

Nathan Goldbaum issues-reply at bitbucket.org
Mon Jul 22 14:02:55 PDT 2013


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).





More information about the yt-dev mailing list