<div dir="ltr">Hello,<div><br></div><div>one of my evolutions produced the attached file.</div><div><br></div><div>What is interesting about this file is that refinement level 10 has a component</div><div>that has shape (1, 53, 105). So, we have only one point along the z direction.</div><div>I am not entirely sure how this came to be, but the problem is that the HDF5<br></div><div>attribute &#39;delta&#39; reports a cell width of 0 along the z direction:</div><div><br></div><div>In [12]: with h5py.File(&quot;alp.xyz.file_0.h5&quot;) as f:<br>    ...:         print(f[&#39;ADMBASE::alp it=55296 tl=0 rl=10 c=0&#39;].shape)<br>(1, 53, 105)<br><br>In [13]: with h5py.File(&quot;alp.xyz.file_0.h5&quot;) as f:<br>    ...:         print(f[&#39;ADMBASE::alp it=55296 tl=0 rl=10 c=0&#39;].attrs[&#39;delta&#39;])<br>[0.072 0.072 0.   ]<br></div><div><br></div><div>I would have expected the value of 0.072, even if the cell has only one grid point.</div><div><br></div><div>Is this a bug?</div><div><br></div><div>Thanks,</div><div>Gabriele</div></div>