gtvarsliceauto.f90

Go to the documentation of this file.

 
 
subroutine gtvarsliceauto(var, compatible)
  use gtdata_types, only: gt_variable
  use gtdata_generic, only: inquire, get_slice, gtvarslice
  type(gt_variable), intent(inout):: var
  type(gt_variable), intent(in), optional:: compatible
  integer:: nd, i
  integer, allocatable:: start(:), count(:), stride(:)
  call inquire(var, alldims=nd)
  if (nd <= 0) return
  allocate(start(nd), count(nd), stride(nd))
  if (present(compatible)) then
    call get_slice(compatible, start, count, stride)
    do, i = 1, nd
      call gtvarslice(var, i, start(i), count(i), stride(i))
    enddo
  else
    call get_slice(var, start, count, stride)
    do, i = 1, nd
      if (count(i) < 1) count(i) = 1
      call gtvarslice(var, i, start(i), count(i), stride(i))
    enddo
  endif
  deallocate(start, count, stride)
end subroutine gtvarsliceauto