scale_file_h.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_file_h
  use scale_precision, only: &
     dp
  use iso_c_binding
  implicit none
  private
 
#include "scale_file_const.h"
 
  !--- character length
  integer, public, parameter :: file_hshort = file_hshort
  integer, public, parameter :: file_hmid   = file_hmid 
  integer, public, parameter :: file_hlong  = file_hlong
 
  !--- data type
  integer, public, parameter :: file_real4    = file_real4
  integer, public, parameter :: file_real8    = file_real8
  integer, public, parameter :: file_integer2 = file_integer2
  integer, public, parameter :: file_integer4 = file_integer4
  integer, public, parameter :: file_integer8 = file_integer8
  integer, public, parameter :: file_text     = file_text
 
  !--- action type
  integer, public, parameter :: file_fread   = file_fread
  integer, public, parameter :: file_fwrite  = file_fwrite
  integer, public, parameter :: file_fappend = file_fappend
 
  !--- return codes
  integer, public, parameter :: file_error_code           = error_code
  integer, public, parameter :: file_success_code         = success_code
  integer, public, parameter :: file_already_closed_code  = already_closed_code
  integer, public, parameter :: file_already_existed_code = already_existed_code
 
  !--- max
  integer, public, parameter :: file_file_max = file_max
  integer, public, parameter :: file_var_max  = var_max
  integer, public, parameter :: file_rank_max = rank_max
 
  !--- missing value
#define DBL(v) v
  real(dp), public, parameter :: file_rmiss = dbl(rmiss)_dp
 
  !--- struct for data infomation
  type, public, Bind(C) :: datainfo
     character(c_char) :: varname(file_hshort)
     character(c_char) :: description(file_hmid)
     character(c_char) :: units(file_hshort)
     character(c_char) :: standard_name(file_hmid)
     integer           :: datatype
     integer           :: rank
     character(c_char) :: dim_name(file_hshort,rank_max)
     integer           :: dim_size(rank_max)
     integer           :: step
     logical(c_bool)   :: has_tdim
     real(dp)          :: time_start
     real(dp)          :: time_end
     character(c_char) :: time_units(file_hmid)
     character(c_char) :: calendar(file_hshort)
     integer           :: natts
     character(c_char) :: att_name(file_hshort,att_max)
     integer           :: att_type(att_max)
     integer           :: att_len (att_max)
     integer           :: fid
  endtype datainfo
 
  integer, public, parameter :: file_preclist(0:3) = (/ 4, 8, 4, 8 /)
 
  character(len=FILE_HSHORT), public :: file_dtypelist(0:4)
 
  data file_dtypelist / "REAL4", "REAL8", "INTEGER2", "INTEGER4", "INTEGER8" /
 
  public :: cstr
  public :: fstr
 
  interface fstr
     module procedure fstr1
     module procedure fstr2
  end interface fstr
 
contains
  
 
  function cstr(str)
    character(*), intent(in) :: str
    character(:,c_char), allocatable, target :: cstr
    cstr = trim(str) // c_null_char
  end function cstr
 
  subroutine fstr1(str)
    character(len=*), intent(inout) :: str
    integer :: i, j
    do i = 1, len(str)
       if ( str(i:i) == c_null_char ) exit
    end do
    do j = i, len(str)
       str(j:j) = " "
    end do
    return
  end subroutine fstr1
 
  subroutine fstr2(fstr, cstr)
    character(len=*),  intent(out) :: fstr
    character(c_char), intent(in)  :: cstr(:)
    integer :: i, j
    do i = 1, len(fstr)
       if ( cstr(i) == c_null_char ) exit
       fstr(i:i) = cstr(i)
    end do
    do j = i, len(fstr)
       fstr(j:j) = " "
    end do
    return
  end subroutine fstr2
 
 
end module scale_file_h