gdncvargetnum.f90

Go to the documentation of this file.

! -*- coding: utf-8; mode: f90 -*-
!-------------------------------------------------------------------------------------
! Copyright (c) 2000-2026 Gtool Development Group. All rights reserved.
!-------------------------------------------------------------------------------------
subroutine gdncvargetreal(var, start, cnt, stride, imap, siz, val, iostat)
  use gtdata_netcdf_types,    only: gd_nc_variable
  use gtdata_netcdf_internal, only: gd_nc_variable_entry, vtable_lookup
  use netcdf, only:   &
    & nf90_noerr,     &
    & nf90_einval,    &
    & nf90_eindefine, &
    & nf90_get_var,   &
    & nf90_redef,     &
    & nf90_enddef
  use dc_types,    only: sp
  implicit none
  type(gd_nc_variable), intent(in):: var
  integer, intent(in):: start(:)
  integer, intent(in):: cnt(:)
  integer, intent(in):: stride(:)
  integer, intent(in):: imap(:)
  ! NetCDF変数と内部データ配列のメモリ内構
  ! 造間のマッピングを指定する整数ベクトル.
  ! 詳しくは NetCDF マニュアル
  ! (NF_PUT_VARM_type 等 を参照のこと)
  integer, intent(in):: siz
  real(SP), intent(out):: val(siz)
  integer, intent(out):: iostat
  integer:: nd ! var が保持する変数が依存する次元変数の数
  type(gd_nc_variable_entry):: ent
  integer, allocatable:: istart(:), istride(:), iimap(:)
  continue
  iostat = vtable_lookup(var, ent)
  if (iostat /= nf90_noerr) goto 999
  ! --- nd check ---
  nd = 0
  if (associated(ent%dimids)) nd = size(ent%dimids)
  if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then
    iostat = nf90_einval
    goto 999
  endif
  if (nd == 0) then
    iostat = nf90_get_var(ent%fileid, ent%varid, val(1), start)
    goto 999
  endif
  ! --- stride ovarwrite buffer ---
  allocate(istart(nd), istride(nd), iimap(nd))
  istart(1:nd) = start(1:nd)
  istride(1:nd) = stride(1:nd)
  iimap(1:nd) = imap(1:nd)
  ! --- do read ---
  iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
  if (iostat == nf90_eindefine) then
    iostat = nf90_enddef(ent%fileid)
    if (iostat /= nf90_noerr) return
    iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
    if (iostat /= nf90_noerr) return
    iostat = nf90_redef(ent%fileid)
    if (iostat /= nf90_noerr) return
  end if
  deallocate(istart, istride, iimap)
999 continue
end subroutine gdncvargetreal
 
subroutine gdncvargetdouble(var, start, cnt, stride, imap, siz, val, iostat)
  use gtdata_netcdf_types,    only: gd_nc_variable
  use gtdata_netcdf_internal, only: gd_nc_variable_entry, vtable_lookup
  use netcdf, only:   &
    & nf90_noerr,     &
    & nf90_einval,    &
    & nf90_eindefine, &
    & nf90_get_var,   &
    & nf90_redef,     &
    & nf90_enddef
  use dc_types,    only: dp
  implicit none
  type(gd_nc_variable), intent(in):: var
  integer, intent(in):: start(:)
  integer, intent(in):: cnt(:)
  integer, intent(in):: stride(:)
  integer, intent(in):: imap(:)
  ! NetCDF変数と内部データ配列のメモリ内構
  ! 造間のマッピングを指定する整数ベクトル.
  ! 詳しくは NetCDF マニュアル
  ! (NF_PUT_VARM_type 等 を参照のこと)
  integer, intent(in):: siz
  real(DP), intent(out):: val(siz)
  integer, intent(out):: iostat
  integer:: nd ! var が保持する変数が依存する次元変数の数
  type(gd_nc_variable_entry):: ent
  integer, allocatable:: istart(:), istride(:), iimap(:)
  continue
  iostat = vtable_lookup(var, ent)
  if (iostat /= nf90_noerr) goto 999
  ! --- nd check ---
  nd = 0
  if (associated(ent%dimids)) nd = size(ent%dimids)
  if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then
    iostat = nf90_einval
    goto 999
  endif
  if (nd == 0) then
    iostat = nf90_get_var(ent%fileid, ent%varid, val(1), start)
    goto 999
  endif
  ! --- stride ovarwrite buffer ---
  allocate(istart(nd), istride(nd), iimap(nd))
  istart(1:nd) = start(1:nd)
  istride(1:nd) = stride(1:nd)
  iimap(1:nd) = imap(1:nd)
  ! --- do read ---
  iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
  if (iostat == nf90_eindefine) then
    iostat = nf90_enddef(ent%fileid)
    if (iostat /= nf90_noerr) return
    iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
    if (iostat /= nf90_noerr) return
    iostat = nf90_redef(ent%fileid)
    if (iostat /= nf90_noerr) return
  end if
  deallocate(istart, istride, iimap)
999 continue
end subroutine gdncvargetdouble
 
subroutine gdncvargetint(var, start, cnt, stride, imap, siz, val, iostat)
  use gtdata_netcdf_types,    only: gd_nc_variable
  use gtdata_netcdf_internal, only: gd_nc_variable_entry, vtable_lookup
  use netcdf, only:   &
    & nf90_noerr,     &
    & nf90_einval,    &
    & nf90_eindefine, &
    & nf90_get_var,   &
    & nf90_redef,     &
    & nf90_enddef
  implicit none
  type(gd_nc_variable), intent(in):: var
  integer, intent(in):: start(:)
  integer, intent(in):: cnt(:)
  integer, intent(in):: stride(:)
  integer, intent(in):: imap(:)
  ! NetCDF変数と内部データ配列のメモリ内構
  ! 造間のマッピングを指定する整数ベクトル.
  ! 詳しくは NetCDF マニュアル
  ! (NF_PUT_VARM_type 等 を参照のこと)
  integer, intent(in):: siz
  integer, intent(out):: val(siz)
  integer, intent(out):: iostat
  integer:: nd ! var が保持する変数が依存する次元変数の数
  type(gd_nc_variable_entry):: ent
  integer, allocatable:: istart(:), istride(:), iimap(:)
  continue
  iostat = vtable_lookup(var, ent)
  if (iostat /= nf90_noerr) goto 999
  ! --- nd check ---
  nd = 0
  if (associated(ent%dimids)) nd = size(ent%dimids)
  if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then
    iostat = nf90_einval
    goto 999
  endif
  if (nd == 0) then
    iostat = nf90_get_var(ent%fileid, ent%varid, val(1), start)
    goto 999
  endif
  ! --- stride ovarwrite buffer ---
  allocate(istart(nd), istride(nd), iimap(nd))
  istart(1:nd) = start(1:nd)
  istride(1:nd) = stride(1:nd)
  iimap(1:nd) = imap(1:nd)
  ! --- do read ---
  iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
  if (iostat == nf90_eindefine) then
    iostat = nf90_enddef(ent%fileid)
    if (iostat /= nf90_noerr) return
    iostat = nf90_get_var(ent%fileid, ent%varid, val, istart, cnt, istride, iimap)
    if (iostat /= nf90_noerr) return
    iostat = nf90_redef(ent%fileid)
    if (iostat /= nf90_noerr) return
  end if
  deallocate(istart, istride, iimap)
999 continue
end subroutine gdncvargetint