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    write(unit=stats_unit, fmt='(a/)') '   var             rh(K)     n    k'
18 
19    write(unit=stats_unit, fmt='(a,i16,4i22)') &
20         '  Number: ', nrh
21 
22    if (nrh < 1) nrh = 1
23    
24    write(unit=stats_unit, fmt='((a,f12.4,2i5))') &
25         ' Minimum(n,k): ', minimum,    &
26         ' Maximum(n,k): ', maximum
27    write(unit=stats_unit, fmt='((a,f12.4,10x))') &
28         ' Average     : ', average/real(nrh),    &
29         '    RMSE     : ', sqrt(rms_err/real(nrh))
30 
31 end subroutine da_print_stats_ssmt2
32 
33