da_print_stats_gpspw.inc

1 subroutine da_print_stats_gpspw(stats_unit, ntpw, Gpspw)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer,           intent(in)    :: stats_unit
10    integer,           intent(inout) :: ntpw
11    type (stats_gpspw_type), intent(in)    :: Gpspw
12 
13    if (ntpw > 0) then
14       write(unit=stats_unit, fmt='(a/)') '   var           tpw(cm)     n'
15  
16       write(unit=stats_unit, fmt='(a,i14)') '  Number: ', ntpw
17 
18       write(unit=stats_unit, fmt='(a, f12.4,i5)') &
19         ' Minimum(n): ', gpspw%minimum%tpw%value, &
20                          gpspw%minimum%tpw%n    , &
21         ' Maximum(n): ', gpspw%maximum%tpw%value, &
22                          gpspw%maximum%tpw%n
23       write(unit=stats_unit, fmt='(a, f12.4,5x)') &
24         ' Average   : ', gpspw%average%tpw/real(ntpw), &
25         '    RMSE   : ', sqrt(gpspw%rms_err%tpw/real(ntpw))
26    end if
27 
28 end subroutine da_print_stats_gpspw
29 
30