da_print_stats_ssmt2.inc

1 subroutine da_print_stats_ssmt2(stats_unit, nrh, minimum, maximum, &
2                                  average, rms_err)
3 
4    !-----------------------------------------------------------------------
5    ! Purpose: TBD
6    !-----------------------------------------------------------------------
7 
8    implicit none
9 
10    integer,            intent(in)    :: stats_unit
11    integer,            intent(inout) :: nrh
12    type (maxmin_type), intent(in)    :: minimum
13    type (maxmin_type), intent(in)    :: maximum
14    real,               intent(in)    :: average
15    real,               intent(in)    :: rms_err
16 
17    if (trace_use_dull) call da_trace_entry("da_print_stats_ssmt2")
18    
19    write(unit=stats_unit, fmt='(a/)') '   var             rh(K)     n    k'
20 
21    write(unit=stats_unit, fmt='(a,i16,4i22)') &
22         '  Number: ', nrh
23 
24    if (nrh < 1) nrh = 1
25    
26    write(unit=stats_unit, fmt='((a,f12.4,2i5))') &
27         ' Minimum(n,k): ', minimum,    &
28         ' Maximum(n,k): ', maximum
29    write(unit=stats_unit, fmt='((a,f12.4,10x))') &
30         ' Average     : ', average/real(nrh),    &
31         '    RMSE     : ', sqrt(rms_err/real(nrh))
32 
33 end subroutine da_print_stats_ssmt2
34 
35