var/da/da_ssmi/da_print_stats_ssmt2.inc

   1 subroutine da_print_stats_ssmt2(stats_unit, nrh, minimum, maximum, &
   2                                  average, rms_err)
   3
   4    !-----------------------------------------------------------------------
   5    ! Purpose: TBD
   6    !-----------------------------------------------------------------------
   7
   8    implicit none
   9
  10    integer,            intent(in)    :: stats_unit
  11    integer,            intent(inout) :: nrh
  12    type (maxmin_type), intent(in)    :: minimum
  13    type (maxmin_type), intent(in)    :: maximum
  14    real,               intent(in)    :: average
  15    real,               intent(in)    :: rms_err
  16
  17    if (trace_use_dull) call da_trace_entry("da_print_stats_ssmt2")
  18
  19    write(unit=stats_unit, fmt='(a/)') '   var             rh(K)     n    k'
  20
  21    write(unit=stats_unit, fmt='(a,i16,4i22)') &
  22         '  Number: ', nrh
  23
  24    if (nrh < 1) nrh = 1
  25
  26    write(unit=stats_unit, fmt='((a,f12.4,2i5))') &
  27         ' Minimum(n,k): ', minimum,    &
  28         ' Maximum(n,k): ', maximum
  29    write(unit=stats_unit, fmt='((a,f12.4,10x))') &
  30         ' Average     : ', average/real(nrh),    &
  31         '    RMSE     : ', sqrt(rms_err/real(nrh))
  32
  33 end subroutine da_print_stats_ssmt2
  34
  35