gtvaropen.f90

Go to the documentation of this file.

 
 
subroutine gtvaropen(var, url, writable, err)
  use dc_string, only: strhead
  use dc_types, only: string
  use gtdata_types, only: gt_variable
  use gtdata_generic, only: limit
  use gtdata_internal_map, only: map_create, vtb_class_netcdf, gtvar_dump
  use gtdata_netcdf_generic, only: open, inquire
  use gtdata_netcdf_types, only: gd_nc_variable
  use dc_url, only: url_chop_iorange
  use dc_error, only: storeerror, dc_noerr, gt_enotvar
  use dc_trace, only: beginsub, endsub
  use dc_present, only: present_and_true
  implicit none
  type(gt_variable), intent(out):: var
  character(*), intent(in):: url
  logical, intent(in), optional:: writable
  logical, intent(out), optional:: err
  integer:: ndims, stat, cause_i
  character(STRING):: cause_c
  integer, allocatable:: dimlen(:)
  type(gd_nc_variable):: gdnc
  character(STRING):: filevar, iorange
  character(*),      parameter:: subname = "GTVarOpen"
  character(*),      parameter:: version = &
    & '$Name:  $' // &
    & '$Id: gtvaropen.f90,v 1.4 2009-05-25 09:55:57 morikawa Exp $'
continue
  call beginsub(subname, fmt='<%c>', c1=trim(url), version=version)
  stat = dc_noerr
  cause_i = 0
  cause_c = ''
  var = gt_variable(-1)
  call url_chop_iorange(url, iorange=iorange, remainder=filevar)
  call open(gdnc, filevar, writable, err)
  if ( present_and_true(err) ) then
    stat = gt_enotvar
    goto 999
  end if
  call inquire(gdnc, ndims=ndims)
  allocate(dimlen(max(1, ndims)))
  call inquire(gdnc, dimlen=dimlen)
  call map_create(var, vtb_class_netcdf, gdnc%id, ndims, dimlen, stat)
  if (stat /= dc_noerr) then
    cause_i = ndims
    goto 999
  end if
  deallocate(dimlen)
  call limit(var, trim(iorange))
  call gtvar_dump(var)
999 continue
  call storeerror(stat, subname, err, cause_c = cause_c, cause_i = cause_i)
  call endsub(subname, 'mapid=%d', i=(/var%mapid/))
end subroutine gtvaropen