var/da/da_update_bc/da_get_times_cdf.inc

   1 subroutine da_get_times_cdf(file, times, n_times, max_times, debug)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9 #include "netcdf.inc"
  10
  11    integer,            intent(in)  ::  max_times
  12    integer,            intent(out) ::  n_times
  13    character (len=*),  intent(in)  :: file
  14    character (len=80), intent(out) :: times(max_times)
  15    logical,            intent(in)  :: debug
  16
  17    integer            :: cdfid, rcode, id_time
  18    character (len=80) :: varnam, time1
  19    integer            :: ndims, natts, idims(10), istart(10),iend(10), dimids(10)
  20    integer            :: i, ivtype
  21
  22    ! if (trace_use) call da_trace_entry("da_get_times_cdf")
  23
  24    cdfid = ncopn(file, NCNOWRIT, rcode)
  25
  26    if (rcode == 0) then
  27       if (debug) write(unit=stdout,fmt=*) ' open netcdf file ', trim(file)
  28    else
  29       write(unit=stdout,fmt=*) ' error openiing netcdf file ', trim(file)
  30       stop
  31    end if
  32
  33    id_time = ncvid(cdfid, 'Times', rcode)
  34
  35    rcode = nf_inq_var(cdfid, id_time, varnam, ivtype, ndims, dimids, natts)
  36    if (debug) then
  37       write(unit=stdout,fmt=*) ' number of dims for Time ',ndims
  38    end if
  39    do i=1,ndims
  40       rcode = nf_inq_dimlen(cdfid, dimids(i), idims(i))
  41       if (debug) write(unit=stdout,fmt=*) ' dimension ',i,idims(i)
  42    end do
  43
  44    !  get the times
  45
  46    n_times = idims(2)
  47    do i=1,idims(2)
  48       istart(1) = 1
  49       iend(1) = idims(1)
  50       istart(2) = i
  51       iend(2) = 1
  52
  53       rcode = NF_GET_VARA_TEXT  (cdfid, id_time,  &
  54                                  istart, iend,    &
  55                                  times(i)         )
  56       time1 = times(i)
  57
  58       if (debug) write(unit=stdout,fmt=*) trim(file), time1(1:19)
  59       times(i) = time1(1:19)
  60    end do
  61
  62    write(unit=stdout,fmt=*) ' exiting get_times_cdf '
  63
  64    call ncclos(cdfid,rcode)
  65
  66    ! if (trace_use) call da_trace_exit("da_get_times_cdf")
  67
  68 end subroutine da_get_times_cdf
  69
  70