var/da/da_rain/da_print_stats_rain.inc

   1 subroutine da_print_stats_rain(stats_unit, nrain, rain)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    integer,                intent(in)    :: stats_unit
  10    integer,                intent(inout) :: nrain
  11    type (stats_rain_type), intent(in)  :: rain
  12
  13    if (trace_use) call da_trace_entry("da_print_stats_rain")
  14
  15    if (nrain > 0) then
  16
  17       write(unit=stats_unit, fmt='(a/)') '   var        rainfall(mm)     n'
  18
  19       write(unit=stats_unit, fmt='(a,i14)') '  Number: ', nrain
  20
  21       write(unit=stats_unit, fmt='(a, f12.4,i5)') &
  22          ' Minimum(n): ', rain%minimum%rain%value, &
  23                           rain%minimum%rain%n    , &
  24          ' Maximum(n): ', rain%maximum%rain%value, &
  25                           rain%maximum%rain%n
  26       write(unit=stats_unit, fmt='(a, f12.4,5x)') &
  27          ' Average   : ', rain%average%rain/real(nrain), &
  28          '    RMSE   : ', sqrt(rain%rms_err%rain/real(nrain))
  29    end if
  30
  31    if (trace_use) call da_trace_exit("da_print_stats_rain")
  32
  33 end subroutine da_print_stats_rain
  34
  35