correct very-long missing feature of reading with stride

Co-authored-by: jvdp1 <jvdp1@users.noreply.github.com>

Author: scivision

src/utils.f90

@@ -332,10 +332,12 @@
 istride = 1
 if(present(stride)) istride = int(stride, HSIZE_T)
 
-!! NOTE: 0-based hyperslab vs. 1-based Fortran
+!! NOTE: indices of 0-based hyperslab vs. 1-based Fortran indices
 i0 = istart - 1
 c_mem_dims = iend - i0
 
+if (any(istride > 1)) c_mem_dims = merge(c_mem_dims/istride, c_mem_dims/istride + 1, mod(c_mem_dims, istride)==0)
+
 if(size(mem_dims) == 0) then
   !! rank(dims(0)) == 1, but size(dims(0)) == 0
   if (sum(c_mem_dims) /= 1) error stop "ERROR:h5fortran:hdf_get_slice: scalar index of array failed " // dset_name

test/test_array.f90

@@ -10,7 +10,7 @@ program array_test
 print *, 'PASSED: array write'
 
 call test_read_slice('test_array.h5')
-print *, 'PASSED: slice read'
+print *, 'PASSED: slice read with or without stride'
 
 call test_write_slice('test_array.h5')
 print *, 'PASSED: slice write'