[Users] 3D Carpet component with one grid point along the z direction

[Gabriele Bozzola]

Hello,

one of my evolutions produced the attached file.

What is interesting about this file is that refinement level 10 has a
component
that has shape (1, 53, 105). So, we have only one point along the z
direction.
I am not entirely sure how this came to be, but the problem is that the HDF5
attribute 'delta' reports a cell width of 0 along the z direction:

In [12]: with h5py.File("alp.xyz.file_0.h5") as f:
    ...:         print(f['ADMBASE::alp it=55296 tl=0 rl=10 c=0'].shape)
(1, 53, 105)

In [13]: with h5py.File("alp.xyz.file_0.h5") as f:
    ...:         print(f['ADMBASE::alp it=55296 tl=0 rl=10
c=0'].attrs['delta'])
[0.072 0.072 0.   ]

I would have expected the value of 0.072, even if the cell has only one
grid point.

Is this a bug?

Thanks,
Gabriele
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20210412/a2e2a189/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: alp.xyz.file_0.h5
Type: application/x-hdf
Size: 1664262 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20210412/a2e2a189/attachment-0001.hdf