[yt-dev] Issue #1191: covering_grid with RAMSES dataset (yt_analysis/yt)

[Enrico Garaldi]

New issue 1191: covering_grid with RAMSES dataset
http://link.bitbucket.org/wf/click?upn=8USRlNyft-2BCzk2l4Ywl6gDx2lD2xxoS9E7MwXb2SMR-2FnGLXjyWbHcdK6NXtIgqd-2FEL171OetN9gKh1z12lSUz-2Bvhdb8R7mb1nh5YEmkxpXqSLVlnWVhtD1QnVsUM5D7Z_JlBWjyZTyKwhM1rBKe-2BcAdFPwdECN4PBooLGVSVl7EKJU04db3QpHMYWBfyuXRxv-2F-2FCx4YQVNGRVJg0iEHZIOpgdERFD1M7OAwBmSONFIJNoY3yBXK1PoHEaGZg5ouvOAcWGa-2Bhi-2BCaOGInl13cu3fk8-2BGSU12EBjk0joxAyU5zLnsdmi-2FLl3Bohr-2BYOZ66LGs-2FDoLwNeBhdU5oW0os4IB6pyNAFtbUptI6idxP-2BE1c-3D

Enrico Garaldi:

When using covering_grid with RAMSES datasets the cells refined beyond the level of the grid are left unfilled. This is not the case with e.g. ENZO datasets.

You can verify this using:

```
#!python
import yt
import matplotlib.pyplot as plt

ds = yt.load('output_00080/info_00080.txt')

cgrid = ds.covering_grid(0, left_edge=ds.domain_left_edge, dims=ds.domain_dimensions)
density_field = cgrid["density"]

print (density_field == 0.0).any()  #should be False but it is True

plt.imshow(numpy.log10(density_field[:,:,48].v), interpolation='none') #log10 to highlight the zeros in density_field
plt.show()    #holes with zero-density corresponding to refined regions
```
resulting in the following:
![cgrid.png](http://link.bitbucket.org/wf/click?upn=8USRlNyft-2BCzk2l4Ywl6gJkmhpCq670QTqJWWjYIuPnNm0cUwoi12sahp1TzmuVyf0s8fYzjWpRcEWx8ACFP4kqIeLBwXTiSEFf-2F1wW53Lc-3D_JlBWjyZTyKwhM1rBKe-2BcAdFPwdECN4PBooLGVSVl7EKJU04db3QpHMYWBfyuXRxv-2F-2FCx4YQVNGRVJg0iEHZIOsXvA6bEH21d1vKb9NSYiUKhfTJROyFt6z1yMXlAfJUBLNLYADd-2BUjCz-2BrsQh3OyJC0VbZGX2nOU5ZzanxK68VOXtalB0FbTa5EUpaB5sTUuZ-2F0qvxCMFieWJ9lmDlauY1WMjiuTE9MpJF5C4Zsbdms-3D

Another (related?) issue is that smoothed_covering_grids fails with RAMSES datasets (raising a RuntimeError because tot != 0, line 880 of construction_data_containers.pyc)


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20160323/d263e6fa/attachment.htm>