writereal1darrh5()

subroutine readwriteh5::writereal1darrh5	(	character(len=1024), intent(in)	outfile,
		character(len=*), intent(in)	dsetname,
		real(dp), dimension(:), intent(in)	arr
	)

Definition at line 55 of file rdwrhdf5.f08.

     
    character(len=1024),intent(in) :: outfile
    real(dp),intent(in) :: arr(:)
    character(len=*),intent(in) :: dsetname
    
    ! filename length must be the same than in dummy variable 
    integer(hid_t)     :: file_id      ! file identifier
    integer(hid_t)     :: dset_id      ! dataset identifier
    integer            :: error        ! error flag
    integer(hid_t)     :: dspace_id
    integer(hsize_t), dimension(1) :: data_dims1
    integer :: rank
       
    !! copmpression stuff
    logical :: avail
    integer ::  filter_info
    integer :: filter_info_both
    integer(hid_t)  ::  dcpl
    integer  :: chunkingfactor
    integer(hsize_t), dimension(1) :: chunk1

    
    chunkingfactor = 10

    !##############################
    !#    initialize hdf5         #
    !##############################
    ! initialize fortran interface.
    call h5open_f(error)
      !-------------------------------------------------------------
    !  check if gzip compression is available and can be used for both
    !  compression and decompression.  normally we do not perform error
    !  checking in these examples for the sake of clarity, but in this
    !  case we will make an exception because this filter is an
    !  optional part of the hdf5 library.
    call h5zfilter_avail_f(h5z_filter_deflate_f, avail, error)
    if (.not.avail) then
       write(*,'("gzip filter not available.",/)')
       stop
    endif
    call h5zget_filter_info_f(h5z_filter_deflate_f, filter_info, error)
    filter_info_both=ior(h5z_filter_encode_enabled_f,h5z_filter_decode_enabled_f)
    if (filter_info .ne. filter_info_both) then
       write(*,'("gzip filter not available for encoding and decoding.",/)')
       stop
    endif

    !-------------------------------------------------------------
    if (h5firsttimeread.eqv..true.) then
       call h5fcreate_f(trim(adjustl(outfile)), h5f_acc_trunc_f,file_id,error)
       h5firsttimeread=.false.
    else
       ! open an existing file using the default properties.
       call h5fopen_f(trim(adjustl(outfile)), h5f_acc_rdwr_f, file_id, error)
    end if

    
    !##############################
    !#  write arr                 #
    !##############################
    rank=1
    data_dims1=shape(arr)
    chunk1 = data_dims1/chunkingfactor
    call checkchunking(chunk1)
    ! create dataspace.  setting maximum size to null sets the maximum
    ! size to be the current size.
    call h5screate_simple_f (rank, data_dims1, dspace_id, error)
    !-------------------------------------------------------------
    ! create the dataset creation property list, add the gzip
    ! compression filter and set the chunk size.
    call h5pcreate_f(h5p_dataset_create_f, dcpl, error)
    call h5pset_deflate_f(dcpl, 9, error)
    call h5pset_chunk_f(dcpl, rank, chunk1, error)
    !-------------------------------------------------------------
    ! create the dataset with default properties.
    call h5dcreate_f(file_id, dsetname, h5t_ieee_f64le, dspace_id, &
         dset_id, error, dcpl)
    !-------------------------------------------------------------
    ! write the data to the dataset.
    call h5dwrite_f(dset_id,h5t_ieee_f64le, arr, data_dims1, error)
    !-------------------------------------------------------------
    ! close and release resources.
    call h5pclose_f(dcpl,error)
    !-------------------------------------------------------------
    call h5sclose_f(dspace_id, error)
    ! close the dataset.
    call h5dclose_f(dset_id, error)


    ! ##############################
    ! #    terminate hdf5 stuff    #
    ! ##############################
    ! close the file.
    call h5fclose_f(file_id, error)
    ! close fortran interface.
    call h5close_f(error)
    
    return

Here is the call graph for this function:

Here is the caller graph for this function: