subroutine da_print_stats_radar(stats_unit, nrv, nrf, radar)

   !-----------------------------------------------------------------------
   ! Purpose: TBD
   !-----------------------------------------------------------------------

   implicit none

   integer,                 intent(in)    :: stats_unit
   integer,                 intent(inout) :: nrv, nrf
   type (stats_radar_type), intent(in)    :: radar

   if (trace_use) call da_trace_entry("da_print_stats_radar")
   
   write(unit=stats_unit, fmt='(a/)') &
      '   var             rv (m/s)        n    k         rf (dBZ)        n    k  '

   write(unit=stats_unit, fmt='(a,i16,i31)') &
      '  Number: ', nrv, nrf

   if (nrv < 1) nrv = 1
   if (nrf < 1) nrf = 1
   
   write(unit=stats_unit, fmt='((a,f12.4,i9,i5, f17.4,i9,i5))')  &
      ' Minimum(n,k): ', radar%minimum%rv, radar%minimum%rf, &
      ' Maximum(n,k): ', radar%maximum%rv, radar%maximum%rf
   write(unit=stats_unit, fmt='((a,2(f12.4,19x)))')              &
      ' Average     : ', radar%average%rv/real(nrv), radar%average%rf/real(nrf), &
      '    RMSE     : ', sqrt(radar%rms_err%rv/real(nrv)), &
      sqrt(radar%rms_err%rf/real(nrf))

   if (trace_use) call da_trace_exit("da_print_stats_radar")

end subroutine da_print_stats_radar