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    if (trace_use) call da_trace_entry("da_print_stats_radar")
14    
15    write(unit=stats_unit, fmt='(a/)') &
16       '   var             rv (m/s)        n    k         rf (dBZ)        n    k  '
17 
18    write(unit=stats_unit, fmt='(a,i16,i31)') &
19       '  Number: ', nrv, nrf
20 
21    if (nrv < 1) nrv = 1
22    if (nrf < 1) nrf = 1
23    
24    write(unit=stats_unit, fmt='((a,f12.4,i9,i5, f17.4,i9,i5))')  &
25       ' Minimum(n,k): ', radar%minimum%rv, radar%minimum%rf, &
26       ' Maximum(n,k): ', radar%maximum%rv, radar%maximum%rf
27    write(unit=stats_unit, fmt='((a,2(f12.4,19x)))')              &
28       ' Average     : ', radar%average%rv/real(nrv), radar%average%rf/real(nrf), &
29       '    RMSE     : ', sqrt(radar%rms_err%rv/real(nrv)), &
30       sqrt(radar%rms_err%rf/real(nrf))
31 
32    if (trace_use) call da_trace_exit("da_print_stats_radar")
33 
34 end subroutine da_print_stats_radar
35 
36