var/da/da_obs/da_obs_proc_station.inc

   1 subroutine da_obs_proc_station(obs,fm,uvq_direct)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: Processing obs data read in at a station:
   5    !
   6    ! (1) Surface correction
   7    ! (2) Vertical checks
   8    ! (3) Missing data check
   9    !-----------------------------------------------------------------------
  10
  11    implicit none
  12
  13    type (multi_level_type), intent(inout) :: obs
  14
  15    logical,optional :: uvq_direct
  16    real    :: po, to, rho, es, qs, qvo
  17    integer :: i,fm
  18
  19    if (trace_use_dull) call da_trace_entry("da_obs_proc_station")
  20
  21    do i = 1, obs % info % levels
  22       !-----------------------------------------------------------------------
  23       ! Gross check for t, p & q
  24       !-----------------------------------------------------------------------
  25
  26       if (obs%each(i)%t %inv <  75.0 .or. obs%each(i)%t %inv > 350.0  ) then
  27          obs%each(i)%t %inv = missing_r
  28          obs%each(i)%t %qc  = missing
  29       end if
  30       if (obs%each(i)%p %inv <=  0.0 .or. obs%each(i)%p %inv > 110000.0) then
  31          obs%each(i)%p %inv = missing_r
  32          obs%each(i)%p %qc  = missing
  33       end if
  34       if (obs%each(i)%rh%inv < 0.0)  then
  35          obs%each(i)%rh%inv = missing_r
  36          obs%each(i)%rh%qc  = missing
  37       end if
  38
  39      if(.not. (present(uvq_direct) .and. uvq_direct) .and. fm.ne.161) then
  40       po   = obs % each(i) % p  % inv
  41       to   = obs % each(i) % t  % inv
  42       rho  = obs % each(i) % rh  % inv
  43
  44       if (ob_format == ob_format_ascii) then ! Calculate q if possible:
  45          if (abs(po -  missing_r) > 1.0 .and. abs(to -  missing_r) > 1.0 .and. &
  46              abs(rho-  missing_r) > 1.0) then
  47
  48             call da_tp_to_qs(to, po, es, qs)
  49
  50             if (rho > 100.0) then
  51                qvo = qs
  52             else
  53                qvo  = rho * qs / 100.0
  54             end if
  55          else
  56             qvo       = missing_r
  57          end if
  58
  59          obs % each(i) % q  % inv = qvo
  60          obs % each(i) % q  % qc = obs % each(i) % rh % qc
  61          obs % each(i) % q  % error = obs % each(i) % rh % error
  62       end if
  63     else
  64          obs % each(i) % q  % inv = obs % each(i) % rh  % inv
  65          obs % each(i) % q  % qc = obs % each(i) % rh % qc
  66          obs % each(i) % q  % error = obs % each(i) % rh % error
  67     end if
  68    end do
  69    if (trace_use_dull) call da_trace_exit("da_obs_proc_station")
  70
  71 end subroutine da_obs_proc_station
  72
  73