New issue 959: Ray provides non-contiguous cell information https://bitbucket.org/yt_analysis/yt/issue/959/ray-provides-non-contiguous-c...
Using the Ray object directly with enzo simulation data that uses AMR seems to give non-contiguous cell information. Mainly, it looks like the higher-resolution cells are appended to the end of the array.
For a simple example of a high density cloud in an ambient medium, you can see a slice of the hydrogen number density in the simulation here with refinement overlaid.
You can see in the Ray's hydrogen number density that the sharp peak around 80 corresponds with the density gradient at the edges of the cloud while the higher resolution information of the cloud interior has been tacked on to the end (i.e. the "top-hat" portion).
Printing of the Ray's fcoords also supports this. You can see that the x-coordinates of the cells between lines 167 and 200 are out of place.
yt-dev mailing list email@example.com http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org