da_ao_stats_ssmi_rv.inc

1 subroutine da_ao_stats_ssmi_rv (stats_unit, iv, re)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer,        intent (in)    :: stats_unit    ! Output unit for stats.
10    type (iv_type), intent (inout) :: iv            ! iv
11    type  (y_type), intent (in)    :: re            ! A - O
12 
13    integer                          :: ntpw, nspeed, n
14    type (stats_ssmi_rv_type) :: stats
15 
16    if (trace_use) call da_trace_entry("da_ao_stats_ssmi_rv")
17 
18    ntpw = 0
19    nspeed = 0
20 
21    stats%maximum%tpw   = maxmin_type (missing_r, 0, 0)
22    stats%maximum%speed = maxmin_type (missing_r, 0, 0)
23    stats%minimum%tpw   = maxmin_type(-missing_r, 0, 0)
24    stats%minimum%speed = maxmin_type(-missing_r, 0, 0)
25    stats%average = residual_ssmi_rv_type(0.0, 0.0)
26    stats%rms_err = stats%average
27 
28    if (iv%info(ssmi_rv)%nlocal > 0) then
29       do n=1, iv%info(ssmi_rv)%nlocal
30          if (iv%info(ssmi_rv)%proc_domain(1,n)) then
31             call da_stats_calculate (n, 0, iv%ssmi_rv(n)%tpw%qc, &
32                re%ssmi_rv(n)%tpw, ntpw,  &
33                stats%minimum%tpw,  stats%maximum%tpw,&
34                stats%average%tpw,  stats%rms_err%tpw)
35 
36             call da_stats_calculate (n, 0, iv%ssmi_rv(n)%speed%qc, &
37                re%ssmi_rv(n)%speed,    nspeed,  &
38                stats%minimum%speed,  stats%maximum%speed,&
39                stats%average%speed,  stats%rms_err%speed)
40          end if
41       end do   
42    end if
43 
44    ! Do inter-processor communication to gather statistics.
45    call da_proc_sum_int (ntpw)
46    call da_proc_sum_int (nspeed)
47    iv%nstats(ssmi_rv) = ntpw + nspeed
48    
49    call da_proc_stats_combine(stats%average%tpw, stats%rms_err%tpw, &
50       stats%minimum%tpw%value, stats%maximum%tpw%value, &
51       stats%minimum%tpw%n, stats%maximum%tpw%n, &
52       stats%minimum%tpw%l, stats%maximum%tpw%l)
53 
54    call da_proc_stats_combine(stats%average%speed, stats%rms_err%speed, &
55       stats%minimum%speed%value, stats%maximum%speed%value, &
56       stats%minimum%speed%n, stats%maximum%speed%n, &
57       stats%minimum%speed%l, stats%maximum%speed%l)
58 
59    if (rootproc) then
60       if (ntpw /= 0) then
61 
62          write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of AO for ssmi_retrieval'
63          write(unit=stats_unit, fmt='(a/)') '   var           tpw(cm)     n'
64          write(unit=stats_unit, fmt='(a,i14)') '  Number: ', ntpw
65          write(unit=stats_unit, fmt='(a, f12.4,i5)') &
66             ' Minimum(n): ', stats%minimum%tpw%value, stats%minimum%tpw%n    , &
67             ' Maximum(n): ', stats%maximum%tpw%value, stats%maximum%tpw%n
68          write(unit=stats_unit, fmt='(a, f12.4,5x)') &
69             ' Average   : ', stats%average%tpw/real(ntpw), &
70             '    RMSE   : ', sqrt(stats%rms_err%tpw/real(ntpw))
71       end if
72 
73       if (nspeed /= 0) then
74          write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of AO for ssmi_retrieval'
75          write(unit=stats_unit, fmt='(a/)') '   var           Speed(m/s)     n'
76          write(unit=stats_unit, fmt='(a,i14)') '  Number: ', nspeed
77          write(unit=stats_unit, fmt='(a, f12.4,i5)') &
78             ' Minimum(n): ', stats%minimum%speed%value, stats%minimum%speed%n    , &
79             ' Maximum(n): ', stats%maximum%speed%value, stats%maximum%speed%n
80          write(unit=stats_unit, fmt='(a, f12.4,5x)') &
81             ' Average   : ', stats%average%Speed/real(nspeed), &
82             '    RMSE   : ', sqrt(stats%rms_err%Speed/real(nspeed))
83       end if
84    end if
85 
86    if (trace_use) call da_trace_exit("da_ao_stats_ssmi_rv")
87 
88 end subroutine da_ao_stats_ssmi_rv
89 
90