
<html><body>

<p>New issue 1191: covering_grid with RAMSES dataset <a href="http://link.bitbucket.org/wf/click?upn=8USRlNyft-2BCzk2l4Ywl6gDx2lD2xxoS9E7MwXb2SMR-2FnGLXjyWbHcdK6NXtIgqd-2FEL171OetN9gKh1z12lSUz-2Bvhdb8R7mb1nh5YEmkxpXqSLVlnWVhtD1QnVsUM5D7Z_JlBWjyZTyKwhM1rBKe-2BcAdFPwdECN4PBooLGVSVl7EKJU04db3QpHMYWBfyuXRxv-2F-2FCx4YQVNGRVJg0iEHZIOieyquKnVZ9TisA4xNkW9dKpoAk0T-2BwOIZdiOi51Ub3G-2FDaNvoAsBu1KdZihCbPx-2FkyU9tfAc0Czo5MFRMgTIsRTpyZ-2BizZpDPPWLmBLqecwReQFPiYGFeBnAbhWOHQmTcQ3w6xd6nIiOUjlHpxO7YQ-3D">https://bitbucket.org/yt_analysis/yt/issues/1191/covering_grid-with-ramses-dataset</a></p>

<p>Enrico Garaldi:</p>

<p>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.</p>

<p>You can verify this using:</p>

<p>``` #!python import yt import matplotlib.pyplot as plt</p>

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

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

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

<p>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](<a href="http://link.bitbucket.org/wf/click?upn=8USRlNyft-2BCzk2l4Ywl6gJkmhpCq670QTqJWWjYIuPnNm0cUwoi12sahp1TzmuVyf0s8fYzjWpRcEWx8ACFP4pdhzZcXiTJX3IE8AzaQboc-3D_JlBWjyZTyKwhM1rBKe-2BcAdFPwdECN4PBooLGVSVl7EKJU04db3QpHMYWBfyuXRxv-2F-2FCx4YQVNGRVJg0iEHZIOmvujXZLlXTi9GcpAVp1cbZlBcX9nwqR-2Br23CFQ1c9fZocRfg-2B248yiE-2BCRVl71gmfRjX7vnSA56h9O8bMdqdNnleBRLsBNUUcxpn1mbLikM5MnHiH-2F9mf69QnASK-2F9g-2Bp6v31Ay0RoV3b5nWfDJJl4-3D">https://bitbucket.org/repo/BXbAb/images/2447880600-cgrid.png</a>)</p>

<p>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)</p>


<img src="http://link.bitbucket.org/wf/open?upn=JlBWjyZTyKwhM1rBKe-2BcAdFPwdECN4PBooLGVSVl7EKJU04db3QpHMYWBfyuXRxv-2F-2FCx4YQVNGRVJg0iEHZIOnwPmRN2tFeoJkYys2W-2FoXzOZ4xLdiB0vZbcItxrp5VefSVLX-2BYEE3RUh2kxpKk4grXwnfHIgb4BuIznxiwNHxw1l9vaRRuruTCfELeVh9xjFz2lx3pcSRaeYE36IA6Ak4nwcETW0FhgQXyoPmjq-2F8M-3D" alt="" width="1" height="1" border="0" style="height:1px !important;width:1px !important;border-width:0 !important;margin-top:0 !important;margin-bottom:0 !important;margin-right:0 !important;margin-left:0 !important;padding-top:0 !important;padding-bottom:0 !important;padding-right:0 !important;padding-left:0 !important;"/>

</body></html>