readreal2darrh5()

subroutine readwriteh5::readreal2darrh5	(	character(len=1024), intent(in)	inpfile,
		character(len=*), intent(in)	dsetname,
		real(dp), dimension(:,:), intent(inout), allocatable	arr
	)

Definition at line 322 of file rdwrhdf5.f08.

    implicit none
    character(len=1024),intent(in) :: inpfile
    real(dp),allocatable,intent(inout) :: arr(:,:)
    character(len=*),intent(in) :: dsetname
    
    ! filename length must be the same than in dummy variable 
    integer(hid_t)     :: file_id,dset_id   
    integer            :: rank,error       
    integer(hid_t)     :: dataspace_id
    integer(hsize_t),allocatable :: data_dims(:),maxdims(:)

    ! Initialize FORTRAN interface.
    CALL h5open_f(error)
    ! Open an existing file.
    !print*, 'Reading ',trim(inpfile)
    CALL h5fopen_f(trim(inpfile), h5f_acc_rdonly_f, file_id, error)
    if (error/=0) then
       write(*,*)
       write(*,*) "Error opening file from hdf5. "
       write(*,*) "filename: ",trim(inpfile)
       stop
    endif
  
    ! Open an existing dataset.
    call h5dopen_f(file_id, trim(dsetname), dset_id, error)
    if (error/=0) then
       write(*,*)
       write(*,*) "Error opening file from hdf5. "
       write(*,*) "err dsetname: ",trim(dsetname)
       stop
    endif
    ! get the dataspace id
    call h5dget_space_f(dset_id, dataspace_id, error) 
    ! get the rank of the dataset
    call h5sget_simple_extent_ndims_f(dataspace_id, rank, error)
    allocate(data_dims(rank),maxdims(rank))
    ! get the dimensions of the dataset
    call h5sget_simple_extent_dims_f(dataspace_id, data_dims, maxdims, error)
    !print*,'rank:',rank,'maxdims:',maxdims
    ! get the datatype
    !call h5dget_type_f(dset_id, datatype_id, error)
    allocate(arr(data_dims(1),data_dims(2)))
    print*,"Reading ",trim(dsetname), " with shape ", data_dims
    ! read data
    CALL h5dread_f(dset_id, h5t_ieee_f64le, arr, data_dims, error)
    ! Close the dataset.
    CALL h5dclose_f(dset_id, error)
    ! close the datatype
  
    ! Close the file.
    CALL h5fclose_f(file_id, error)
    ! Close FORTRAN interface.
    CALL h5close_f(error)

    return

Here is the caller graph for this function: