da_print_stats_radar.inc

1 subroutine da_print_stats_radar(stats_unit, nrv, nrf, Radar)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer,                 intent(in)    :: stats_unit
10    integer,                 intent(inout) :: nrv, nrf
11    type (stats_Radar_type), intent(in)    :: Radar
12    
13    write(unit=stats_unit, fmt='(a/)') &
14       '   var             rv (m/s)        n    k         rf (dBZ)        n    k  '
15 
16    write(unit=stats_unit, fmt='(a,i16,i31)') &
17       '  Number: ', nrv, nrf
18 
19    if (nrv < 1) nrv = 1
20    if (nrf < 1) nrf = 1
21    
22    write(unit=stats_unit, fmt='((a,f12.4,i9,i5, f17.4,i9,i5))')  &
23       ' Minimum(n,k): ', Radar%minimum%rv, Radar%minimum%rf, &
24       ' Maximum(n,k): ', Radar%maximum%rv, Radar%maximum%rf
25    write(unit=stats_unit, fmt='((a,2(f12.4,19x)))')              &
26       ' Average     : ', Radar%average%rv/real(nrv), Radar%average%rf/real(nrf), &
27       '    RMSE     : ', sqrt(Radar%rms_err%rv/real(nrv)), &
28       sqrt(Radar%rms_err%rf/real(nrf))
29 
30 end subroutine da_print_stats_radar
31 
32