gtool_file.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
! Warning: This file was generated from gtool_file.f90.erb.
!          Do not edit this file.
!-------------------------------------------------------------------------------
module gtool_file
  !-----------------------------------------------------------------------------
  !
  !++ Used modules
  !
  use gtool_file_h
  use dc_log, only: &
       log, &
       log_lmsg
  use dc_types, only: &
       dp, &
       sp
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedures
  !
  public :: fileopen
  public :: filecreate
  public :: filesetoption
  public :: filedefaxis
  public :: fileputaxis
  public :: filewriteaxis
  public :: filedefassociatedcoordinates
  public :: fileputassociatedcoordinates
  public :: filewriteassociatedcoordinates
  public :: fileaddvariable
  public :: filedefinevariable
  public :: filesettattr
  public :: filegetshape
  public :: filegetdatainfo
  public :: filegetalldatainfo
  public :: fileread
  public :: filewrite
  public :: filegetglobalattribute
  public :: filesetglobalattribute
  public :: fileenddef
  public :: fileflush
  public :: fileclose
  public :: filecloseall
  public :: filemakefname
  public :: fileattachbuffer
  public :: filedetachbuffer

  interface filegetdatainfo
     module procedure filegetdatainfofid
     module procedure filegetdatainfofname
  end interface filegetdatainfo
  interface filegetalldatainfo
     module procedure filegetalldatainfofid
     module procedure filegetalldatainfofname
  end interface filegetalldatainfo

  interface fileputaxis
     module procedure fileputaxisrealsp
     module procedure fileputaxisrealdp
  end interface fileputaxis
  interface filewriteaxis
     module procedure filewriteaxisrealsp
     module procedure filewriteaxisrealdp
  end interface filewriteaxis
  interface fileputassociatedcoordinates
    module procedure fileput1dassociatedcoordinatesrealsp
    module procedure fileput1dassociatedcoordinatesrealdp
    module procedure fileput2dassociatedcoordinatesrealsp
    module procedure fileput2dassociatedcoordinatesrealdp
    module procedure fileput3dassociatedcoordinatesrealsp
    module procedure fileput3dassociatedcoordinatesrealdp
    module procedure fileput4dassociatedcoordinatesrealsp
    module procedure fileput4dassociatedcoordinatesrealdp
  end interface fileputassociatedcoordinates
  interface filewriteassociatedcoordinates
    module procedure filewrite1dassociatedcoordinatesrealsp
    module procedure filewrite1dassociatedcoordinatesrealdp
    module procedure filewrite2dassociatedcoordinatesrealsp
    module procedure filewrite2dassociatedcoordinatesrealdp
    module procedure filewrite3dassociatedcoordinatesrealsp
    module procedure filewrite3dassociatedcoordinatesrealdp
    module procedure filewrite4dassociatedcoordinatesrealsp
    module procedure filewrite4dassociatedcoordinatesrealdp
  end interface filewriteassociatedcoordinates
  interface fileaddvariable
     module procedure fileaddvariablenot
     module procedure fileaddvariablerealsp
     module procedure fileaddvariablerealdp
  end interface fileaddvariable
  interface fileread
    module procedure fileread1drealsp
    module procedure fileread1drealdp
    module procedure fileread2drealsp
    module procedure fileread2drealdp
    module procedure fileread3drealsp
    module procedure fileread3drealdp
    module procedure fileread4drealsp
    module procedure fileread4drealdp
    module procedure filereadvar1drealsp
    module procedure filereadvar1drealdp
    module procedure filereadvar2drealsp
    module procedure filereadvar2drealdp
    module procedure filereadvar3drealsp
    module procedure filereadvar3drealdp
    module procedure filereadvar4drealsp
    module procedure filereadvar4drealdp
  end interface fileread
  interface filewrite
    module procedure filewrite1drealsp
    module procedure filewrite1drealdp
    module procedure filewrite2drealsp
    module procedure filewrite2drealdp
    module procedure filewrite3drealsp
    module procedure filewrite3drealdp
    module procedure filewrite4drealsp
    module procedure filewrite4drealdp
  end interface filewrite
  interface filegetglobalattribute
     module procedure filegetglobalattributetext
     module procedure filegetglobalattributeint
     module procedure filegetglobalattributefloat
     module procedure filegetglobalattributedouble
  end interface filegetglobalattribute
  interface filesetglobalattribute
     module procedure filesetglobalattributetext
     module procedure filesetglobalattributeint
     module procedure filesetglobalattributefloat
     module procedure filesetglobalattributedouble
  end interface filesetglobalattribute

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(DP), parameter, public :: rmiss = -9.9999e+30
  !-----------------------------------------------------------------------------
  !
  !++ Private procedures
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  integer,                   private, parameter :: file_nfile_max = 512   ! number limit of file
                                                 ! Keep consistency with "FILE_MAX" in gtool_netcdf.c
  integer,                   private, parameter :: file_nvar_max  = 40960 ! number limit of variables
                                                 ! Keep consistency with "VAR_MAX" in gtool_netcdf.c

  character(len=File_HLONG), private,      save :: file_fname_list(file_nfile_max)
  integer,                   private,      save :: file_fid_list  (file_nfile_max)
  integer,                   private,      save :: file_fid_count = 1
  character(len=File_HLONG), private,      save :: file_vname_list  (file_nvar_max)
  integer,                   private,      save :: file_vid_fid_list(file_nvar_max)
  integer,                   private,      save :: file_vid_list    (file_nvar_max)
  integer,                   private,      save :: file_vid_count = 1
  integer,                   private,      save :: mpi_myrank

  character(len=LOG_LMSG),   private            :: message

