gdncvarcreate.f90 File Reference

netCDF ファイルへ変数作成 More...

Go to the source code of this file.

Functions/Subroutines
subroutine	gdncvarcreate (var, url, xtype, dims, overwrite, err)

Detailed Description

netCDF ファイルへ変数作成

Author

Eizi TOYODA, Yasuhiro MORIKAWA

Copyright

Copyright (C) GFD Dennou Club, 2000-2026. All rights reserved.
License is BSD-2-Clause. See COPYRIGHT in detail

Definition in file gdncvarcreate.f90.

Function/Subroutine Documentation

gdncvarcreate()

subroutine gdncvarcreate	(	type(gd_nc_variable), intent(out)	var,
		character(len = *), intent(in)	url,
		character(len = *), intent(in)	xtype,
		type(gd_nc_variable), dimension(:), intent(in)	dims,
		logical, intent(in), optional	overwrite,
		logical, intent(out), optional	err
	)

変数を作成

変数 URL url に変数を作成します。 変数が依存する次元を dims に与えます。 返される引数 var には変数 ID などの情報が格納されます。

overwrite に .true. を設定すると上書き可能モード (デフォルト: 上書き不可)。 err を与えた場合、エラー時はプログラム終了せず .false. が返されます。

Parameters

[out]	var	作成した変数ハンドル
[in]	url	変数 URL
[in]	xtype

iline 48 データ型 (float, double, int, char)

Parameters

[in]	dims	次元変数の配列
[in]	overwrite

iline 51 上書きモードを有効化 (省略可能)

Parameters

[out]

err

エラーフラグ (省略可能)

Definition at line 54 of file gdncvarcreate.f90.

  use gtdata_netcdf_types, only: gd_nc_variable, gd_nc_variable_entry, gd_nc_variable_search
  use dc_types, only: string
  use dc_string, only: strieq
  use gtdata_netcdf_internal, only: vtable_add, vtable_lookup
  use gtdata_netcdf_file_generic, only: gdncfileopen, gdncfiledefinemode
  use dc_url, only: urlsplit
  use dc_trace, only: beginsub, endsub, dbgmessage
  use gtdata_netcdf_generic, only: tostring ! for debug
  use netcdf, only: &
    & nf90_noerr, nf90_float, nf90_double, nf90_int, nf90_char, nf90_ebaddim, nf90_def_var
  use dc_error, only: storeerror, gt_enomem, gt_eotherfile, &
    & gt_edimmultidim
  implicit none
  type(GD_NC_VARIABLE), intent(out):: var
  character(len = *), intent(in):: url
  character(len = *), intent(in):: xtype
  type(GD_NC_VARIABLE), intent(in):: dims(:)
  logical, intent(in), optional:: overwrite
  logical, intent(out), optional:: err
  type(GD_NC_VARIABLE_SEARCH):: ent
  type(GD_NC_VARIABLE_ENTRY):: ent_dim
  character(len = string):: filename, varname
  integer, allocatable:: dimids(:)
  integer:: stat, nvdims, i
  integer:: nc_xtype
  logical:: clobber
  intrinsic trim
  character(len = *), parameter:: subnam = "GDNcVarCreate"
continue
  clobber = .false.
  if (present(overwrite)) clobber = overwrite
  call beginsub(subnam)
  call dbgmessage('url=%c', c1=trim(url))
  call dbgmessage('xtype=%c', c1=trim(xtype))
  call dbgmessage('dims=(/%*d/)', i=(/dims(:)%id/), n=(/size(dims)/))
  call dbgmessage('ovwr=%y', l=(/clobber/))
 
  ! もし必要ならファイル作成
  call urlsplit(url, filename, varname)
  call gdncfileopen(ent%fileid, filename, stat=stat, writable=.true., &
    & overwrite=clobber)
  if (stat /= nf90_noerr) goto 999
 
  ! 次元にまつわる準備
  nvdims = size(dims)
  allocate(dimids(max(1, nvdims)), stat=stat)
  if (stat /= 0) then
    stat = gt_enomem
    goto 999
  end if
  do, i = 1, nvdims
    stat = vtable_lookup(dims(i), ent_dim)
    if (stat /= nf90_noerr) then
      stat = nf90_ebaddim
      goto 999
    endif
    if (ent%fileid /= ent_dim%fileid) then
      stat = gt_eotherfile
      goto 999
    endif
    if (ent_dim%dimid <= 0) then
      stat = gt_edimmultidim
      goto 999
    endif
    dimids(i) = ent_dim%dimid
  enddo
  ent%dimid = 0
 
  ! 変数の型の判定
  nc_xtype = nf90_float
  if (strieq(xtype, "double")) then
    nc_xtype = nf90_double
  else if (strieq(xtype, "DOUBLEPRECISION")) then
    nc_xtype = nf90_double
  end if
  if (strieq(xtype, "int")) then
    nc_xtype = nf90_int
  else if (strieq(xtype, "INTEGER")) then
    nc_xtype = nf90_int
  end if
  if (strieq(xtype, "char")) then
    nc_xtype = nf90_char
  else if (strieq(xtype, "CHARACTER")) then
    nc_xtype = nf90_char
  end if
 
  ! 本当の変数作成操作
  stat = gdncfiledefinemode(ent%fileid)
  if (stat /= nf90_noerr) goto 999
  if ( nvdims == 0 ) then
    stat = nf90_def_var(ent%fileid, name = trim(varname), &
      & xtype = nc_xtype, varid=ent%varid)
  else
    stat = nf90_def_var(ent%fileid, name = trim(varname), &
      & xtype = nc_xtype, dimids = dimids, varid=ent%varid)
  end if
  if (stat /= nf90_noerr) goto 999
 
  ! 登録
  stat = vtable_add(var, ent)
 
999 continue
  if (allocated(dimids)) deallocate(dimids)
  if (stat /= nf90_noerr) var % id = -1
  call storeerror(stat, subnam, err, cause_c=url)
  call endsub(subnam, 'stat=%d, var.id=%d', i=(/stat, var % id/))

References dc_trace::beginsub(), dc_trace::dbgmessage(), dc_trace::endsub(), dc_error::gt_edimmultidim, dc_error::gt_enomem, dc_error::gt_eotherfile, dc_error::storeerror(), dc_types::string, gtdata_netcdf_internal::vtable_add(), and gtdata_netcdf_internal::vtable_lookup().

Here is the call graph for this function: