I am getting the following error while trying to read a file I generated with numpy and h5py.
This is the structure of the dataset: it is 4d array
I am calling these within a Fortran subroutine, below are the calls I use:
rank = 4
count = (/ 1, 1, 1, 1 /)
offset = (/ 0, 0, 0, dr_meshMe /)
stride = (/ 1, 1, 1, 1 /)
dsetname = “/u”
data_dims = (/ nxb, nyb, nzb, 1 /)
block1 = (/ nxb, nyb, nzb, 1 /)
allocate(sdata(nxb, nyb, nzb, 1))
! ! read file
call h5dopen_f(file_id, dsetname, dset_id, error)
call h5dget_space_f(dset_id, dataspace, error)
call h5sselect_hyperslab_f(dataspace, H5S_SELECT_SET_F, offset, count, error, stride, block1)
call h5screate_simple_f(rank, data_dims, memspace, error)
call h5dread_f(dset_id, H5T_IEEE_F64LE, sdata, data_dims, error, memspace, dataspace)
call h5dclose_f(dset_id, error)
Thanks in advance for the help!