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