gtvardeldim.f90 File Reference

Delete dimension. More...

Go to the source code of this file.

Functions/Subroutines
subroutine	gtvardeldim (var, dimord, err)

Detailed Description

Delete 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 gtvardeldim.f90.

Function/Subroutine Documentation

gtvardeldim()

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

Delete dimension from variable

Deletes dimension dimord from variable var. Since this only lowers the rank in the dimension correspondence table and decrements the valid dimension count, I/O is still possible after this operation if the dimension is already degenerate.

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

Note

In NetCDF implementation, variable is not deleted but renamed. This is because netCDF API lacks variable deletion.

Parameters

[in]	var	Variable handle
[in]	dimord	Dimension order number to delete
[out]	err	Error flag

Definition at line 58 of file gtvardeldim.f90.

  use gtdata_types, only: gt_variable
  use gtdata_internal_map, only: map_lookup, gt_dimmap, map_set_ndims, map_set
  use dc_trace, only: beginsub, endsub, dbgmessage
  implicit none
  type(gt_variable), intent(in):: var
  integer, intent(in):: dimord
  logical, intent(out):: err
  type(gt_dimmap), allocatable:: map(:)
  type(gt_dimmap):: tmpmap
  integer:: ndimsp, stat
  character(*), parameter:: subname = 'GTVarDelDim'
continue
  err = .true.
  call beginsub(subname)
  if (dimord < 1) then
    call endsub(subname, "negative dimord=%d invalid", i=(/dimord/))
    return
  endif
  call map_lookup(var, ndims=ndimsp)
  if (ndimsp <= 0) then
    call endsub(subname, "variable invalid")
    return
  else if (dimord > ndimsp) then
    call endsub(subname, "dimord=%d not exist", i=(/dimord/))
    return
  endif
 
  allocate(map(ndimsp))
  call map_lookup(var, map=map)
  tmpmap = map(dimord)
  map(dimord: ndimsp-1) = map(dimord+1: ndimsp)
  map(ndimsp) = tmpmap
  call map_set(var, map, stat)
  deallocate(map)
 
  call map_set_ndims(var, ndims = ndimsp - 1, stat=stat)
  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_set(), and gtdata_internal_map::map_set_ndims().

Here is the call graph for this function: