da_print_stats_airep.inc

1 subroutine da_print_stats_airep(stats_unit, nu, nv, nt, Airep)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7     implicit none
8 
9     integer,           intent(in)    :: stats_unit
10     integer,           intent(inout) :: nu, nv, nt
11     type (stats_airep_type), intent(in):: Airep
12 
13     write(unit=stats_unit, fmt='(4a/)') &
14        '   var             ', &
15        'u (m/s)     n    k    ', &
16        'v (m/s)     n    k    ', &
17        't (K)       n    k    '
18 
19     write(unit=stats_unit, fmt='(a,i16,3i22)') &
20        '  Number: ', nu, nv, nt
21 
22     if (nu < 1) nu = 1
23     if (nv < 1) nv = 1
24     if (nt < 1) nt = 1
25 
26     write(unit=stats_unit, fmt='((a,3(f12.4,2i5)))') &
27        ' Minimum(n,k): ', Airep%minimum%u, Airep%minimum%v, Airep%minimum%t, &
28        ' Maximum(n,k): ', Airep%maximum%u, Airep%maximum%v, Airep%maximum%t
29     write(unit=stats_unit, fmt='((a,3(f12.4,10x)))') &
30        ' Average     : ', Airep%average%u/real(nu), &
31                           Airep%average%v/real(nv), &
32                           Airep%average%t/real(nt), &
33        '    RMSE     : ', sqrt(Airep%rms_err%u/real(nu)), &
34                          sqrt(Airep%rms_err%v/real(nv)), &
35                          sqrt(Airep%rms_err%t/real(nt))
36 
37 end subroutine da_print_stats_airep
38 
39