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