var/da/da_lightning/da_print_stats_lightning.inc

   1 subroutine da_print_stats_lightning(stats_unit, nw, ndiv, nqv, lightning)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    ! Authors: Z Chen (zchen@fjnu.edu.cn), Jenny Sun (NCAR), X Qie (IAP)
   6    !-----------------------------------------------------------------------
   7
   8    implicit none
   9
  10    integer,                 intent(in)    :: stats_unit
  11    integer,                 intent(inout) :: nw, ndiv, nqv
  12    type (stats_lightning_type), intent(in):: lightning
  13
  14    if (trace_use_dull) call da_trace_entry("da_print_stats_lightning")
  15
  16    write(unit=stats_unit, fmt='(a/)') &
  17       '   var              w (m/s)        n    k        div (1/s)        n    k        qv (kg/kg)       n    k'
  18
  19    write(unit=stats_unit, fmt='(a,(i16,2i31))') &
  20       '  Number: ', nw, ndiv, nqv
  21
  22    if (nw < 1)   nw = 1
  23    if (ndiv < 1) ndiv = 1
  24    if (nqv < 1)  nqv = 1
  25
  26    write(unit=stats_unit, fmt='((a,f12.4,i9,i5, 2(f17.4,i9,i5)))')  &
  27         ' Minimum(n,k): ', lightning%minimum%w, lightning%minimum%div, lightning%minimum%qv
  28    write(unit=stats_unit, fmt='((a,f12.4,i9,i5, 2(f17.4,i9,i5)))')  &
  29         ' Maximum(n,k): ', lightning%maximum%w, lightning%maximum%div, lightning%maximum%qv
  30    write(unit=stats_unit, fmt='((a,3(f12.4,19x)))') &
  31         ' Average     : ', lightning%average%w/real(nw), &
  32                            lightning%average%div/real(ndiv), &
  33                            lightning%average%qv/real(nqv)
  34    write(unit=stats_unit, fmt='((a,3(f12.4,19x)))') &
  35         '    RMSE     : ', sqrt(lightning%rms_err%w/real(nw)), &
  36                            sqrt(lightning%rms_err%div/real(ndiv)), &
  37                            sqrt(lightning%rms_err%qv/real(nqv))
  38
  39    if (trace_use_dull) call da_trace_exit("da_print_stats_lightning")
  40
  41 end subroutine da_print_stats_lightning