var/da/da_ssmi/da_check_max_iv_ssmt1.inc

   1 subroutine da_check_max_iv_ssmt1(iv, it, num_qcstat_conv)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    ! Update:
   6    !    Removed Outerloop check as it is done in da_get_innov
   7    !    Author: Syed RH Rizvi,  MMM/NESL/NCAR,  Date: 07/12/2009
   8    !-----------------------------------------------------------------------
   9
  10    implicit none
  11
  12    type(iv_type), intent(inout) :: iv
  13    integer,       intent(in)    ::  it      ! outer iteration
  14    integer,       intent(inout) :: num_qcstat_conv(:,:,:,:)
  15
  16
  17
  18    integer ::  k,n, ipr
  19    logical :: failed
  20
  21    if (trace_use_dull) call da_trace_entry("da_check_max_iv_ssmt1")
  22
  23    !---------------------------------------------------------------------------
  24    ! [1.0] Perform maximum innovation vector check:
  25    !---------------------------------------------------------------------------
  26
  27    do n=iv%info(ssmt1)%n1,iv%info(ssmt1)%n2
  28       do k = 1, iv%info(ssmt1)%levels(n)
  29          call da_get_print_lvl(iv%ssmt1(n)%p(k),ipr)
  30
  31         failed= .false.
  32         if( iv%ssmt1(n)%t(k)%qc >= obs_qc_pointer ) then
  33         call da_max_error_qc(it, iv%info(ssmt1), n, iv%ssmt1(n)%t(k), max_error_t, failed)
  34         if( iv%info(ssmt1)%proc_domain(k,n) ) then
  35          num_qcstat_conv(1,ssmt1,3,ipr) = num_qcstat_conv(1,ssmt1,3,ipr) + 1
  36         if(failed) then
  37         num_qcstat_conv(2,ssmt1,3,ipr) = num_qcstat_conv(2,ssmt1,3,ipr) + 1
  38         if ( write_rej_obs_conv ) then
  39         write(qcstat_conv_unit,'(2x,a10,2x,a4,3f12.2)')&
  40         'ssmt1',ob_vars(3),iv%info(ssmt1)%lat(k,n),iv%info(ssmt1)%lon(k,n),0.01*iv%ssmt1(n)%p(k)
  41         endif
  42         endif
  43         endif
  44         endif
  45
  46       end do
  47    end do
  48
  49    if (trace_use_dull) call da_trace_exit("da_check_max_iv_ssmt1")
  50
  51 end subroutine da_check_max_iv_ssmt1
  52
  53