var/da/da_update_bc/da_get_bdytimestr_cdf.inc

   1 subroutine da_get_bdytimestr_cdf(file, time_flag, bdytimestr, max_times, debug)
   2
   3    implicit none
   4
   5 #include "netcdf.inc"
   6
   7    integer,           intent(in) :: max_times
   8    logical,           intent(in) :: debug
   9    character(len=*),  intent(in) :: file
  10    character(len=*),  intent(in) :: time_flag
  11    character(len=80), intent(out) :: bdytimestr(max_times)
  12
  13    character(len=80) :: varnam, time1
  14
  15    integer :: cdfid, rcode, id_time
  16    integer :: ndims, natts, idims(10), istart(10),iend(10), dimids(10)
  17    integer :: i, ivtype
  18
  19    ! if (trace_use) call da_trace_entry("da_get_bdytimestr_cdf")
  20
  21    cdfid = ncopn(file, NCNOWRIT, rcode)
  22
  23    if (rcode == 0) then
  24       if (debug) write(unit=stdout,fmt=*) ' open netcdf file ', trim(file)
  25    else
  26       write(unit=stdout,fmt=*) ' error openiing netcdf file ', trim(file)
  27       stop
  28    end if
  29
  30    varnam = 'md___' // trim(time_flag) // 'e_x_t_d_o_m_a_i_n_m_e_t_a_data_'
  31
  32    id_time = ncvid(cdfid, varnam, rcode)
  33
  34    rcode = nf_inq_var(cdfid, id_time, varnam, ivtype, ndims, dimids, natts)
  35
  36    if (debug) then
  37      write(unit=stdout,fmt=*) ' number of dims for bdytime ',ndims
  38    end if
  39
  40    do i=1,ndims
  41       rcode = nf_inq_dimlen(cdfid, dimids(i), idims(i))
  42       if (debug) write(unit=stdout,fmt=*) ' dimension ',i,idims(i)
  43    end do
  44
  45    ! get the bdytime
  46
  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                                  bdytimestr(i))
  56       time1 = bdytimestr(i)
  57
  58       if (debug) write(unit=stdout,fmt=*) trim(file), time1(1:19)
  59       bdytimestr(i) = time1(1:19)
  60    end do
  61
  62    call ncclos(cdfid,rcode)
  63
  64    ! if (trace_use) call da_trace_exit("da_get_bdytimestr_cdf")
  65
  66 end subroutine da_get_bdytimestr_cdf
  67
  68