[yt-users] covering_grid with RAMSES dataset
Enrico Garaldi
egaraldi at uni-bonn.de
Wed Mar 23 10:24:07 PDT 2016
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 <nathan12343 at gmail.com> 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
><egaraldi at uni-bonn.de>
> 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 at 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
More information about the yt-users
mailing list