da_check_max_iv_gpspw.inc

1 subroutine da_check_max_iv_gpspw(iv, it, ittpw,ittpwf)
2 
3    !-------------------------------------------------------------------------
4    ! Purpose: TBD
5    !-------------------------------------------------------------------------
6 
7    implicit none
8 
9 
10    type(iv_type), intent(inout) :: iv
11    integer,       intent(in)    :: it      ! Outer iteration
12    integer,       intent(inout) :: ittpw,ittpwf
13 
14    logical :: failed 
15    integer :: n
16 
17    if (trace_use_dull) call da_trace_entry("da_check_max_iv_gpspw")
18 
19    if (rootproc .and. check_max_iv_print) then
20       write (check_max_iv_unit,'(A)')  &
21          '----------------------------------------------------------------'
22       write (unit = check_max_iv_unit, fmt = '(A,/)') 'MAX ERROR TEST QC:'
23       write (unit = check_max_iv_unit, fmt = '(/,9(A,F3.0,/))')  &
24          'Error max test ratio for gps_tpw = ',max_error_pw
25    end if
26 
27    !---------------------------------------------------------------------------
28    ! [2.0] Perform maximum innovation vector check:
29    !---------------------------------------------------------------------------
30 
31    failed = .false.
32 
33    ! TPW:
34 
35    do n=iv%info(gpspw)%n1,iv%info(gpspw)%n2
36       call da_max_error_qc(it, iv%info(gpspw), n, iv%gpspw(n)%tpw, max_error_pw, ittpw, &
37         ittpwf, check_max_iv_unit, 'TPW  ', failed, check_max_iv_print)
38    end do
39 
40    if (trace_use_dull) call da_trace_exit("da_check_max_iv_gpspw")
41 
42 end subroutine da_check_max_iv_gpspw
43 
44