<div dir="ltr">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.<div><br></div><div>I'm going to try debugging, but I've never touched this code before so I may be useless ;)</div>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jul 22, 2013 at 2:10 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">Hi Nathan,</p>
<p dir="ltr">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.</p>
<p dir="ltr">Matt</p><div class="HOEnZb"><div class="h5">
<div class="gmail_quote">On Jul 22, 2013 5:02 PM, "Nathan Goldbaum" <<a href="mailto:issues-reply@bitbucket.org" target="_blank">issues-reply@bitbucket.org</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
New issue 610: Trouble plotting a simple ramses dataset<br>
<a href="https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ramses-dataset" target="_blank">https://bitbucket.org/yt_analysis/yt/issue/610/trouble-plotting-a-simple-ramses-dataset</a><br>
<br>
Nathan Goldbaum:<br>
<br>
Using the following script:<br>
<br>
<br>
```<br>
#!python<br>
<br>
from yt.mods import *<br>
pf = load('/pfs/hennebel/KH/output_00037/info_00037.txt')<br>
slc = SlicePlot(pf, 'z', 'Density')<br>
```<br>
<br>
I get the following traceback:<br>
<br>
<br>
```<br>
#!python traceback<br>
<br>
Traceback (most recent call last):<br>
File "test.py", line 3, in <module><br>
slc = SlicePlot(pf, 'z', 'Density')<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/visualization/plot_window.py", line 1254, in __init__<br>
slc = pf.h.slice(axis, center[axis],<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/data_objects/static_output.py", line 218, in hierarchy<br>
self, data_style=self.data_style)<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py", line 303, in __init__<br>
super(RAMSESGeometryHandler, self).__init__(pf, data_style)<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/geometry/geometry_handler.py", line 68, in __init__<br>
self._setup_geometry()<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/geometry/oct_geometry_handler.py", line 50, in _setup_geometry<br>
self._initialize_oct_handler()<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py", line 308, in _initialize_oct_handler<br>
for i in range(self.parameter_file['ncpu'])]<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py", line 72, in __init__<br>
self._read_amr()<br>
File "/home/goldbaum/yt-x86_64/src/yt-3.0/yt/frontends/ramses/data_structures.py", line 250, in _read_amr<br>
assert(n == ng)<br>
AssertionError<br>
```<br>
<br>
The dataset is a simple 3D kelvin helmholtz test. I've uploaded a copy to my website (warning, 5GB): <a href="http://ucolick.org/~goldbaum/files/kh.tar" target="_blank">http://ucolick.org/~goldbaum/files/kh.tar</a><br>
<br>
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).<br>
<br>
<br>
_______________________________________________<br>
yt-dev mailing list<br>
<a href="mailto:yt-dev@lists.spacepope.org" target="_blank">yt-dev@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
</blockquote></div>
</div></div><br>_______________________________________________<br>
yt-dev mailing list<br>
<a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
<br></blockquote></div><br></div>