da_oi_stats_satem.inc

1 subroutine da_oi_stats_satem (stats_unit, oi)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer,        intent (in)      :: stats_unit    ! Output unit for stats.
10    type (ob_type), intent (in)      :: oi            ! OI
11 
12    type (stats_satem_type)          :: Satem
13    integer                          :: nthickness
14    integer                          :: n, k
15 
16    nthickness = 0
17 
18    Satem%maximum%thickness = maxmin_type(missing_r, 0, 0)
19    Satem%minimum%thickness = maxmin_type(-missing_r, 0, 0)
20    Satem%average = residual_satem1_type(0.0)
21    Satem%rms_err = Satem%average
22 
23    if (oi%num_Satem > 0) then
24       do n=1, oi%num_Satem
25          if (oi%satem(n)%loc%proc_domain) then
26             do k=1, oi%satem(n)%info%levels
27                call da_stats_calculate(oi%Satem(n)%loc%obs_global_index, &
28                           k, oi%Satem(n)%thickness(k)%qc, &
29                           oi%Satem(n)%thickness(k)%inv, nthickness, &
30                           Satem%minimum%thickness, Satem%maximum%thickness, &
31                           Satem%average%thickness, Satem%rms_err%thickness)
32             end do
33          end if    ! end if (oi%satem(n)%loc%proc_domain)
34       end do
35    end if
36 
37    ! Do inter-processor communication to gather statistics.
38    call da_proc_sum_int(nthickness)
39 
40    call da_proc_stats_combine(Satem%average%thickness, Satem%rms_err%thickness, &
41              Satem%minimum%thickness%value, Satem%maximum%thickness%value, &
42              Satem%minimum%thickness%n, Satem%maximum%thickness%n, &
43              Satem%minimum%thickness%l, Satem%maximum%thickness%l)
44 
45    if (rootproc) then
46       if (nthickness /= 0) then
47          write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of OI for satem'
48          call da_print_stats_satem(stats_unit, nthickness, Satem)
49       end if
50    end if
51 
52 end subroutine da_oi_stats_satem
53 
54