gtvaradddim.f90 File Reference

Add dimension. More...

Go to the source code of this file.

Functions/Subroutines
subroutine	gtvaradddim (var, dimord, dimvar, err)

Detailed Description

Add dimension.

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 gtvaradddim.f90.

Function/Subroutine Documentation

gtvaradddim()

subroutine gtvaradddim	(	type(gt_variable), intent(in)	var,
		integer, intent(in)	dimord,
		type(gt_variable), intent(in)	dimvar,
		logical, intent(out)	err
	)

Add dimension to variable

Adds dimension dimvar at position dimord in variable var. Dimensions at dimord and after shift one position back. If dimord is greater than the number of valid dimensions of var, it is treated as if (valid_dimensions + 1) was given.

If an error occurs, outputs a message and terminates the program. If err is provided, returns .true. and program does not terminate.

Parameters

[in]	var	Variable handle
[in]	dimord	Dimension order number for insertion
[in]	dimvar	Dimension variable to add
[out]	err	Error flag

Definition at line 55 of file gtvaradddim.f90.

  use gtdata_types, only: gt_variable
  use gtdata_generic, only: inquire
  use gtdata_internal_map, only: map_lookup, gt_dimmap, map_set_ndims, map_set, map_resize
  use dc_trace, only: beginsub, endsub, dbgmessage
  implicit none
  type(GT_VARIABLE), intent(in):: var
  type(GT_VARIABLE), intent(in):: dimvar
  integer, intent(in):: dimord
  logical, intent(out):: err
  type(gt_dimmap), pointer:: map(:)
  type(gt_dimmap):: tmpmap
  integer:: id, nd, ndimsp, stat, vid
  character(*), parameter:: subname = 'GTVarAddDim'
continue
  err = .true.
  call beginsub(subname)
 
  if (dimord < 1) then
    call endsub(subname, "negative dimord=%d invalid", i=(/dimord/))
    return
  endif
 
  ! dimvar をチェックしマップ設定を tmpmap に保存
  call map_lookup(dimvar, vid=vid, ndims=nd)
  if (vid < 0) then
    call endsub(subname, "dimvar invalid")
    return
  endif
  if (nd <= 0) then
    call endsub(subname, "dimvar nondimensional")
    return
  else if (nd > 1) then
    call endsub(subname, "dimvar multidimensional")
    return
  endif
  allocate(map(nd))
  call map_lookup(dimvar, map=map)
  tmpmap = map(1)
  deallocate(map)
 
  ! dimord 番目 (ただし ndimsp + 1 を越えない) に挿入する隙間をあける
  call map_lookup(var, ndims=ndimsp)
  if (dimord > ndimsp + 1) then
    id = ndimsp + 1
  else
    id = dimord
  endif
  allocate(map(nd + 1))
  call map_resize(var, nd + 1)
  call map_lookup(var, map=map)
  map(id+1: nd+1) = map(id: nd)
 
  ! 新しい次元への参照を挿入
  map(id)%dimno = -1
  call inquire(dimvar, url=map(id)%url)
  map(id)%allcount = tmpmap%allcount
  map(id)%offset = tmpmap%offset
  map(id)%step = tmpmap%step
  map(id)%start = tmpmap%start
  map(id)%count = tmpmap%count
  map(id)%stride = tmpmap%stride
 
  ! 登録
  call map_set(var, map=map, stat=stat)
  if (stat /= 0) goto 999
  call map_set_ndims(var, ndims=ndimsp + 1, stat=stat)
 
999 continue
  err = (stat /= 0)
  call endsub(subname)

References dc_trace::beginsub(), dc_trace::dbgmessage(), dc_trace::endsub(), gtdata_internal_map::map_lookup(), gtdata_internal_map::map_resize(), gtdata_internal_map::map_set(), and gtdata_internal_map::map_set_ndims().

Here is the call graph for this function: