da_check_max_iv_ssmi_rv.inc
References to this file elsewhere.
1 subroutine da_check_max_iv_ssmi_rv(iv, it, itpw, itpwf, ispeed, ispeedf)
2
3 !-----------------------------------------------------------------------
4 ! Purpose: TBD
5 !-----------------------------------------------------------------------
6
7 implicit none
8
9 type(iv_type), intent(inout) :: iv
10 integer, intent(in) :: it ! Outer loop
11 integer, intent(inout) :: itpw, itpwf, ispeed, ispeedf
12
13 logical :: failed
14 integer :: n
15
16 if (trace_use) call da_trace_entry("da_check_max_iv_ssmi_rv")
17
18 !---------------------------------------------------------------------------
19 ! [1.0] Open diagnostic file:
20 !---------------------------------------------------------------------------
21
22 if (rootproc .and. check_max_iv_print) then
23 write (check_max_iv_unit,'(A)') &
24 '----------------------------------------------------------------'
25 write (unit = check_max_iv_unit, fmt = '(A,/)') 'MAX ERROR TEST QC:'
26
27 write (unit = check_max_iv_unit, fmt = '(/,9(A,F3.0,/))') &
28 'Error max test ratio for SSMI_rv_TPW = ',max_error_pw, &
29 'Error max test ratio for SSMI_rv_SPEED = ',max_error_uv
30 end if
31
32 !---------------------------------------------------------------------------
33 ! [2.0] Perform maximum innovation vector check:
34 !---------------------------------------------------------------------------
35
36 failed = .false.
37
38 do n=iv%info(ssmi_rv)%n1,iv%info(ssmi_rv)%n2
39 call da_max_error_qc (it, iv%info(ssmi_rv), n, iv%ssmi_rv(n)%tpw, &
40 max_error_pw, itpw, itpwf, check_max_iv_unit, 'TPW ', failed, check_max_iv_print)
41 call da_max_error_qc (it, iv%info(ssmi_rv), n, iv%ssmi_rv(n)%speed, &
42 max_error_uv, ispeed, ispeedf, check_max_iv_unit, 'SPEED', failed, check_max_iv_print)
43 end do
44
45 if (trace_use) call da_trace_exit("da_check_max_iv_ssmi_rv")
46
47 end subroutine da_check_max_iv_ssmi_rv
48
49