Issue #1191: covering_grid with RAMSES dataset (yt_analysis/yt)
New issue 1191: covering_grid with RAMSES dataset http://link.bitbucket.org/wf/click?upn=8USRlNyft-2BCzk2l4Ywl6gDx2lD2xxoS9E7M... 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-2BCzk2l4Ywl6gJkmhpCq670QTqJ... 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)
participants (1)
-
Enrico Garaldi