Hi, I think Nathan's approach is the right one -- more specifically because `max_level` is not a property of the values inside a data container, but rather an instruction to the selection machinery not to select above that level. On Sun, Sep 8, 2019 at 2:17 PM Nathan <nathan.goldbaum@gmail.com> wrote:
Maybe you just want to find the extrema of the (“index” “grid_level”) field? That field gives the AMR level of every leaf zone in the data object.
On Sun, Sep 8, 2019 at 1:02 PM Guido granda muñoz <guidogranda@gmail.com> wrote:
Hello Yt users, I'm trying to obtain the maximum and minimum level of refinement of a data container in a simulation. What I am doing is the following:
ds=yt.load(bname) sp_open=ds.disk(center,normal,radius,(res,'pc')) # smallest disk possible
Where bname is the name of the snapshot, center, normal, radius, and res, are the desired center, normal, radius and height of the disk. I chose the resolution of the simulation as the height of the disk because I want it to represent a disk-like differential area and I the disk is located at the region of highest level of refinement.
I tried to obtain the minimum and maximum levels of refinement over the sp_open subregion, however only max_level is implemented doing:
sp_open.index.max_level
Do you know how I can obtain the minimum level of refinement of a data container of a sub-region?
Besides, when I print the level statistics of the subregion and the ones of the whole computational volume I get the same output. I print them using:
sp_open.index.level_stats # subregion ds.index.level_stats # computational volume And I get: rec.array([( 1, 512, 0), ( 8, 4096, 1), ( 64, 32768, 2), ( 512, 262144, 3), ( 4096, 2097152, 4), (32768, 16777216, 5), ( 0, 0, 6), ( 0, 0, 7), ( 0, 0, 8), ( 0, 0, 9), ( 0, 0, 10), ( 0, 0, 11), ( 0, 0, 12), ( 0, 0, 13), ( 0, 0, 14), ( 0, 0, 15), ( 0, 0, 16), ( 0, 0, 17), ( 0, 0, 18), ( 0, 0, 19), ( 0, 0, 20), ( 0, 0, 21), ( 0, 0, 22), ( 0, 0, 23), ( 0, 0, 24), ( 0, 0, 25), ( 0, 0, 26), ( 0, 0, 27), ( 0, 0, 28), ( 0, 0, 29), ( 0, 0, 30), ( 0, 0, 31), ( 0, 0, 32), ( 0, 0, 33), ( 0, 0, 34), ( 0, 0, 35), ( 0, 0, 36), ( 0, 0, 37), ( 0, 0, 38), ( 0, 0, 39), ( 0, 0, 40), ( 0, 0, 41), ( 0, 0, 42), ( 0, 0, 43), ( 0, 0, 44), ( 0, 0, 45), ( 0, 0, 46), ( 0, 0, 47)], dtype=[('numgrids', '<i8'), ('numcells', '<i8'), ('level', '<i8')])
Which does not make sense to me, because the sub-region level statistics should not include coarse grids, even though this sub-region is in fact inside lower refinement regions. Shouldn't this statistics just print the lowest and highest refinement level inside the sub-region ?
Thank you, -- Guido _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org