da_ao_stats_sound.inc

1 subroutine da_ao_stats_sound (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    type (stats_sound_type) :: stats
14    integer                 :: nu, nv, nt, nq
15    integer                 :: n, k
16 
17    if (trace_use_dull) call da_trace_entry("da_ao_stats_sound")
18 
19    nu = 0
20    nv = 0
21    nt = 0
22    nq = 0
23 
24    stats%maximum%u = maxmin_type (missing_r, 0, 0)
25    stats%maximum%v = maxmin_type (missing_r, 0, 0)
26    stats%maximum%t = maxmin_type (missing_r, 0, 0)
27    stats%maximum%q = maxmin_type (missing_r, 0, 0)
28    stats%minimum%u = maxmin_type(-missing_r, 0, 0)
29    stats%minimum%v = maxmin_type(-missing_r, 0, 0)
30    stats%minimum%t = maxmin_type(-missing_r, 0, 0)
31    stats%minimum%q = maxmin_type(-missing_r, 0, 0)
32 
33    stats%average = residual_sound1_type(0.0, 0.0, 0.0, 0.0)
34    stats%rms_err = stats%average
35 
36    do n=1, iv%info(sound)%nlocal
37       do k=1, iv%info(sound)%levels(n)
38          if (iv%info(sound)%proc_domain(1,n)) then
39             call da_stats_calculate (n, k, iv%sound(n)%u(k)%qc, & 
40                re%sound(n)%u(k), nu, &
41                stats%minimum%u, stats%maximum%u, &
42                stats%average%u, stats%rms_err%u)
43             call da_stats_calculate (n, k, iv%sound(n)%v(k)%qc, & 
44                re%sound(n)%v(k), nv, &
45                stats%minimum%v, stats%maximum%v, &
46                stats%average%v, stats%rms_err%v)
47             call da_stats_calculate (n, k, iv%sound(n)%t(k)%qc, & 
48                re%sound(n)%t(k), nt, &
49                stats%minimum%t, stats%maximum%t, &
50                stats%average%t, stats%rms_err%t)
51             call da_stats_calculate (n, k, iv%sound(n)%q(k)%qc, & 
52                re%sound(n)%q(k), nq, &
53                stats%minimum%q, stats%maximum%q, &
54                stats%average%q, stats%rms_err%q)
55          end if
56       end do
57    end do
58 
59    ! Do inter-processor communication to gather statistics.
60    call da_proc_sum_int (nu)
61    call da_proc_sum_int (nv)
62    call da_proc_sum_int (nt)
63    call da_proc_sum_int (nq)
64    iv%nstats(sound) = nu + nv + nt + nq
65 
66    call da_proc_stats_combine(stats%average%u, stats%rms_err%u, &
67       stats%minimum%u%value, stats%maximum%u%value, &
68       stats%minimum%u%n, stats%maximum%u%n, &
69       stats%minimum%u%l, stats%maximum%u%l)
70    call da_proc_stats_combine(stats%average%v, stats%rms_err%v, &
71       stats%minimum%v%value, stats%maximum%v%value, &
72       stats%minimum%v%n, stats%maximum%v%n, &
73       stats%minimum%v%l, stats%maximum%v%l)
74    call da_proc_stats_combine(stats%average%t, stats%rms_err%t, &
75       stats%minimum%t%value, stats%maximum%t%value, &
76       stats%minimum%t%n, stats%maximum%t%n, &
77       stats%minimum%t%l, stats%maximum%t%l)
78    call da_proc_stats_combine(stats%average%q, stats%rms_err%q, &
79       stats%minimum%q%value, stats%maximum%q%value, &
80       stats%minimum%q%n, stats%maximum%q%n, &
81       stats%minimum%q%l, stats%maximum%q%l)
82 
83    if (rootproc) then
84       if (nu /= 0 .or. nv /= 0 .or. nt /= 0 .or. nq /= 0) then
85          write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of AO for sound'
86          call da_print_stats_sound(stats_unit, nu, nv, nt, nq, stats)
87       end if
88    end if
89 
90    if (trace_use_dull) call da_trace_exit("da_ao_stats_sound")
91 
92 end subroutine da_ao_stats_sound
93 
94