subroutine da_print_stats_rain(stats_unit, nrain, rain)

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

   implicit none

   integer,                intent(in)    :: stats_unit
   integer,                intent(inout) :: nrain 
   type (stats_rain_type), intent(in)  :: rain

   if (trace_use) call da_trace_entry("da_print_stats_rain")

   if (nrain > 0) then

      write(unit=stats_unit, fmt='(a/)') '   var        rainfall(mm)     n'

      write(unit=stats_unit, fmt='(a,i14)') '  Number: ', nrain

      write(unit=stats_unit, fmt='(a, f12.4,i5)') &
         ' Minimum(n): ', rain%minimum%rain%value, &
                          rain%minimum%rain%n    , &
         ' Maximum(n): ', rain%maximum%rain%value, &
                          rain%maximum%rain%n
      write(unit=stats_unit, fmt='(a, f12.4,5x)') &
         ' Average   : ', rain%average%rain/real(nrain), &
         '    RMSE   : ', sqrt(rain%rms_err%rain/real(nrain))
   end if

   if (trace_use) call da_trace_exit("da_print_stats_rain")

end subroutine da_print_stats_rain