contains
  !-----------------------------------------------------------------------------
  subroutine filecreate( &
       fid,         & ! (out)
       existed,     & ! (out)
       basename,    & ! (in)
       title,       & ! (in)
       source,      & ! (in)
       institution, & ! (in)
       master,      & ! (in)
       myrank,      & ! (in)
       rankidx,     & ! (in)
       single,      & ! (in) optional
       time_units,  & ! (in) optional
       append,      & ! (in) optional
       comm         ) ! (in) optional
    use mpi, only : mpi_comm_null
    implicit none

    integer,          intent(out)           :: fid
    logical,          intent(out)           :: existed
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: title
    character(len=*), intent( in)           :: source
    character(len=*), intent( in)           :: institution
    integer,          intent( in)           :: master
    integer,          intent( in)           :: myrank
    integer,          intent( in)           :: rankidx(:)
    character(len=*), intent( in), optional :: time_units
    logical,          intent( in), optional :: single
    logical,          intent( in), optional :: append
    integer,          intent( in), optional :: comm ! MPI communicator

    character(len=File_HMID) :: time_units_
    logical :: single_
    integer :: mode
    integer :: error

    intrinsic size

    if ( present(time_units) ) then
       time_units_ = time_units
    else
       time_units_ = 'seconds'
    end if

    mpi_myrank = myrank

    if ( present(single) ) then
       if ( single .and. (myrank .ne. master) ) return
       single_ = single
    else
       single_ = .false.
    endif

    mode = file_fwrite
    if ( present(append) ) then
       if ( append ) mode = file_fappend
    end if

    call filegetfid( &
         fid,        & ! (out)
         existed,    & ! (out)
         basename,   & ! (in)
         mode,       & ! (in)
         single_,    & ! (in)
         comm        & ! (in)
         )

    if ( existed ) return

    !--- append package header to the file
    call filesetglobalattribute( fid, & ! (in)
         "title", title               ) ! (in)
    call filesetglobalattribute( fid, & ! (in)
         "source", source             ) ! (in)
    call filesetglobalattribute( fid, & ! (in)
         "institution", institution   ) ! (in)

    if ( .NOT. present(comm) .OR. comm == mpi_comm_null ) then
       ! for shared-file parallel I/O, skip attributes related to MPI processes
       call filesetglobalattribute( fid, & ! (in)
            "myrank", (/myrank/)         ) ! (in)
       call filesetglobalattribute( fid, & ! (in)
            "rankidx", rankidx           ) ! (in)
    end if

    call file_set_tunits( fid, & ! (in)
         time_units_,          & ! (in)
         error                 ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set time units')
    end if

    return
  end subroutine filecreate

  !-----------------------------------------------------------------------------
  subroutine filegetglobalattributetext( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (out)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    character(len=*), intent(out) :: val

    integer error

    intrinsic size

    call file_get_global_attribute_text( & ! (in)
         fid, key,                       & ! (in)
         val, error                      ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get text global attribute: '//trim(key))
    end if

    return
  end subroutine filegetglobalattributetext

  !-----------------------------------------------------------------------------
  subroutine filegetglobalattributeint( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (out)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    integer,          intent(out) :: val(:)

    integer error

    intrinsic size

    call file_get_global_attribute_int( & ! (in)
         fid, key, size(val),           & ! (in)
         val, error                     ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get integer global attribute: '//trim(key))
    end if

    return
  end subroutine filegetglobalattributeint

  !-----------------------------------------------------------------------------
  subroutine filegetglobalattributefloat( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (out)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    real(SP),    intent(out) :: val(:)

    integer error

    intrinsic size

    call file_get_global_attribute_float( & ! (in)
         fid, key, size(val),                          & ! (in)
         val, error                                    ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get float global attribute: '//trim(key))
    end if

    return
  end subroutine filegetglobalattributefloat

  !-----------------------------------------------------------------------------
  subroutine filegetglobalattributedouble( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (out)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    real(DP),    intent(out) :: val(:)

    integer error

    intrinsic size

    call file_get_global_attribute_double( & ! (in)
         fid, key, size(val),                          & ! (in)
         val, error                                    ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get double global attribute: '//trim(key))
    end if

    return
  end subroutine filegetglobalattributedouble


  !-----------------------------------------------------------------------------
  subroutine filesetglobalattributetext( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (in)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    character(len=*), intent(in) :: val

    integer error

    call file_set_global_attribute_text( fid, & ! (in)
         key, val,                            & ! (in)
         error                                ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set text global attribute: '//trim(key))
    end if

    return
  end subroutine filesetglobalattributetext

  !-----------------------------------------------------------------------------
  subroutine filesetglobalattributeint( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (in)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    integer,          intent(in) :: val(:)

    integer error

    intrinsic size

    call file_set_global_attribute_int( fid, & ! (in)
         key, val, size(val),                & ! (in)
         error                               ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set integer global attribute: '//trim(key))
    end if

    return
  end subroutine filesetglobalattributeint

  !-----------------------------------------------------------------------------
  subroutine filesetglobalattributefloat( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (in)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    real(SP),    intent(in) :: val(:)

    integer error

    intrinsic size

    call file_set_global_attribute_float( fid, & ! (in)
         key, val, size(val),                               & ! (in)
         error                                              ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set float global attribute: '//trim(key))
    end if

    return
  end subroutine filesetglobalattributefloat

  !-----------------------------------------------------------------------------
  subroutine filesetglobalattributedouble( &
       fid,      & ! (in)
       key,      & ! (in)
       val       & ! (in)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: key
    real(DP),    intent(in) :: val(:)

    integer error

    intrinsic size

    call file_set_global_attribute_double( fid, & ! (in)
         key, val, size(val),                               & ! (in)
         error                                              ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set double global attribute: '//trim(key))
    end if

    return
  end subroutine filesetglobalattributedouble

  !-----------------------------------------------------------------------------
  subroutine filesetoption( &
       fid,      & ! (in)
       filetype, & ! (in)
       key,      & ! (in)
       val       & ! (in)
       )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: filetype
    character(len=*), intent(in) :: key
    character(len=*), intent(in) :: val

    integer error

    call file_set_option( fid, filetype, key, val, & ! (in)
                          error                    ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to set option')
    end if

    return
  end subroutine filesetoption

  !-----------------------------------------------------------------------------
  subroutine fileopen( &
      fid,       & ! (out)
      basename,  & ! (in)
      mode,      & ! (in)
      single,    & ! (in) optional
      comm,      & ! (in) optional
      myrank     & ! (in) optional
      )
    implicit none

    integer,          intent(out) :: fid
    character(len=*), intent( in) :: basename
    integer,          intent( in) :: mode
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: comm
    integer,          intent( in), optional :: myrank

    logical :: existed
    logical :: single_

    single_ = .false.

    if ( present(single) ) single_ = single
    if ( present(myrank) ) mpi_myrank = myrank

    call filegetfid( fid, existed,     & ! (out)
         basename, mode, single_, comm ) ! (in)

    return
  end subroutine fileopen

  !-----------------------------------------------------------------------------
  ! interface FilePutAxis
  !-----------------------------------------------------------------------------
  subroutine fileputaxisrealsp( &
       fid,      & ! (in)
       name,     & ! (in)
       desc,     & ! (in)
       units,    & ! (in)
       dim_name, & ! (in)
       dtype,    & ! (in)
       val       ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_name
    integer,          intent(in) :: dtype
    real(SP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_put_axis( fid, name, desc, units, dim_name, dtype, val, size(val), sp, & ! (in)
         error                                                                            ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put axis')
    end if

    return
  end subroutine fileputaxisrealsp
  subroutine fileputaxisrealdp( &
       fid,      & ! (in)
       name,     & ! (in)
       desc,     & ! (in)
       units,    & ! (in)
       dim_name, & ! (in)
       dtype,    & ! (in)
       val       ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_name
    integer,          intent(in) :: dtype
    real(DP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_put_axis( fid, name, desc, units, dim_name, dtype, val, size(val), dp, & ! (in)
         error                                                                            ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put axis')
    end if

    return
  end subroutine fileputaxisrealdp

  subroutine filedefaxis( &
       fid,      & ! (in)
       name,     & ! (in)
       desc,     & ! (in)
       units,    & ! (in)
       dim_name, & ! (in)
       dtype,    & ! (in)
       dim_size  ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_name
    integer,          intent(in) :: dtype
    integer,          intent(in) :: dim_size

    integer error
    intrinsic size

    call file_def_axis( fid, name, desc, units, dim_name, dtype, dim_size, & ! (in)
         error                                                             ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to define axis')
    end if

    return
  end subroutine filedefaxis

  !-----------------------------------------------------------------------------
  ! interface FileWriteAxis
  !-----------------------------------------------------------------------------
  subroutine filewriteaxisrealsp( &
       fid,      & ! (in)
       name,     & ! (in)
       val,      & ! (in)
       start     ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(SP),    intent(in)           :: val(:)
    integer,          intent(in), optional :: start(:)

    integer error
    intrinsic shape

    if ( present(start) ) then
       call file_write_axis( fid, name, val, sp, start, shape(val), & ! (in)
            error                                                        ) ! (out)
    else
       call file_write_axis( fid, name, val, sp, (/1/), shape(val), & ! (in)
            error                                                        ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to write axis')
    end if

    return
  end subroutine filewriteaxisrealsp
  subroutine filewriteaxisrealdp( &
       fid,      & ! (in)
       name,     & ! (in)
       val,      & ! (in)
       start     ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(DP),    intent(in)           :: val(:)
    integer,          intent(in), optional :: start(:)

    integer error
    intrinsic shape

    if ( present(start) ) then
       call file_write_axis( fid, name, val, dp, start, shape(val), & ! (in)
            error                                                        ) ! (out)
    else
       call file_write_axis( fid, name, val, dp, (/1/), shape(val), & ! (in)
            error                                                        ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to write axis')
    end if

    return
  end subroutine filewriteaxisrealdp

  !-----------------------------------------------------------------------------
  ! interface FilePutAssociatedCoordinates
  !-----------------------------------------------------------------------------
  subroutine fileput1dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(SP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, sp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput1dassociatedcoordinatesrealsp
  subroutine fileput1dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(DP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, dp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput1dassociatedcoordinatesrealdp
  subroutine fileput2dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(SP),    intent(in) :: val(:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, sp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput2dassociatedcoordinatesrealsp
  subroutine fileput2dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(DP),    intent(in) :: val(:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, dp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput2dassociatedcoordinatesrealdp
  subroutine fileput3dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(SP),    intent(in) :: val(:,:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, sp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput3dassociatedcoordinatesrealsp
  subroutine fileput3dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(DP),    intent(in) :: val(:,:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, dp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput3dassociatedcoordinatesrealdp
  subroutine fileput4dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(SP),    intent(in) :: val(:,:,:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, sp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput4dassociatedcoordinatesrealsp
  subroutine fileput4dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype,     & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype
    real(DP),    intent(in) :: val(:,:,:,:)

    integer error
    intrinsic size

    call file_put_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         val, dp,                                         & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine fileput4dassociatedcoordinatesrealdp

  subroutine filedefassociatedcoordinates( &
       fid,       & ! (in)
       name,      & ! (in)
       desc,      & ! (in)
       units,     & ! (in)
       dim_names, & ! (in)
       dtype      ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    character(len=*), intent(in) :: desc
    character(len=*), intent(in) :: units
    character(len=*), intent(in) :: dim_names(:)
    integer,          intent(in) :: dtype

    integer error
    intrinsic size

    call file_def_associated_coordinates( fid,                 & ! (in)
         name, desc, units, dim_names, size(dim_names), dtype, & ! (in)
         error                                                 ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates')
    end if

    return
  end subroutine filedefassociatedcoordinates

  !-----------------------------------------------------------------------------
  ! interface FileWriteAssociatedCoordinates
  !-----------------------------------------------------------------------------
  subroutine filewrite1dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(SP),    intent(in)           :: val(:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            1, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            1, (/1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite1dassociatedcoordinatesrealsp
  subroutine filewrite1dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(DP),    intent(in)           :: val(:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            1, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            1, (/1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite1dassociatedcoordinatesrealdp
  subroutine filewrite2dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(SP),    intent(in)           :: val(:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            2, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            2, (/1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite2dassociatedcoordinatesrealsp
  subroutine filewrite2dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(DP),    intent(in)           :: val(:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            2, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            2, (/1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite2dassociatedcoordinatesrealdp
  subroutine filewrite3dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(SP),    intent(in)           :: val(:,:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            3, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            3, (/1,1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite3dassociatedcoordinatesrealsp
  subroutine filewrite3dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(DP),    intent(in)           :: val(:,:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            3, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            3, (/1,1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite3dassociatedcoordinatesrealdp
  subroutine filewrite4dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(SP),    intent(in)           :: val(:,:,:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            4, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, sp, & ! (in)
            4, (/1,1,1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite4dassociatedcoordinatesrealsp
  subroutine filewrite4dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val,       & ! (in)
       start,     & ! (in)
       count,     & ! (in)
       ndims      ) ! (in)
    integer,          intent(in)           :: fid
    character(len=*), intent(in)           :: name
    real(DP),    intent(in)           :: val(:,:,:,:)
    integer,          intent(in), optional :: start(:)
    integer,          intent(in), optional :: count(:)  ! in case val has been reshaped
    integer,          intent(in), optional :: ndims     ! in case val has been reshaped

    integer error
    intrinsic size, shape

    if ( present(ndims) ) then
       ! Note this is called for history coordinates which have been reshaped
       ! from 2D/3D into 1D array. In this case, start and count must be also present
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            ndims, start, count,                                        & ! (in)
            error                                                       ) ! (out)
    else if ( present(start) ) then
       ! Note this is called for restart coordinates
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            4, start, shape(val),                                  & ! (in)
            error                                                       ) ! (out)
    else
       ! Note this is for the one-file-per-process I/O method
       call file_write_associated_coordinates( fid, name, val, dp, & ! (in)
            4, (/1,1,1,1/), shape(val),             & ! (in)
            error                                                       ) ! (out)
    end if
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to put associated coordinates: '//trim(name))
    end if

    return
  end subroutine filewrite4dassociatedcoordinatesrealdp

  !-----------------------------------------------------------------------------
  ! interface FileAddVariable
  !-----------------------------------------------------------------------------
  subroutine fileaddvariablenot( &
       vid,     & ! (out)
       fid,     & ! (in)
       varname, & ! (in)
       desc,    & ! (in)
       units,   & ! (in)
       dims,    & ! (in)
       dtype,   & ! (in)
       tavg     & ! (in) optional
       )
    integer,          intent(out) :: vid
    integer,          intent( in) :: fid
    character(len=*), intent( in) :: varname
    character(len=*), intent( in) :: desc
    character(len=*), intent( in) :: units
    character(len=*), intent( in) :: dims(:)
    integer,          intent( in) :: dtype
    logical,          intent( in), optional :: tavg

    call fileaddvariablerealdp(vid, fid, varname, desc, units, dims, dtype, &
         -1.0_dp, tavg )

    return
  end subroutine fileaddvariablenot
  subroutine fileaddvariablerealsp( &
       vid,     & ! (out)
       fid,     & ! (in)
       varname, & ! (in)
       desc,    & ! (in)
       units,   & ! (in)
       dims,    & ! (in)
       dtype,   & ! (in)
       tint,    & ! (in)
       tavg     & ! (in) optional
       )
    integer,          intent(out) :: vid
    integer,          intent( in) :: fid
    character(len=*), intent( in) :: varname
    character(len=*), intent( in) :: desc
    character(len=*), intent( in) :: units
    character(len=*), intent( in) :: dims(:)
    integer,          intent( in) :: dtype
    real(SP),    intent( in) :: tint
    logical,          intent( in), optional :: tavg

    real(DP) :: tint8
    integer  :: itavg
    integer  :: error
    integer  :: n

    intrinsic size
    !---------------------------------------------------------------------------

    vid = -1
    do n = 1, file_vid_count
       if ( file_vid_fid_list(n) == fid .and. &
            varname == file_vname_list(n) ) then
          vid = file_vid_list(n)
       end if
    enddo

    if ( vid < 0 ) then ! variable registration
       !--- register new variable
       write(message,'(2A)') '###### Variable registration : name = ', trim(varname)
       call log("I",message)

       tint8 = real(tint,dp)

       if ( present(tavg) ) then
          if ( tavg ) then
             itavg = 1
          else
             itavg = 0
          end if
       else
          itavg = 0
       end if

       call file_add_variable( vid,                             & ! (out)
            fid, varname, desc, units, dims, size(dims), dtype, & ! (in)
            tint8, itavg,                                       & ! (in)
            error                                               ) ! (out)
       if ( error /= success_code ) then
          call log('E', 'xxx failed to add variable: '//trim(varname))
       end if

       file_vname_list(file_vid_count) = trim(varname)
       file_vid_list(file_vid_count) = vid
       file_vid_fid_list(file_vid_count) = fid
       file_vid_count                    = file_vid_count + 1
    endif

    return
  end subroutine fileaddvariablerealsp
  subroutine fileaddvariablerealdp( &
       vid,     & ! (out)
       fid,     & ! (in)
       varname, & ! (in)
       desc,    & ! (in)
       units,   & ! (in)
       dims,    & ! (in)
       dtype,   & ! (in)
       tint,    & ! (in)
       tavg     & ! (in) optional
       )
    integer,          intent(out) :: vid
    integer,          intent( in) :: fid
    character(len=*), intent( in) :: varname
    character(len=*), intent( in) :: desc
    character(len=*), intent( in) :: units
    character(len=*), intent( in) :: dims(:)
    integer,          intent( in) :: dtype
    real(DP),    intent( in) :: tint
    logical,          intent( in), optional :: tavg

    real(DP) :: tint8
    integer  :: itavg
    integer  :: error
    integer  :: n

    intrinsic size
    !---------------------------------------------------------------------------

    vid = -1
    do n = 1, file_vid_count
       if ( file_vid_fid_list(n) == fid .and. &
            varname == file_vname_list(n) ) then
          vid = file_vid_list(n)
       end if
    enddo

    if ( vid < 0 ) then ! variable registration
       !--- register new variable
       write(message,'(2A)') '###### Variable registration : name = ', trim(varname)
       call log("I",message)

       tint8 = real(tint,dp)

       if ( present(tavg) ) then
          if ( tavg ) then
             itavg = 1
          else
             itavg = 0
          end if
       else
          itavg = 0
       end if

       call file_add_variable( vid,                             & ! (out)
            fid, varname, desc, units, dims, size(dims), dtype, & ! (in)
            tint8, itavg,                                       & ! (in)
            error                                               ) ! (out)
       if ( error /= success_code ) then
          call log('E', 'xxx failed to add variable: '//trim(varname))
       end if

       file_vname_list(file_vid_count) = trim(varname)
       file_vid_list(file_vid_count) = vid
       file_vid_fid_list(file_vid_count) = fid
       file_vid_count                    = file_vid_count + 1
    endif

    return
  end subroutine fileaddvariablerealdp

  subroutine filedefinevariable( &
       fid,     & ! (in)
       vid,     & ! (out)
       varname, & ! (in)
       desc,    & ! (in)
       units,   & ! (in)
       ndims,   & ! (in)
       dims,    & ! (in)
       dtype,   & ! (in)
       tint,    & ! (in) optional
       tavg     & ! (in) optional
       )
    integer,          intent(out) :: vid
    integer,          intent( in) :: fid
    character(len=*), intent( in) :: varname
    character(len=*), intent( in) :: desc
    character(len=*), intent( in) :: units
    integer,          intent( in) :: ndims
    character(len=*), intent( in) :: dims(:)
    integer,          intent( in) :: dtype
    real(DP),         intent( in), optional :: tint
    logical,          intent( in), optional :: tavg

    real(DP) :: tint_
    integer  :: itavg
    integer  :: error
    integer  :: n

    intrinsic size
    !---------------------------------------------------------------------------

    vid = -1
    do n = 1, file_vid_count
       if ( file_vid_fid_list(n) == fid .and. &
            varname == file_vname_list(n) ) then
          vid = file_vid_list(n)
       end if
    enddo

    if ( vid < 0 ) then ! variable registration
       !--- register new variable
       write(message,'(2A)') '###### Variable registration : name = ', trim(varname)
       call log("I",message)

       if ( present(tint) ) then
          tint_ = tint
       else
          tint_ = -1.0_dp
       endif

       if ( present(tavg) ) then
          if ( tavg ) then
             itavg = 1
          else
             itavg = 0
          end if
       else
          itavg = 0
       end if

       call file_add_variable( vid,                        & ! (out)
            fid, varname, desc, units, dims, ndims, dtype, & ! (in)
            tint_, itavg,                                  & ! (in)
            error                                          ) ! (out)
       if ( error /= success_code ) then
          call log('E', 'xxx failed to add variable: '//trim(varname))
       end if

       file_vname_list(file_vid_count) = trim(varname)
       file_vid_list(file_vid_count) = vid
       file_vid_fid_list(file_vid_count) = fid
       file_vid_count                    = file_vid_count + 1
    endif

    return
  end subroutine filedefinevariable

  !-----------------------------------------------------------------------------
  ! FileSetTAttr
  !-----------------------------------------------------------------------------
  subroutine filesettattr( &
     fid,   & ! (in)
     vname, & ! (in)
     key,   & ! (in)
     val    & ! (in)
     )
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: vname
    character(len=*), intent(in) :: key
    character(len=*), intent(in) :: val

    integer :: error

    call file_set_tattr( &
         fid, vname, & ! (in)
         key, val,   & ! (in)
         error       ) ! (out)
    if ( error /= success_code .and. error /= already_existed_code ) then
       call log('E', 'xxx failed to set attr for axis')
    end if

    return
  end subroutine filesettattr

  !-----------------------------------------------------------------------------
  ! FileGetShape
  !-----------------------------------------------------------------------------
  subroutine filegetshape( &
      dims,          & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      myrank,        & ! (in)
      single,        & ! (in) optional
      error          ) ! (out) optional
    implicit none

    integer,          intent(out)           :: dims(:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: single
    logical,          intent(out), optional :: error

    integer :: fid
    type(datainfo) :: dinfo
    integer :: ierror
    integer :: n

    logical :: single_
    logical :: suppress

    intrinsic size
    intrinsic shape
    !---------------------------------------------------------------------------

    mpi_myrank = myrank

    if ( present(single) ) then
       single_ = single
    else
       single_ = .false.
    end if

    if ( present(error) ) then
       suppress = .true.
    else
       suppress = .false.
    end if

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,  & ! (out)
         fid, varname, 1, suppress, & ! (in)
         ierror                     ) ! (out)

    !--- verify
    if ( ierror /= success_code ) then
       if ( present(error) ) then
          error = .true.
          return
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= size(dims) ) then
       write(message,*) 'xxx rank is different, ', size(dims), dinfo%rank
       call log('E', message)
    end if
    do n = 1, size(dims)
       dims(n) = dinfo%dim_size(n)
    end do

    if ( present(error) ) error = .false.

    return
  end subroutine filegetshape

  !-----------------------------------------------------------------------------
  ! FileGetData
  !-----------------------------------------------------------------------------
  subroutine filegetdatainfofname( &
      basename,    &
      varname,     &
      myrank,      &
      istep,       &
      single,      &
      description, &
      units,       &
      datatype,    &
      dim_rank,    &
      dim_name,    &
      dim_size,    &
      time_start,  &
      time_end,    &
      time_units   )
    implicit none
    character(len=*),           intent(in)  :: basename
    character(len=*),           intent(in)  :: varname
    integer,                    intent(in)  :: myrank
    integer,                    intent(in)  :: istep
    logical,                    intent(in), optional :: single

    character(len=File_HMID),   intent(out), optional :: description
    character(len=File_HSHORT), intent(out), optional :: units
    integer,                    intent(out), optional :: datatype
    integer,                    intent(out), optional :: dim_rank
    character(len=File_HSHORT), intent(out), optional :: dim_name(:)
    integer,                    intent(out), optional :: dim_size(:)
    real(DP),                   intent(out), optional :: time_start
    real(DP),                   intent(out), optional :: time_end
    character(len=File_HMID),   intent(out), optional :: time_units

    logical :: single_
    integer :: fid

    if ( present(single) ) then
       single_ = single
    else
       single_ = .false.
    end if

    mpi_myrank = myrank

    !--- search/register file
    call fileopen( fid,        & ! [OUT]
                   basename,   & ! [IN]
                   file_fread, & ! [IN]
                   single_     ) ! [IN]

    call filegetdatainfofid( &
         fid,         &
         varname,     &
         istep,       &
         description, &
         units,       &
         datatype,    &
         dim_rank,    &
         dim_name,    &
         dim_size,    &
         time_start,  &
         time_end,    &
         time_units   )

    return
  end subroutine filegetdatainfofname
  subroutine filegetdatainfofid( &
      fid,         &
      varname,     &
      istep,       &
      description, &
      units,       &
      datatype,    &
      dim_rank,    &
      dim_name,    &
      dim_size,    &
      time_start,  &
      time_end,    &
      time_units   )
    implicit none
    integer,                    intent(in)  :: fid
    character(len=*),           intent(in)  :: varname
    integer,                    intent(in)  :: istep

    character(len=File_HMID),   intent(out), optional :: description
    character(len=File_HSHORT), intent(out), optional :: units
    integer,                    intent(out), optional :: datatype
    integer,                    intent(out), optional :: dim_rank
    character(len=File_HSHORT), intent(out), optional :: dim_name(:)
    integer,                    intent(out), optional :: dim_size(:)
    real(DP),                   intent(out), optional :: time_start
    real(DP),                   intent(out), optional :: time_end
    character(len=File_HMID),   intent(out), optional :: time_units

    type(datainfo) :: dinfo

    integer :: ndim, idim
    real(DP):: time(1)

    integer :: error

    intrinsic size
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,   & ! [OUT]
                            fid,     & ! [IN]
                            varname, & ! [IN]
                            istep,   & ! [IN]
                            .false., & ! [IN]
                            error    ) ! [OUT]

    !--- verify and exit
    if ( error /= success_code ) then
       call log('E', 'xxx data info not found')
    endif

    if ( present(description) ) description = dinfo%description
    if ( present(units)       ) units       = dinfo%units
    if ( present(datatype)    ) datatype    = dinfo%datatype
    if ( present(dim_rank)    ) dim_rank    = dinfo%rank

    if ( present(dim_name) ) then
       ndim = min( dinfo%rank, size(dim_name) ) ! limit dimension rank
       do idim = 1, ndim
          dim_name(idim) = dinfo%dim_name(idim)
       enddo
    end if
    if ( present(dim_size) ) then
       ndim = min( dinfo%rank, size(dim_size) ) ! limit dimension rank
       do idim = 1, ndim
          dim_size(idim) = dinfo%dim_size(idim)
       enddo
    end if

    if ( present(time_units)  ) then
       if ( dinfo%time_units == "" ) then
          call filegetglobalattribute( fid, "time_units", time_units )
       else
          time_units  = dinfo%time_units
       end if
    end if
    if ( present(time_start)  ) then
       if ( dinfo%time_units == "" ) then
          call filegetglobalattribute( fid, "time", time )
          time_start = time(1)
       else
          time_start  = dinfo%time_start
       end if
    end if
    if ( present(time_end)  ) then
       if ( dinfo%time_units == "" ) then
          call filegetglobalattribute( fid, "time", time )
          time_end = time(1)
       else
          time_end  = dinfo%time_end
       end if
    end if

    return
  end subroutine filegetdatainfofid

  !-----------------------------------------------------------------------------
  ! FileGetData
  !-----------------------------------------------------------------------------
  subroutine filegetalldatainfofname( &
      step_limit,  &
      dim_limit,   &
      basename,    &
      varname,     &
      myrank,      &
      step_nmax,   &
      description, &
      units,       &
      datatype,    &
      dim_rank,    &
      dim_name,    &
      dim_size,    &
      time_start,  &
      time_end,    &
      time_units,  &
      single       )
    implicit none

    integer,                    intent(in)  :: step_limit
    integer,                    intent(in)  :: dim_limit
    character(len=*),           intent(in)  :: basename
    character(len=*),           intent(in)  :: varname
    integer,                    intent(in)  :: myrank
    integer,                    intent(out) :: step_nmax
    character(len=File_HMID),   intent(out) :: description
    character(len=File_HSHORT), intent(out) :: units
    integer,                    intent(out) :: datatype
    integer,                    intent(out) :: dim_rank
    character(len=File_HSHORT), intent(out) :: dim_name  (dim_limit)
    integer,                    intent(out) :: dim_size  (dim_limit)
    real(DP),                   intent(out) :: time_start(step_limit)
    real(DP),                   intent(out) :: time_end  (step_limit)
    character(len=File_HMID),   intent(out) :: time_units

    logical,                    intent(in), optional :: single


    integer :: fid
    logical :: single_

    mpi_myrank = myrank

    if ( present(single) ) then
       single_ = single
    else
       single_ = .false.
    end if

    !--- search/register file
    call fileopen( fid,        & ! [OUT]
                   basename,   & ! [IN]
                   file_fread, & ! [IN]
                   single_     ) ! [IN]

    call filegetalldatainfofid( &
      step_limit,  &
      dim_limit,   &
      fid,         &
      varname,     &
      step_nmax,   &
      description, &
      units,       &
      datatype,    &
      dim_rank,    &
      dim_name,    &
      dim_size,    &
      time_start,  &
      time_end,    &
      time_units   )

    return
  end subroutine filegetalldatainfofname
  subroutine filegetalldatainfofid( &
      step_limit,  &
      dim_limit,   &
      fid,         &
      varname,     &
      step_nmax,   &
      description, &
      units,       &
      datatype,    &
      dim_rank,    &
      dim_name,    &
      dim_size,    &
      time_start,  &
      time_end,    &
      time_units   )
    implicit none

    integer,                    intent(in)  :: step_limit
    integer,                    intent(in)  :: dim_limit
    integer,                    intent(in)  :: fid
    character(len=*),           intent(in)  :: varname
    integer,                    intent(out) :: step_nmax
    character(len=File_HMID),   intent(out) :: description
    character(len=File_HSHORT), intent(out) :: units
    integer,                    intent(out) :: datatype
    integer,                    intent(out) :: dim_rank
    character(len=File_HSHORT), intent(out) :: dim_name  (dim_limit)
    integer,                    intent(out) :: dim_size  (dim_limit)
    real(DP),                   intent(out) :: time_start(step_limit)
    real(DP),                   intent(out) :: time_end  (step_limit)
    character(len=File_HMID),   intent(out) :: time_units

    type(datainfo) :: dinfo

    integer :: ndim
    integer :: istep, idim

    integer :: error
    !---------------------------------------------------------------------------

    ! initialize
    description   = ""
    units         = ""
    datatype      = -1
    dim_rank      = -1
    dim_name(:) = ""
    dim_size(:) = -1
    time_start(:) = rmiss
    time_end(:) = rmiss

    do istep = 1, step_limit
       !--- get data information
       call file_get_datainfo( dinfo,   & ! [OUT]
                               fid,     & ! [IN]
                               varname, & ! [IN]
                               istep,   & ! [IN]
                               .true.,  & ! [IN]
                               error    ) ! [OUT]

       !--- verify and exit
       if ( error /= success_code ) then
          step_nmax = istep - 1
          exit
       endif

       if ( istep == 1 ) then
          description = dinfo%description
          units       = dinfo%units
          datatype    = dinfo%datatype
          dim_rank    = dinfo%rank

          ndim = min( dinfo%rank, dim_limit ) ! limit dimension rank
          do idim = 1, ndim
             dim_name(idim) = dinfo%dim_name(idim)
             dim_size(idim) = dinfo%dim_size(idim)
          enddo

          time_units = dinfo%time_units
       endif

       time_start(istep) = dinfo%time_start
       time_end(istep) = dinfo%time_end
    enddo

    return
  end subroutine filegetalldatainfofid

  !-----------------------------------------------------------------------------
  ! interface File_read
  !-----------------------------------------------------------------------------
  subroutine fileread1drealsp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(SP),    intent(out)           :: var(:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(1)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 1 ) then
       write(message,*) 'xxx rank is not 1', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 1
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:), & ! (out)
         dinfo, sp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread1drealsp
  subroutine fileread1drealdp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(DP),    intent(out)           :: var(:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(1)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 1 ) then
       write(message,*) 'xxx rank is not 1', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 1
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:), & ! (out)
         dinfo, dp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread1drealdp
  subroutine fileread2drealsp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(SP),    intent(out)           :: var(:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(2)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 2 ) then
       write(message,*) 'xxx rank is not 2', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 2
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:), & ! (out)
         dinfo, sp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread2drealsp
  subroutine fileread2drealdp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(DP),    intent(out)           :: var(:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(2)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 2 ) then
       write(message,*) 'xxx rank is not 2', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 2
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:), & ! (out)
         dinfo, dp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread2drealdp
  subroutine fileread3drealsp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(SP),    intent(out)           :: var(:,:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(3)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 3 ) then
       write(message,*) 'xxx rank is not 3', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 3
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:,:), & ! (out)
         dinfo, sp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread3drealsp
  subroutine fileread3drealdp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(DP),    intent(out)           :: var(:,:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(3)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 3 ) then
       write(message,*) 'xxx rank is not 3', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 3
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:,:), & ! (out)
         dinfo, dp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread3drealdp
  subroutine fileread4drealsp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(SP),    intent(out)           :: var(:,:,:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(4)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 4 ) then
       write(message,*) 'xxx rank is not 4', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 4
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:,:,:), & ! (out)
         dinfo, sp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread4drealsp
  subroutine fileread4drealdp( &
      var,           & ! (out)
      basename,      & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      myrank,        & ! (in)
      allow_missing, & ! (in) optional
      single         & ! (in) optional
      )
    implicit none

    real(DP),    intent(out)           :: var(:,:,:,:)
    character(len=*), intent( in)           :: basename
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    integer,          intent( in)           :: myrank
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single

    integer :: fid
    type(datainfo) :: dinfo
    integer :: dim_size(4)
    integer :: error
    integer :: n

    logical :: single_

    intrinsic shape
    !---------------------------------------------------------------------------

    single_ = .false.

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

    !--- search/register file
    call fileopen( fid,                & ! (out)
         basename, file_fread, single_ ) ! (in)

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 4 ) then
       write(message,*) 'xxx rank is not 4', dinfo%rank
       call log('E', message)
    end if
    dim_size(:) = shape(var)
    do n = 1, 4
       if ( dinfo%dim_size(n) /= dim_size(n) ) then
          write(message,*) 'xxx shape is different: ', varname, n, dinfo%dim_size(n), dim_size(n)
          call log('E', message)
       end if
    end do

    call file_read_data( var(:,:,:,:), & ! (out)
         dinfo, dp,                & ! (in)
         error                          ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine fileread4drealdp

  subroutine filereadvar1drealsp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(SP),    intent(out)           :: var(:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 1 ) then
       write(message,*) 'xxx rank is not 1', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:), & ! (out)
            dinfo, sp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar1drealsp
  subroutine filereadvar1drealdp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(DP),    intent(out)           :: var(:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 1 ) then
       write(message,*) 'xxx rank is not 1', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:), & ! (out)
            dinfo, dp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar1drealdp
  subroutine filereadvar2drealsp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(SP),    intent(out)           :: var(:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 2 ) then
       write(message,*) 'xxx rank is not 2', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:), & ! (out)
            dinfo, sp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar2drealsp
  subroutine filereadvar2drealdp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(DP),    intent(out)           :: var(:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 2 ) then
       write(message,*) 'xxx rank is not 2', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:), & ! (out)
            dinfo, dp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar2drealdp
  subroutine filereadvar3drealsp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(SP),    intent(out)           :: var(:,:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 3 ) then
       write(message,*) 'xxx rank is not 3', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:,:), & ! (out)
            dinfo, sp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar3drealsp
  subroutine filereadvar3drealdp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(DP),    intent(out)           :: var(:,:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 3 ) then
       write(message,*) 'xxx rank is not 3', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:,:), & ! (out)
            dinfo, dp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar3drealdp
  subroutine filereadvar4drealsp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(SP),    intent(out)           :: var(:,:,:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:,:) = 0.0_sp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 4 ) then
       write(message,*) 'xxx rank is not 4', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:,:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:,:,:), & ! (out)
            dinfo, sp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar4drealsp
  subroutine filereadvar4drealdp( &
      var,           & ! (out)
      fid,           & ! (in)
      varname,       & ! (in)
      step,          & ! (in)
      allow_missing, & ! (in) optional
      single,        & ! (in) optional
      ntypes,        & ! (in)
      dtype,         & ! (in)
      start,         & ! (in)
      count          & ! (in)
      )
    use mpi, only : mpi_comm_null
    implicit none

    real(DP),    intent(out)           :: var(:,:,:,:)
    integer,          intent( in)           :: fid
    character(len=*), intent( in)           :: varname
    integer,          intent( in)           :: step
    logical,          intent( in), optional :: allow_missing !--- if data is missing, set value to zero
    logical,          intent( in), optional :: single
    integer,          intent( in), optional :: ntypes      ! number of dtypes
    integer,          intent( in), optional :: dtype       ! MPI derived datatype for read buffer
    integer,          intent( in), optional :: start(:)    ! request starts to global variable
    integer,          intent( in), optional :: count(:)    ! request sizes to global variable

    type(datainfo) :: dinfo
    integer :: error

    intrinsic size, shape
    !---------------------------------------------------------------------------

    !--- get data information
    call file_get_datainfo( dinfo,    & ! (out)
         fid, varname, step, .false., & ! (in)
         error                        ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       if ( present(allow_missing) ) then
          if ( allow_missing ) then
             write(message,*) 'xxx [INPUT]/[File] data not found! : ', &
                  'varname= ',trim(varname),', step=',step
             call log('I', message)
             call log('I', 'xxx [INPUT]/[File] Value is set to 0.')
             var(:,:,:,:) = 0.0_dp
          else
             call log('E', 'xxx failed to get data information :'//trim(varname))
          end if
       else
          call log('E', 'xxx failed to get data information :'//trim(varname))
       end if
    end if

    if ( dinfo%rank /= 4 ) then
       write(message,*) 'xxx rank is not 4', dinfo%rank
       call log('E', message)
    end if

    if (present(ntypes) ) then
       call file_read_data_par( var(:,:,:,:),                    & ! (out)
            dinfo, size(shape(var)), ntypes, dtype, start, count, & ! (in)
            error                                                 ) ! (out)
    else
       call file_read_data( var(:,:,:,:), & ! (out)
            dinfo, dp,                & ! (in)
            error                          ) ! (out)
    end if
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data value')
    end if

    return
  end subroutine filereadvar4drealdp

  !-----------------------------------------------------------------------------
  ! interface FileWrite
  !-----------------------------------------------------------------------------
  subroutine filewrite1drealsp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(1)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:), ts, te, sp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:), ts, te, sp, & ! (in)
            1, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite1drealsp
  subroutine filewrite1drealdp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(1)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:), ts, te, dp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:), ts, te, dp, & ! (in)
            1, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite1drealdp
  subroutine filewrite2drealsp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(2)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:), ts, te, sp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:), ts, te, sp, & ! (in)
            2, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite2drealsp
  subroutine filewrite2drealdp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(2)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:), ts, te, dp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:), ts, te, dp, & ! (in)
            2, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite2drealdp
  subroutine filewrite3drealsp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(3)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:,:), ts, te, sp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:,:), ts, te, sp, & ! (in)
            3, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite3drealsp
  subroutine filewrite3drealdp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(3)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:,:), ts, te, dp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:,:), ts, te, dp, & ! (in)
            3, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite3drealdp
  subroutine filewrite4drealsp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:,:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(4)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:,:,:), ts, te, sp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:,:,:), ts, te, sp, & ! (in)
            4, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite4drealsp
  subroutine filewrite4drealdp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end,   & ! (in)
      start,   & ! (in)
      count,   & ! (in)
      ndims    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:,:,:)
    integer,  intent(in)           :: fid
    integer,  intent(in)           :: vid
    real(DP), intent(in)           :: t_start
    real(DP), intent(in)           :: t_end
    integer,  intent(in), optional :: start(:)
    integer,  intent(in), optional :: count(:) ! when var has been reshaped to 1D
    integer,  intent(in), optional :: ndims    ! when var has been reshaped to 1D

    real(DP) :: ts, te

    integer :: start_(4)

    integer :: error, n
    character(len=100) :: str

    intrinsic shape
    !---------------------------------------------------------------------------

    ts = t_start
    te = t_end

    if ( present(ndims) ) then
       ! history variable has been reshaped to 1D
       ! In this case, start and count must be present

       if ( .not. present(start) ) then
          call log('E', 'start argument is neccessary when ndims is specified')
       end if
       if ( .not. present(count) ) then
          call log('E', 'count argument is neccessary when ndims is specified')
       end if

       call file_write_data( fid, vid, var(:,:,:,:), ts, te, dp, & ! (in)
            ndims, start, count,                                       & ! (in)
            error                                                      ) ! (out)
    else
       ! this is for restart variable which keeps its original shape
       if ( present(start) ) then
          start_(:) = start(:)
       else
          start_(:) = 1
       end if
       call file_write_data( fid, vid, var(:,:,:,:), ts, te, dp, & ! (in)
            4, start_, shape(var),                                & ! (in)
            error                                                      ) ! (out)
    end if
    if ( error /= success_code ) then
       do n = 1, file_vid_count
          if ( file_vid_list(n) == vid ) then
             write(str,*) 'xxx failed to write data: ', trim(file_vname_list(n)), mpi_myrank
             exit
          end if
       enddo
       call log('E', trim(str))
    end if

    return
  end subroutine filewrite4drealdp

  !-----------------------------------------------------------------------------
  ! exit netCDF define mode and enter data mode
  subroutine fileenddef( &
       fid & ! (in)
       )
    implicit none

    integer, intent(in) :: fid

    integer :: error, n
    !---------------------------------------------------------------------------

    if ( fid < 0 ) return

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) exit
    end do
    if ( fid .NE. file_fid_list(n) ) then
       write(message,*) 'xxx in FileEndDef invalid fid' , fid
       call log('E', message)
    end if

    call file_enddef( fid, error )

    if ( error == success_code ) then

       call log("I",'')
       write(message,'(A,I3.3,2A)') '###### File end define mode : No.', file_fid_list(n), &
                                                  ', name = ', trim(file_fname_list(n))
       call log("I",message)
       call log("I",'')

    else

       call log('E', 'xxx failed to exit define mode')

    end if

    return
  end subroutine fileenddef

  !-----------------------------------------------------------------------------
  ! This subroutine is used when PnetCDF I/O method is enabled
  subroutine fileattachbuffer( &
       fid,       & ! (in)
       buf_amount ) ! (in)
    implicit none

    integer, intent(in) :: fid
    integer, intent(in) :: buf_amount

    integer :: error, n
    !---------------------------------------------------------------------------

    if ( fid < 0 ) return

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) exit
    end do
    if ( fid .NE. file_fid_list(n) ) then
       write(message,*) 'xxx in FileAttachBuffer invalid fid' , fid
       call log('E', message)
    end if

    call file_attach_buffer( fid, buf_amount, error )

    if ( error == success_code ) then

       write(message,'(A,I3.3,3A,I10)') '###### File attach buffer : No.', file_fid_list(n),  &
                                                       ', name = ', trim(file_fname_list(n)), &
                                                       ', size = ', buf_amount
       call log("I",message)

    else

       call log('E', 'xxx failed to attach buffer in PnetCDF')

    end if

    return
  end subroutine fileattachbuffer

  !-----------------------------------------------------------------------------
  ! This subroutine is used when PnetCDF I/O method is enabled
  subroutine filedetachbuffer( &
       fid        ) ! (in)
    implicit none

    integer, intent(in) :: fid

    integer :: error, n
    !---------------------------------------------------------------------------

    if ( fid < 0 ) return

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) exit
    end do
    if ( n == file_fid_count ) return  ! already closed

    if ( fid .NE. file_fid_list(n) ) then
       write(message,*) 'xxx in FileDetachBuffer invalid fid' , fid
       call log('E', message)
    end if

    call file_detach_buffer( fid, error )

    if ( error == success_code ) then

       write(message,'(A,I3.3,2A)') '###### File detach buffer : No.', file_fid_list(n), &
                                                   ', name = ', trim(file_fname_list(n))
       call log("I",message)

    else

       call log('E', 'xxx failed to detach buffer in PnetCDF')

    end if

    return
  end subroutine filedetachbuffer

  !-----------------------------------------------------------------------------
  ! This subroutine is used when PnetCDF I/O method is enabled
  subroutine fileflush( &
       fid & ! (in)
       )
    implicit none

    integer, intent(in) :: fid

    integer :: error, n
    !---------------------------------------------------------------------------

    if ( fid < 0 ) return
    if ( file_fid_count == 1 ) return

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) exit
    end do
    if ( n == file_fid_count ) return  ! already closed

    if ( fid .NE. file_fid_list(n) ) then
       write(message,*) 'xxx in FileFlush invalid fid' , fid
       call log('E', message)
    end if

    call file_flush( fid, error )

    if ( error == success_code ) then

       write(message,'(A,I3.3,2A)') '###### File flush : No.', file_fid_list(n), &
                                                  ', name = ', trim(file_fname_list(n))
       call log("I",message)
       call log("I",'')

    else

       call log('E', 'xxx failed to flush PnetCDF pending requests')

    end if

    return
  end subroutine fileflush

  !-----------------------------------------------------------------------------
  subroutine fileclose( &
       fid & ! (in)
       )
    implicit none

    integer, intent(in) :: fid

    character(len=File_HLONG) :: fname
    integer                   :: error
    integer                   :: n
    !---------------------------------------------------------------------------

    if ( fid < 0 ) return

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) exit
    end do
    if ( n == file_fid_count ) return  ! already closed

    if ( fid /= file_fid_list(n) ) then
       write(message,*) 'xxx in FileClose invalid fid ', fid
       call log('E', message)
    end if

    call file_close( fid, error )

    if ( error == success_code ) then

       write(message,'(A,I3.3,2A)') '###### File close : No.', file_fid_list(n), &
                                                  ', name = ', trim(file_fname_list(n))
       call log("I",message)
       call log("I",'')

    elseif( error /= already_closed_code ) then

       call log('E', 'xxx failed to close file')

    end if

    do n = 1, file_fid_count-1
       if ( file_fid_list(n) == fid ) then
          file_fid_list(n) = -1
          file_fname_list(n) = ''
       end if
    end do

    return
  end subroutine fileclose
  !-----------------------------------------------------------------------------
  subroutine filecloseall
    implicit none

    integer n
    !---------------------------------------------------------------------------

    do n = 1, file_fid_count-1
       call fileclose( file_fid_list(n) )
    enddo

    return
  end subroutine filecloseall

  !-----------------------------------------------------------------------------
  ! private
  !-----------------------------------------------------------------------------
  subroutine filemakefname( &
       fname,    & ! (out)
       basename, & ! (in)
       prefix,   & ! (in)
       myrank,   & ! (in)
       len       ) ! (in)
    character(len=*), intent(out) :: fname
    character(len=*), intent( in) :: basename
    character(len=*), intent( in) :: prefix
    integer,          intent( in) :: myrank
    integer,          intent( in) :: len

    !                           12345678901234567
    character(len=17) :: fmt = "(A, '.', A, I*.*)"
    !---------------------------------------------------------------------------

    if ( len < 1 .or. len > 9 ) then
       call log('E', 'xxx len is invalid')
    end if

    write(fmt(14:14),'(I1)') len
    write(fmt(16:16),'(I1)') len
    write(fname, fmt) trim(basename), trim(prefix), myrank

    return
  end subroutine filemakefname
  !-----------------------------------------------------------------------------
  subroutine filegetfid( &
      fid,        &
      existed,    &
      basename,   &
      mode,       &
      single,     &
      comm        )
    use mpi, only : mpi_comm_null, mpi_comm_self
    implicit none

    integer,          intent(out) :: fid
    logical,          intent(out) :: existed
    character(len=*), intent( in) :: basename
    integer,          intent( in) :: mode
    logical,          intent( in) :: single
    integer,          intent( in), optional :: comm

    character(len=File_HSHORT) :: rwname(0:2)
    data rwname / 'READ','WRITE','APPEND' /

    character(len=File_HLONG) :: fname
    integer                   :: n

    integer :: error
    integer :: comm_
    !---------------------------------------------------------------------------

    !--- register new file and open
    comm_ = mpi_comm_null
    if ( present(comm) ) comm_ = comm
    if ( comm_ .NE. mpi_comm_null ) then
       ! parallel I/O on a single shared netCDF file
       fname = basename
       comm_ = comm
    elseif ( single ) then
       fname = trim(basename)//'.peall'
    else
       call filemakefname(fname,trim(basename),'pe',mpi_myrank,6)
    endif

    !--- search existing file
    fid = -1
    do n = 1, file_fid_count-1
       if ( fname==file_fname_list(n) ) fid = file_fid_list(n)
    enddo

    if ( fid >= 0 ) then
       existed = .true.
       return
    end if

    call file_open( fid,     & ! (out)
         fname, mode, comm_, & ! (in)
         error               ) ! (out)

    if ( error /= success_code ) then
       call log('E', 'xxx failed to open file :'//trim(fname)//'.nc')
    end if

    call log("I",'')
    write(message,'(3A,I3.3,2A)') '###### File registration (', &
                                  trim(rwname(mode)), ') : No.', fid, ', name = ', trim(fname)
    call log("I",message)

    file_fname_list(file_fid_count) = trim(fname)
    file_fid_list(file_fid_count) = fid
    file_fid_count                  = file_fid_count + 1

    existed = .false.

    return
  end subroutine filegetfid

end module gtool_file
!-------------------------------------------------------------------------------


!--
! vi:set readonly sw=4 ts=8
!
!Local Variables:
!mode: f90
!buffer-read-only: t
!End:
!
!++