Initial commit from Polytechnique Montreal

1 files changed, 1410 insertions, 0 deletions

diff --git a/Ganlib/src/hdf5_wrap.f90 b/Ganlib/src/hdf5_wrap.f90
new file mode 100644
index 0000000..e1e7026
--- /dev/null
+++ b/Ganlib/src/hdf5_wrap.f90
@@ -0,0 +1,1410 @@
+!
+!-----------------------------------------------------------------------
+!
+!Purpose:
+! Fortran-2003 bindings for hdf5.
+!
+!Copyright:
+! Copyright (C) 2021 Ecole Polytechnique de Montreal
+! This library is free software; you can redistribute it and/or
+! modify it under the terms of the GNU Lesser General Public
+! License as published by the Free Software Foundation; either
+! version 2.1 of the License, or (at your option) any later version.
+!
+!Author(s): A. Hebert
+!
+!-----------------------------------------------------------------------
+!
+module hdf5_wrap
+   use, intrinsic :: iso_c_binding
+   private
+   integer, parameter :: MAX_NAME = 1024
+   public :: hdf5_open_file,hdf5_close_file,hdf5_get_dimensions,hdf5_get_shape, &
+             hdf5_list,hdf5_info,hdf5_read_data,hdf5_write_data,hdf5_list_datasets, &
+             hdf5_list_groups,hdf5_group_exists,hdf5_create_group,hdf5_delete, &
+             hdf5_copy
+
+   interface hdf5_read_data
+     module procedure hdf5_read_data_0d_int4,  hdf5_read_data_1d_int4,  &
+                      hdf5_read_data_2d_int4,  hdf5_read_data_3d_int4,  &
+                      hdf5_read_data_0d_real4, hdf5_read_data_1d_real4, &
+                      hdf5_read_data_2d_real4, hdf5_read_data_3d_real4, &
+                      hdf5_read_data_4d_real4, hdf5_read_data_0d_real8, &
+                      hdf5_read_data_1d_real8, hdf5_read_data_2d_real8, &
+                      hdf5_read_data_3d_real8, hdf5_read_data_4d_real8, &
+                      hdf5_read_data_0d_string,hdf5_read_data_1d_string
+   end interface hdf5_read_data
+   interface hdf5_write_data
+     module procedure hdf5_write_data_0d_int4,  hdf5_write_data_1d_int4,  &
+                      hdf5_write_data_2d_int4,  hdf5_write_data_3d_int4,  &
+                      hdf5_write_data_0d_real4, hdf5_write_data_1d_real4, &
+                      hdf5_write_data_2d_real4, hdf5_write_data_3d_real4, &
+                      hdf5_write_data_4d_real4, hdf5_write_data_0d_real8, &
+                      hdf5_write_data_1d_real8, hdf5_write_data_2d_real8, &
+                      hdf5_write_data_3d_real8, hdf5_write_data_4d_real8, &
+                      hdf5_write_data_0d_string,hdf5_write_data_1d_string
+   end interface hdf5_write_data
+
+   character(len=131) :: hsmg
+   character(kind=c_char), dimension(MAX_NAME) :: name1024
+
+contains
+subroutine STRCUT(name1, name2)
+   ! transform a Fortran string into a C null-terminated string
+   character(kind=c_char), dimension(*) :: name1
+   character(len=*) :: name2
+   integer :: ilong
+   interface
+      subroutine strcut_c (s, ct, n) bind(c)
+         use, intrinsic :: iso_c_binding
+         character(kind=c_char), dimension(*) :: s, ct
+         integer(c_int), value :: n
+      end subroutine strcut_c
+   end interface
+   ilong=len(trim(adjustl(name2)))
+   call strcut_c(name1, trim(adjustl(name2)), ilong)
+end subroutine STRCUT
+!
+subroutine STRFIL(name1, name2, nbyte)
+   ! transform a C null-terminated string into a Fortran string
+   character(len=*) :: name1
+   character(kind=c_char), dimension(*) :: name2
+   integer :: nbyte, ilong
+   interface
+      subroutine strfil_c (s, ct, n) bind(c)
+         use, intrinsic :: iso_c_binding
+         character(kind=c_char), dimension(*) :: s, ct
+         integer(c_int), value :: n
+      end subroutine strfil_c
+   end interface
+   ilong=len(name1)
+   name1=' '
+   call strfil_c(name1, name2, min(ilong,nbyte))
+end subroutine STRFIL
+!
+subroutine STRFIL1D(name1, name2, nbyte)
+   ! transform a C null-terminated string into a Fortran string array
+   character(len=*),dimension(:) :: name1
+   character(kind=c_char), dimension(*) :: name2
+   integer :: nbyte, ilong
+   character(len=MAX_NAME) :: text1024
+   interface
+      subroutine strfil_c (s, ct, n) bind(c)
+         use, intrinsic :: iso_c_binding
+         character(kind=c_char), dimension(*) :: s, ct
+         integer(c_int), value :: n
+      end subroutine strfil_c
+   end interface
+   ilong=len(name1)
+   idim=size(name1,1)
+   iof=1
+   do i=1,idim
+     name1(i)=' '
+     call strfil_c(text1024, name2(iof), nbyte)
+     name1(i)=text1024(:min(ilong,nbyte))
+     iof=iof+nbyte
+   enddo
+end subroutine STRFIL1D
+!
+subroutine hdf5_open_file(fname, ifile, rdonly)
+   ! open a HDF5 file
+   character(len=*), intent(in)  :: fname
+   type(c_ptr), intent(out) :: ifile
+   logical, optional :: rdonly
+   interface
+      subroutine hdf5_open_file_c (fname, ifile, irdonly) bind(c)
+         use, intrinsic :: iso_c_binding
+         character(kind=c_char), dimension(*) :: fname
+         type(c_ptr) :: ifile
+         integer(c_int), value :: irdonly
+      end subroutine hdf5_open_file_c
+   end interface
+   !
+   irdonly=0
+   if(present(rdonly)) then
+     if(rdonly) irdonly=1
+   endif
+   call STRCUT(name1024, fname)
+   call hdf5_open_file_c (name1024, ifile, irdonly)
+end subroutine hdf5_open_file
+!
+subroutine hdf5_close_file(ifile)
+   ! close a HDF5 file
+   type(c_ptr),intent(in) :: ifile
+   interface
+      subroutine hdf5_close_file_c (ifile) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+      end subroutine hdf5_close_file_c
+   end interface
+   call hdf5_close_file_c(ifile)
+end subroutine hdf5_close_file
+!
+subroutine hdf5_list(ifile, name)
+   ! table of contents
+   type(c_ptr) :: ifile
+   character(len=*),intent(in) :: name
+   !
+   interface
+      subroutine hdf5_list_c (ifile, namp) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+      end subroutine hdf5_list_c
+   end interface
+   call STRCUT(name1024, name)
+   call hdf5_list_c(ifile, name1024)
+end subroutine hdf5_list
+!
+subroutine hdf5_info(ifile, name, rank, type, nbyte, dimsr)
+   ! find dataset info
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer,intent(out) :: rank, type, nbyte
+   integer,target,dimension(5),intent(out) :: dimsr
+   !
+   type(c_ptr) :: pt_dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   dimsr(:5)=0
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+end subroutine hdf5_info
+!
+integer function hdf5_get_dimensions(ifile, name)
+   ! find dataset rank (number of dimensions
+   type(c_ptr),intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   integer :: rank
+   !
+   interface
+      subroutine hdf5_get_dimensions_c(ifile, namp, rank) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank
+      end subroutine hdf5_get_dimensions_c
+   end interface
+   call STRCUT(name1024, name)
+   call hdf5_get_dimensions_c(ifile, name1024, rank)
+   hdf5_get_dimensions = rank
+end function hdf5_get_dimensions
+!
+subroutine hdf5_get_shape(ifile, name, dimsr)
+   ! find dataset shape
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   integer, allocatable, target :: dimsr(:)
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr
+   !
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   rank=hdf5_get_dimensions(ifile, name)
+   allocate(dimsr(rank))
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+end subroutine hdf5_get_shape
+!
+subroutine hdf5_list_datasets(ifile, name, dsets)
+   ! collect daughter dataset names in a group
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   character(len=*), allocatable, dimension(:) :: dsets
+   integer :: nbobj, ndsets
+   !
+   character(kind=c_char), allocatable, dimension(:) :: pt_strim
+   interface
+      subroutine hdf5_get_num_group_c(ifile, namp, nbobj) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: nbobj
+      end subroutine hdf5_get_num_group_c
+   end interface
+   interface
+      subroutine hdf5_list_datasets_c(ifile, namp, ndsets, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp, idata
+         integer(c_int) :: ndsets
+      end subroutine hdf5_list_datasets_c
+   end interface
+   !
+   call STRCUT(name1024, name)
+   call hdf5_get_num_group_c (ifile, name1024, nbobj)
+   if (nbobj .eq. 0) return
+   allocate(pt_strim(nbobj*MAX_NAME))
+   call hdf5_list_datasets_c (ifile, name1024, ndsets, pt_strim)
+   allocate(dsets(ndsets))
+   call STRFIL1D(dsets,pt_strim,MAX_NAME)
+   deallocate(pt_strim)
+end subroutine hdf5_list_datasets
+!
+subroutine hdf5_list_groups(ifile, name, groups)
+   ! collect daughter group names in a group
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   character(len=*), allocatable, dimension(:) :: groups
+   integer :: nbobj, ngroups
+   !
+   character(kind=c_char), allocatable, dimension(:) :: pt_strim
+   interface
+      subroutine hdf5_get_num_group_c(ifile, namp, nbobj) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: nbobj
+      end subroutine hdf5_get_num_group_c
+   end interface
+   interface
+      subroutine hdf5_list_groups_c(ifile, namp, ngroups, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp, idata
+         integer(c_int) :: ngroups
+      end subroutine hdf5_list_groups_c
+   end interface
+   !
+   call STRCUT(name1024, name)
+   call hdf5_get_num_group_c (ifile, name1024, nbobj)
+   if (nbobj .eq. 0) return
+   allocate(pt_strim(nbobj*MAX_NAME))
+   call hdf5_list_groups_c (ifile, name1024, ngroups, pt_strim)
+   allocate(groups(ngroups))
+   call STRFIL1D(groups,pt_strim,MAX_NAME)
+   deallocate(pt_strim)
+end subroutine hdf5_list_groups
+!
+function hdf5_group_exists(ifile, name) result(lexist)
+   ! test for existence of a group
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   logical :: lexist
+   interface
+      function hdf5_group_exists_c (ifile, namp) bind(c)
+         use, intrinsic :: iso_c_binding
+         integer(c_int) :: hdf5_group_exists_c
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+      end function hdf5_group_exists_c
+   end interface
+   call STRCUT(name1024, name)
+   lexist = (hdf5_group_exists_c(ifile, name1024) /= 0)
+end function hdf5_group_exists
+!
+subroutine hdf5_create_group(ifile, name)
+   ! create a group
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   interface
+      subroutine hdf5_create_group_c (ifile, namp) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+      end subroutine hdf5_create_group_c
+   end interface
+   call STRCUT(name1024, name)
+   call hdf5_create_group_c(ifile, name1024)
+end subroutine hdf5_create_group
+!
+subroutine hdf5_delete(ifile, name)
+   ! delete a group or a dataset
+   type(c_ptr), intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   interface
+      subroutine hdf5_delete_c (ifile, namp) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+      end subroutine hdf5_delete_c
+   end interface
+   call STRCUT(name1024, name)
+   call hdf5_delete_c(ifile, name1024)
+end subroutine hdf5_delete
+!
+subroutine hdf5_copy(ifile_s, name_s, ifile_d, name_d)
+   ! copy a group or a dataset from one location to another
+   type(c_ptr), intent(in) :: ifile_s, ifile_d
+   character(len=*), intent(in) :: name_s, name_d
+   character(kind=c_char), dimension(MAX_NAME) :: name1024_d
+   interface
+      subroutine hdf5_copy_c (ifile_s, namp_s, ifile_d, namp_d) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile_s, ifile_d
+         character(kind=c_char), dimension(*) :: namp_s, namp_d
+      end subroutine hdf5_copy_c
+   end interface
+   call STRCUT(name1024, name_s)
+   call STRCUT(name1024_d, name_d)
+   call hdf5_copy_c(ifile_s, name1024, ifile_d, name1024_d)
+end subroutine hdf5_copy
+!
+subroutine hdf5_read_data_0d_int4(ifile, name, idata)
+   ! read a rank 0 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer, target :: idata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_int_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: idata
+      end subroutine hdf5_read_data_int_c
+   end interface
+   pt_data=c_loc(idata)
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if((rank.ne.1).or.(dimsr(1).ne.1)) then
+     write(hsmg,'(49hhdf5_read_data_0d_int4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   call hdf5_read_data_int_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_0d_int4
+!
+subroutine hdf5_read_data_1d_int4(ifile, name, idata)
+   ! read a rank 1 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer, allocatable, dimension(:), target :: idata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_int_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: idata
+      end subroutine hdf5_read_data_int_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.1) then
+     write(hsmg,'(49hhdf5_read_data_1d_int4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(idata(dimsr(1)))
+   pt_data=c_loc(idata)
+   call hdf5_read_data_int_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_1d_int4
+!
+subroutine hdf5_read_data_2d_int4(ifile, name, idata)
+   ! read a rank 2 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer, allocatable, dimension(:,:), target :: idata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_int_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: idata
+      end subroutine hdf5_read_data_int_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.2) then
+     write(hsmg,'(49hhdf5_read_data_2d_int4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(idata(dimsr(1),dimsr(2)))
+   pt_data=c_loc(idata)
+   call hdf5_read_data_int_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_2d_int4
+!
+subroutine hdf5_read_data_3d_int4(ifile, name, idata)
+   ! read a rank 3 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer, allocatable, dimension(:,:,:), target :: idata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_int_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: idata
+      end subroutine hdf5_read_data_int_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.3) then
+     write(hsmg,'(49hhdf5_read_data_3d_int4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(idata(dimsr(1),dimsr(2),dimsr(3)))
+   pt_data=c_loc(idata)
+   call hdf5_read_data_int_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_3d_int4
+!
+subroutine hdf5_read_data_0d_real4(ifile, name, rdata)
+   ! read a rank 0 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real4_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real4_c
+   end interface
+   pt_data=c_loc(rdata)
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if((rank.ne.1).or.(dimsr(1).ne.1)) then
+     write(hsmg,'(50hhdf5_read_data_0d_real4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   call hdf5_read_data_real4_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_0d_real4
+!
+subroutine hdf5_read_data_1d_real4(ifile, name, rdata)
+   ! read a rank 1 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4), allocatable, dimension(:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real4_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real4_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.1) then
+     write(hsmg,'(50hhdf5_read_data_1d_real4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real4_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_1d_real4
+!
+subroutine hdf5_read_data_2d_real4(ifile, name, rdata)
+   ! read a rank 2 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4), allocatable, dimension(:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real4_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real4_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.2) then
+     write(hsmg,'(50hhdf5_read_data_2d_real4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real4_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_2d_real4
+!
+subroutine hdf5_read_data_3d_real4(ifile, name, rdata)
+   ! read a rank 3 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4), allocatable, dimension(:,:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real4_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real4_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.3) then
+     write(hsmg,'(50hhdf5_read_data_3d_real4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2),dimsr(3)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real4_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_3d_real4
+!
+subroutine hdf5_read_data_4d_real4(ifile, name, rdata)
+   ! read a rank 4 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4), allocatable, dimension(:,:,:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real4_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real4_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.4) then
+     write(hsmg,'(50hhdf5_read_data_4d_real4: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2),dimsr(3),dimsr(4)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real4_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_4d_real4
+!
+subroutine hdf5_read_data_0d_real8(ifile, name, rdata)
+   ! read a rank 0 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real8_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real8_c
+   end interface
+   pt_data=c_loc(rdata)
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if((rank.ne.1).or.(dimsr(1).ne.1)) then
+     write(hsmg,'(50hhdf5_read_data_0d_real8: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   call hdf5_read_data_real8_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_0d_real8
+!
+subroutine hdf5_read_data_1d_real8(ifile, name, rdata)
+   ! read a rank 1 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8), allocatable, dimension(:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real8_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real8_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.1) then
+     write(hsmg,'(50hhdf5_read_data_1d_real8: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real8_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_1d_real8
+!
+subroutine hdf5_read_data_2d_real8(ifile, name, rdata)
+   ! read a rank 2 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8), allocatable, dimension(:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real8_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real8_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.2) then
+     write(hsmg,'(50hhdf5_read_data_2d_real8: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real8_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_2d_real8
+!
+subroutine hdf5_read_data_3d_real8(ifile, name, rdata)
+   ! read a rank 3 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8), allocatable, dimension(:,:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real8_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real8_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.3) then
+     write(hsmg,'(50hhdf5_read_data_3d_real8: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2),dimsr(3)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real8_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_3d_real8
+!
+subroutine hdf5_read_data_4d_real8(ifile, name, rdata)
+   ! read a rank 4 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8), allocatable, dimension(:,:,:,:), target :: rdata
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_real8_c(ifile, namp, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         type(c_ptr), value :: rdata
+      end subroutine hdf5_read_data_real8_c
+   end interface
+   pt_dimsr=c_loc(dimsr)
+   call STRCUT(name1024, name)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.4) then
+     write(hsmg,'(50hhdf5_read_data_4d_real8: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(rdata(dimsr(1),dimsr(2),dimsr(3),dimsr(4)))
+   pt_data=c_loc(rdata)
+   call hdf5_read_data_real8_c(ifile, name1024, pt_data)
+end subroutine hdf5_read_data_4d_real8
+!
+subroutine hdf5_read_data_0d_string(ifile, name, idata)
+   ! read a rank 0 string dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   character(len=*), target :: idata
+   character(kind=c_char), allocatable, dimension(:) :: pt_data
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_string_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp, idata
+      end subroutine hdf5_read_data_string_c
+   end interface
+   call STRCUT(name1024, name)
+   pt_dimsr=c_loc(dimsr)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if((rank.ne.1).or.(dimsr(1).ne.1)) then
+     write(hsmg,'(51hhdf5_read_data_0d_string: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   allocate(pt_data(nbyte+1)) ! add one byte for null termination
+   call hdf5_read_data_string_c(ifile, name1024, pt_data)
+   call STRFIL(idata,pt_data,nbyte)
+   deallocate(pt_data)
+end subroutine hdf5_read_data_0d_string
+!
+subroutine hdf5_read_data_1d_string(ifile, name, idata)
+   ! read a rank 1 string dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*), intent(in) :: name
+   character(len=*), allocatable, dimension(:), target :: idata
+   character(kind=c_char), allocatable, dimension(:) :: pt_data
+   !
+   integer :: rank, type, nbyte
+   type(c_ptr) :: pt_dimsr
+   integer,target,dimension(5) :: dimsr
+   interface
+      subroutine hdf5_info_c(ifile, namp, rank, type, nbyte, dimsr) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int) :: rank, type, nbyte
+         type(c_ptr), value :: dimsr
+      end subroutine hdf5_info_c
+   end interface
+   interface
+      subroutine hdf5_read_data_string_c(ifile, namp, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp, idata
+      end subroutine hdf5_read_data_string_c
+   end interface
+   call STRCUT(name1024, name)
+   pt_dimsr=c_loc(dimsr)
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   if(rank.ne.1) then
+     write(hsmg,'(51hhdf5_read_data_1d_string: invalid info for dataset ,a)')  name
+     call XABORT(hsmg)
+   endif
+   call hdf5_info_c(ifile, name1024, rank, type, nbyte, pt_dimsr)
+   allocate(idata(dimsr(1)))
+   allocate(pt_data(dimsr(1)*nbyte+1)) ! add one byte for null termination
+   call hdf5_read_data_string_c(ifile, name1024, pt_data)
+   call STRFIL1D(idata,pt_data,nbyte)
+   deallocate(pt_data)
+end subroutine hdf5_read_data_1d_string
+!
+subroutine hdf5_write_data_0d_int4(ifile, name, idata)
+   ! write a rank 1 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer,intent(in), target :: idata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   interface
+      subroutine hdf5_write_data_int_c(ifile, namp, rank, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_int_c
+   end interface
+   dimsr(1)=1
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_int_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_0d_int4
+!
+subroutine hdf5_write_data_1d_int4(ifile, name, idata)
+   ! write a rank 1 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer,dimension(:),intent(in), target :: idata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer, pointer :: idata_p
+   interface
+      subroutine hdf5_write_data_int_c(ifile, namp, rank, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_int_c
+   end interface
+   dimsr(1)=size(idata,1)
+   idata_p => idata(1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_int_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_1d_int4
+!
+subroutine hdf5_write_data_2d_int4(ifile, name, idata)
+   ! write a rank 2 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer,dimension(:,:),intent(in), target :: idata
+   !
+   integer, parameter :: rank = 2
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer, pointer :: idata_p
+   interface
+      subroutine hdf5_write_data_int_c(ifile, namp, rank, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_int_c
+   end interface
+   dimsr(2)=size(idata,1)
+   dimsr(1)=size(idata,2)
+   idata_p => idata(1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_int_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_2d_int4
+!
+subroutine hdf5_write_data_3d_int4(ifile, name, idata)
+   ! write a rank 3 integer dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   integer,dimension(:,:,:),intent(in), target :: idata
+   !
+   integer, parameter :: rank = 3
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   integer, pointer :: idata_p
+   interface
+      subroutine hdf5_write_data_int_c(ifile, namp, rank, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_int_c
+   end interface
+   dimsr(3)=size(idata,1)
+   dimsr(2)=size(idata,2)
+   dimsr(1)=size(idata,3)
+   idata_p => idata(1,1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_int_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_3d_int4
+!
+subroutine hdf5_write_data_0d_real4(ifile, name, rdata)
+   ! write a rank 0 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   interface
+      subroutine hdf5_write_data_real4_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real4_c
+   end interface
+   dimsr(1)=1
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real4_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_0d_real4
+!
+subroutine hdf5_write_data_1d_real4(ifile, name, rdata)
+   ! write a rank 1 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4),dimension(:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real, pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real4_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real4_c
+   end interface
+   dimsr(1)=size(rdata,1)
+   rdata_p => rdata(1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real4_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_1d_real4
+!
+subroutine hdf5_write_data_2d_real4(ifile, name, rdata)
+   ! write a rank 2 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4),dimension(:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 2
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real, pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real4_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real4_c
+   end interface
+   dimsr(2)=size(rdata,1)
+   dimsr(1)=size(rdata,2)
+   rdata_p => rdata(1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real4_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_2d_real4
+!
+subroutine hdf5_write_data_3d_real4(ifile, name, rdata)
+   ! write a rank 3 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4),dimension(:,:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 3
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real, pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real4_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real4_c
+   end interface
+   dimsr(3)=size(rdata,1)
+   dimsr(2)=size(rdata,2)
+   dimsr(1)=size(rdata,3)
+   rdata_p => rdata(1,1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real4_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_3d_real4
+!
+subroutine hdf5_write_data_4d_real4(ifile, name, rdata)
+   ! write a rank 4 real4 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(4),dimension(:,:,:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 4
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real, pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real4_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real4_c
+   end interface
+   dimsr(4)=size(rdata,1)
+   dimsr(3)=size(rdata,2)
+   dimsr(2)=size(rdata,3)
+   dimsr(1)=size(rdata,4)
+   rdata_p => rdata(1,1,1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real4_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_4d_real4
+!
+subroutine hdf5_write_data_0d_real8(ifile, name, rdata)
+   ! write a rank 0 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   interface
+      subroutine hdf5_write_data_real8_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real8_c
+   end interface
+   dimsr(1)=1
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real8_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_0d_real8
+!
+subroutine hdf5_write_data_1d_real8(ifile, name, rdata)
+   ! write a rank 1 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8),dimension(:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real(8), pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real8_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real8_c
+   end interface
+   dimsr(1)=size(rdata,1)
+   rdata_p => rdata(1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real8_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_1d_real8
+!
+subroutine hdf5_write_data_2d_real8(ifile, name, rdata)
+   ! write a rank 2 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8),dimension(:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 2
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real(8), pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real8_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real8_c
+   end interface
+   dimsr(2)=size(rdata,1)
+   dimsr(1)=size(rdata,2)
+   rdata_p => rdata(1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real8_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_2d_real8
+!
+subroutine hdf5_write_data_3d_real8(ifile, name, rdata)
+   ! write a rank 3 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8),dimension(:,:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 3
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real(8), pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real8_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real8_c
+   end interface
+   dimsr(3)=size(rdata,1)
+   dimsr(2)=size(rdata,2)
+   dimsr(1)=size(rdata,3)
+   rdata_p => rdata(1,1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real8_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_3d_real8
+!
+subroutine hdf5_write_data_4d_real8(ifile, name, rdata)
+   ! write a rank 4 real8 dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   real(8),dimension(:,:,:,:),intent(in), target :: rdata
+   !
+   integer, parameter :: rank = 4
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   real(8), pointer :: rdata_p
+   interface
+      subroutine hdf5_write_data_real8_c(ifile, namp, rank, dimsf, rdata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank
+         type(c_ptr), value :: dimsf, rdata
+      end subroutine hdf5_write_data_real8_c
+   end interface
+   dimsr(4)=size(rdata,1)
+   dimsr(3)=size(rdata,2)
+   dimsr(2)=size(rdata,3)
+   dimsr(1)=size(rdata,4)
+   rdata_p => rdata(1,1,1,1)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(rdata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_real8_c(ifile, name1024, rank, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_4d_real8
+!
+subroutine hdf5_write_data_0d_string(ifile, name, idata)
+   ! write a rank 1 string dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   character(len=*),intent(in), target :: idata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   character(len=1), pointer :: idata_p
+   interface
+      subroutine hdf5_write_data_string_c(ifile, namp, rank, lenstr, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank, lenstr
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_string_c
+   end interface
+   dimsr(1)=1
+   idata_p => idata
+   lenstr=len(idata)
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_string_c(ifile, name1024, rank, lenstr, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_0d_string
+!
+subroutine hdf5_write_data_1d_string(ifile, name, idata)
+   ! write a rank 1 string dataset
+   type(c_ptr),intent(in) :: ifile
+   character(len=*),intent(in) :: name
+   character(len=*),intent(in), dimension(:), target :: idata
+   !
+   integer, parameter :: rank = 1
+   integer,dimension(rank),target :: dimsr
+   type(c_ptr) :: pt_dimsr, pt_data
+   character(len=1), pointer :: idata_p
+   interface
+      subroutine hdf5_write_data_string_c(ifile, namp, rank, lenstr, dimsf, idata) bind(c)
+         use, intrinsic :: iso_c_binding
+         type(c_ptr) :: ifile
+         character(kind=c_char), dimension(*) :: namp
+         integer(c_int), value :: rank, lenstr
+         type(c_ptr), value :: dimsf, idata
+      end subroutine hdf5_write_data_string_c
+   end interface
+   dimsr(1)=size(idata,1)
+   idata_p => idata(1)
+   lenstr=len(idata(1))
+   pt_dimsr=c_loc(dimsr)
+   pt_data=c_loc(idata_p)
+   call STRCUT(name1024, name)
+   call hdf5_write_data_string_c(ifile, name1024, rank, lenstr, pt_dimsr, pt_data)
+end subroutine hdf5_write_data_1d_string
+end module hdf5_wrap