Re: [yt-users] covering_grid with RAMSES dataset

Hi Nathan, thanks for the feedback. I will create the issue right away. Enrico On Wed, 23 Mar 2016 09:47:16 -0700 Nathan Goldbaum wrote:

Hi Enrico,

I don't think you're doing anything wrong - this is a bug, likely a generic issue with the way covering grids get filled out for octree data.

This issue with smoothed_covering_grid is likely another, different but related bug.

Can you create an issue to track this?

https://bitbucket.org/yt_analysis/yt/issues/new

-Nathan

On Wed, Mar 23, 2016 at 9:36 AM, Enrico Garaldi wrote:

...

Hi, I am having troubles using the smoothed_covering_grid and covering_grid with RAMSES dataset. I have yt-project version 3.2.3 and I am using the "output_00080" dataset from the yt sample datasets.

smoothed_covering_grid fails (raising a RuntimeError because tot != 0, line 880 of construction_data_containers.pyc), while covering_grid seems to leave unfilled the refined regions. What I mean is that if a cell is refined beyond the maximum allowed level for the grid, it is left unfilled (contrary to what happens with ENZO datasets).

I am using this code:

import yt import matplotlib.pyplot as plt

ds = yt.load('/users/egaraldi/Documents/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(density_field[:,:,48].v, interpolation='none') plt.show() #holes with zero-density corresponding to refined regions

Am I doing something wrong?

Thanks in advance, Enrico

------------------------------------- Enrico Garaldi - Ph.D. student Argelander-Institut für Astronomie Address: Auf dem Hügel 71, 53121 Bonn, Germany Room: 1.024 Phone: +49 (0) 228 73 3433 egaraldi[at]uni-bonn[dot]de _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

------------------------------------- Enrico Garaldi - Ph.D. student Argelander-Institut für Astronomie Address: Auf dem Hügel 71, 53121 Bonn, Germany Room: 1.024 Phone: +49 (0) 228 73 3433 egaraldi[at]uni-bonn[dot]de