gdncvarcreate.f90 File Reference

Create variable in netCDF file. More...

Go to the source code of this file.

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

Detailed Description

Create variable in netCDF file.

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

These subroutines and functions are provided as gtdata_netcdf_generic#Create through gtdata_netcdf_generic.

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 )

Create a variable

Creates a variable at the specified URL. The dimensions the variable depends on are given in dims. The returned var contains the variable ID and other information.

Setting overwrite to .true. enables overwrite mode (default: no overwrite). If err is provided, errors return .false. instead of terminating.

Parameters

[out]	var	Created variable handle
[in]	url	Variable URL
[in]	xtype	Data type (float, double, int, char)
[in]	dims	Array of dimension variables
[in]	overwrite	Enable overwrite mode (optional)
[out]	err	Error flag (optional)

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: