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 :: filewritevar
  public :: filegetglobalattribute
  public :: filesetglobalattribute
  public :: fileenddef
  public :: fileclose
  public :: filecloseall
  public :: filemakefname

  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
  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 filewritevar
    module procedure filewritevar1drealsp
    module procedure filewritevar1drealdp
    module procedure filewritevar2drealsp
    module procedure filewritevar2drealdp
    module procedure filewritevar3drealsp
    module procedure filewritevar3drealdp
    module procedure filewritevar4drealsp
    module procedure filewritevar4drealdp
  end interface filewritevar
  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
    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

    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)
         )

    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)
    call filesetglobalattribute( fid, & ! (in)
         "myrank", (/myrank/)         ) ! (in)
    call filesetglobalattribute( fid, & ! (in)
         "rankidx", rankidx           ) ! (in)

    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,                & ! (in)
                          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
      )
    implicit none

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

    logical :: existed
    logical :: single_ = .false.

    if ( present(single) ) single_ = single

    call filegetfid( fid,        & ! (out)
         existed,                & ! (out)
         basename, mode, single_ ) ! (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,                                          & ! (in)
         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,                                          & ! (in)
         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)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(SP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_write_axis( fid, name, val, sp, & ! (in)
         error                                     ) ! (out)
    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)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(DP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_write_axis( fid, name, val, dp, & ! (in)
         error                                     ) ! (out)
    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)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(SP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite1dassociatedcoordinatesrealsp
  subroutine filewrite1dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(DP),    intent(in) :: val(:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite1dassociatedcoordinatesrealdp
  subroutine filewrite2dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(SP),    intent(in) :: val(:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite2dassociatedcoordinatesrealsp
  subroutine filewrite2dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(DP),    intent(in) :: val(:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite2dassociatedcoordinatesrealdp
  subroutine filewrite3dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(SP),    intent(in) :: val(:,:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite3dassociatedcoordinatesrealsp
  subroutine filewrite3dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(DP),    intent(in) :: val(:,:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite3dassociatedcoordinatesrealdp
  subroutine filewrite4dassociatedcoordinatesrealsp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(SP),    intent(in) :: val(:,:,:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 filewrite4dassociatedcoordinatesrealsp
  subroutine filewrite4dassociatedcoordinatesrealdp( &
       fid,       & ! (in)
       name,      & ! (in)
       val        ) ! (in)
    integer,          intent(in) :: fid
    character(len=*), intent(in) :: name
    real(DP),    intent(in) :: val(:,:,:,:)

    integer error
    intrinsic size

    call file_write_associated_coordinates( fid, name, 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 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,*) '*** [File] Var registration'
       call log("I", message)
       write(message,*) '*** variable 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,*) '*** [File] Var registration'
       call log("I", message)
       write(message,*) '*** variable 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,*) '*** [File] Var registration'
       call log("I", message)
       write(message,*) '*** variable name: ', trim(varname)
       call log("I", message)

       if ( .NOT. present(tint) ) then
          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
      )
    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

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

    logical :: single_ = .false.

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

    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, 1, .false., & ! (in)
         error                     ) ! (out)

    !--- verify
    if ( error /= success_code ) then
       call log('E', 'xxx failed to get data information :'//trim(varname))
    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

    return
  end subroutine filegetshape

  !-----------------------------------------------------------------------------
  ! FileGetData
  !-----------------------------------------------------------------------------
  subroutine filegetdatainfo( &
      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

    integer        :: fid
    type(datainfo) :: dinfo

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

    integer :: error
    logical :: single_ = .false.

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

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

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

    !--- 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 in '//trim(basename))
    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 filegetdatainfo

  !-----------------------------------------------------------------------------
  ! FileGetData
  !-----------------------------------------------------------------------------
  subroutine filegetalldatainfo( &
      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
    type(datainfo) :: dinfo

    integer :: ndim
    integer :: istep, idim
    logical :: flag_first = .true.

    integer :: error
    logical :: single_ = .false.
    !---------------------------------------------------------------------------

    mpi_myrank = myrank

    if ( present(single) ) single_ = single

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

    ! 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 ( flag_first ) then
          flag_first = .false.

          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 filegetalldatainfo

  !-----------------------------------------------------------------------------
  ! 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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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_ = .false.

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

    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

  !-----------------------------------------------------------------------------
  ! interface FileWrite
  !-----------------------------------------------------------------------------
  subroutine filewrite1drealsp( &
      fid,     & ! (in)
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:,:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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)
      )
    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

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_data( fid, vid, var(:,:,:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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

  !-----------------------------------------------------------------------------
  ! interface FileWriteVar
  !-----------------------------------------------------------------------------
  subroutine filewritevar1drealsp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar1drealsp
  subroutine filewritevar1drealdp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar1drealdp
  subroutine filewritevar2drealsp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar2drealsp
  subroutine filewritevar2drealdp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar2drealdp
  subroutine filewritevar3drealsp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar3drealsp
  subroutine filewritevar3drealdp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar3drealdp
  subroutine filewritevar4drealsp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(SP), intent(in) :: var(:,:,:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:,:,:), ts, te, sp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar4drealsp
  subroutine filewritevar4drealdp( &
      vid,     & ! (in)
      var,     & ! (in)
      t_start, & ! (in)
      t_end    & ! (in)
      )
    implicit none

    real(DP), intent(in) :: var(:,:,:,:)
    integer,  intent(in) :: vid
    real(DP), intent(in) :: t_start
    real(DP), intent(in) :: t_end

    real(DP) :: ts, te

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

    ts = t_start
    te = t_end
    call file_write_var( vid, var(:,:,:,:), ts, te, dp, & ! (in)
         error                                     ) ! (out)
    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 filewritevar4drealdp

  !-----------------------------------------------------------------------------
  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 , & ! (in)
         error             ) ! (out)
    if ( error .EQ. success_code ) then
       write(message, '(1x,A,i3)') '*** [File] File enddef : NO.', n
       call log('I', message)
       call log('I', '*** enddef filename: ' // trim(file_fname_list(n)))
    else
       call log('E', 'xxx failed to exit define mode')
    end if

    return
  end subroutine fileenddef

  !-----------------------------------------------------------------------------
  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 .EQ. 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 , & ! (in)
         error             ) ! (out)
    if ( error == success_code ) then
       write(message, '(1x,A,i3)') '*** [File] File Close : NO.', n
       call log('I', message)
       call log('I', '*** closed filename: ' // trim(file_fname_list(n)))
    else if ( 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      )
    implicit none

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


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

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

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

    !--- register new file and open
    if ( 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,    & ! (in)
         error           ) ! (out)
    if ( error /= success_code ) then
       call log('E', 'xxx failed to open file :'//trim(fname)//'.nc')
    end if

    write(message,*) '*** [File] File registration : ',trim(rwname(mode)),' -', fid
    call log("I", message)
    write(message,*) '*** filename: ', 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:
!
!++