sll_k_hdf5_par_write_array.F90

Go to the documentation of this file.

integer(hid_t)   :: plist_id
integer(hid_t)   :: dset_id
integer(hid_t)   :: memspace
integer(hid_t)   :: filespace
integer(hsize_t) :: dimsfi(dspace_dims)
integer(hsize_t) :: block(dspace_dims)
integer(hsize_t) :: start(dspace_dims)
integer(hsize_t) :: count(dspace_dims)
integer(hsize_t) :: stride(dspace_dims)
integer          :: rank
 
! Basic dataset parameters
rank = dspace_dims
count = 1
stride = 1
 
! Integer kind type conversions
dimsfi = int(global_size, hsize_t)
block = int(shape(array), hsize_t)
start = int(offset, hsize_t)
 
!
! Create dataspaces ('file space' and 'memory space').
!
call h5screate_simple_f(rank, dimsfi, filespace, error)
sll_assert_always(error == 0)
 
call h5screate_simple_f(rank, block, memspace, error)
sll_assert_always(error == 0)
 
!
! Create dataset.
!
call h5pcreate_f(h5p_dataset_create_f, plist_id, error)
sll_assert_always(error == 0)
 
!
! If required, use chunked layout
!
if (present(chunk_dims)) then
   call h5pset_chunk_f(plist_id, rank, chunk_dims, error)
   sll_assert_always(error == 0)
end if
 
call h5dcreate_f(handle%file_id, dsetname, datatype, filespace, &
                 dset_id, error, plist_id)
sll_assert_always(error == 0)
 
call h5pclose_f(plist_id, error)
sll_assert_always(error == 0)
 
call h5sclose_f(filespace, error)
sll_assert_always(error == 0)
 
!
! Each process defines dataset in memory and writes it to the hyperslab
! in the file.
!
! Select hyperslab in the file.
!
call h5dget_space_f(dset_id, filespace, error)
sll_assert_always(error == 0)
 
call h5sselect_hyperslab_f(filespace, h5s_select_set_f, &
                           start, count, error, stride, block)
sll_assert_always(error == 0)
 
!
! Initialize data buffer with trivial data.
!
!
! Create property list for collective dataset write
!
call h5pcreate_f(h5p_dataset_xfer_f, plist_id, error)
sll_assert_always(error == 0)
 
call h5pset_dxpl_mpio_f(plist_id, h5fd_mpio_collective_f, error)
sll_assert_always(error == 0)
 
!
! Write the dataset collectively.
!
call h5dwrite_f(dset_id, datatype, array, dimsfi, error, &
                file_space_id=filespace, mem_space_id=memspace, &
                xfer_prp=plist_id)
sll_assert_always(error == 0)
 
!
! Close the property list.
!
call h5pclose_f(plist_id, error)
sll_assert_always(error == 0)
 
!
! Close dataspaces.
!
call h5sclose_f(filespace, error)
sll_assert_always(error == 0)
 
call h5sclose_f(memspace, error)
sll_assert_always(error == 0)
 
!
! Close the dataset.
!
call h5dclose_f(dset_id, error)
sll_assert_always(error == 0)