var/da/da_radar/da_max_error_qc_radar.inc

   1 subroutine da_max_error_qc_radar (it, info, n,field, max_error, ix, ixf, iunit, var, failed, print_details)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    integer          ,   intent (in)   :: it
  10    type (infa_type) ,   intent(in)    :: info
  11    integer,             intent(in)    :: n
  12    type (field_type),   intent(inout) :: field
  13    real             ,   intent(in)    :: max_error
  14    integer          ,   intent(inout) :: ix, ixf
  15    integer          ,   intent(in)    :: iunit
  16    character (len=*),   intent(in)    :: var
  17    logical          ,   intent(out)   :: failed
  18    logical          ,   intent(in)    :: print_details
  19
  20    real                               :: err, err_max
  21    integer                            :: qc_flag
  22
  23    if (trace_use) call da_trace_entry("da_max_error_qc_radar")
  24
  25    failed = .false.
  26
  27    qc_flag = field % qc
  28    err_max = field % error * max_error
  29    err     = field % inv
  30    err     = ABS (err)
  31
  32    ix     = ix + 1
  33    if (it > 1 .and. qc_flag == fails_error_max) field%qc = 0
  34    if (err > err_max) then
  35       if (field % qc > fails_error_max) field % qc = fails_error_max
  36
  37       ixf = ixf + 1
  38       failed = .true.
  39
  40       if (print_details .and. failed) then
  41          if (err_max .LE. 0.0) then
  42             write (iunit , fmt = '(A,3(F12.1,1X),A,A,A,A,A,3f10.2)')   &
  43                "Err_max < 0 ==> ",err,err_max,max_error, " for ", var, &
  44                " OBS ID: ", info%platform(n),     &
  45                " LA/LON/ELV:", info%lat(1,n), info%lon(1,n), info%elv(n)
  46             ! call da_error(__FILE__,__LINE__,(/"Erk"/))
  47          end if
  48
  49          write (iunit , fmt = '(A,A,A,I5,A,I5,A,F4.1,A,A,A,2F12.1)') &
  50             "Err_max failed:ID=", info%platform(n),&
  51             "  Ix=", ix, "  Ixf=", ixf, " Err_max ratio =",err/err_max, &
  52             " for ", var, " inv, error:",field % inv, field % error
  53       end if
  54       field % inv = 0.0
  55    end if
  56
  57    if (trace_use) call da_trace_exit("da_max_error_qc_radar")
  58
  59 end subroutine da_max_error_qc_radar
  60
  61