[yt-dev] Issue #959: Ray provides non-contiguous cell information (yt_analysis/yt)

[Nicholas Earl]

New issue 959: Ray provides non-contiguous cell information
https://bitbucket.org/yt_analysis/yt/issue/959/ray-provides-non-contiguous-cell

Nicholas Earl:

Using the [Ray object](http://yt-project.org/doc/reference/api/generated/yt.data_objects.selection_data_containers.YTRayBase.html?highlight=ray#yt.data_objects.selection_data_containers.YTRayBase) 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](https://www.dropbox.com/s/xhjjqntc64h0lfk/DD0001_Slice_z_H_number_density.png?dl=0) with refinement overlaid.

You can see in the [Ray's hydrogen number density](https://www.dropbox.com/s/a3dp3ip7wjmnbh7/h_num_density_ray.png?dl=0) 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](http://paste.yt-project.org/show/5269/) 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
yt-dev at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org