gdncsearchdim.f90

Go to the documentation of this file.

 
 
integer function gdncsearchdim(var, dimname) result(result)
  use gtdata_netcdf_types, only: gd_nc_variable, gd_nc_variable_entry, gd_nc_variable_search
  use gtdata_netcdf_internal, only: vtable_lookup
  use gtdata_netcdf_file_generic
  use netcdf, only: nf90_ebaddim, nf90_noerr, nf90_inq_dimid
  use dc_error
  use dc_trace, only: beginsub, endsub
  implicit none
  type(gd_nc_variable), intent(in):: var
  character(*), intent(in):: dimname
  type(gd_nc_variable_entry):: ent
  integer:: stat, dimid, i
  character(*), parameter:: subname = 'GDNcSearchDim'
continue
  call beginsub(subname, 'var=%d dimname=%c', i=(/var%id/), c1=trim(dimname))
  result = nf90_ebaddim
  stat = vtable_lookup(var, ent)
  if (stat /= nf90_noerr) goto 999
 
  stat = nf90_inq_dimid(ent%fileid, dimname, dimid)
  if (stat /= nf90_noerr) goto 999
 
  if (dimid == ent%dimid) then
    result = 1
    goto 999
  endif
 
  if (.not. associated(ent%dimids)) then
    stat = gt_enomoredims
    goto 999
  endif
  do, i = 1, size(ent%dimids)
    if (ent%dimids(i) == dimid) then
      result = i
      goto 999
    endif
  enddo
 
999 continue
  call endsub(subname, 'result=%d', i=(/result/))
  return
end function gdncsearchdim