var/da/da_chem_sfc/da_ao_stats_chem_sfc.inc

   1 subroutine da_ao_stats_chem_sfc (stats_unit, iv, re, ob)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    integer,        intent (in)    :: stats_unit    ! Output unit for stats.
  10    type (iv_type), intent (inout) :: iv            ! iv
  11    type  (y_type), intent (in)    :: re            ! A - O
  12    type(y_type),   intent (in)    :: ob            ! Observation structure.
  13
  14    type (stats_chem_sfc_type) :: stats
  15    integer                 :: nchem
  16    integer                 :: n, ichem
  17    real                    :: chem_inc, v_inc
  18
  19    if (trace_use_dull) call da_trace_entry("da_ao_stats_chem_sfc")
  20
  21
  22    allocate(stats%maximum%chem(num_chemic_surf))
  23    allocate(stats%minimum%chem(num_chemic_surf))
  24    allocate(stats%average%chem(num_chemic_surf))
  25    allocate(stats%rms_err%chem(num_chemic_surf))
  26
  27    do ichem=PARAM_FIRST_SCALAR, num_chemic_surf
  28     stats%maximum%chem(ichem) = maxmin_type (0.0, 0, 0)
  29     stats%minimum%chem(ichem) = maxmin_type (0.0, 0, 0)
  30     stats%average%chem(ichem) = 0.0
  31     stats%rms_err%chem(ichem) = 0.0
  32    end do
  33
  34    iv%nstats(chemic_surf) = 0
  35
  36    do ichem=PARAM_FIRST_SCALAR, num_chemic_surf
  37    nchem=0
  38
  39    do n=1, iv%info(chemic_surf)%nlocal
  40       if (iv%info(chemic_surf)%proc_domain(1,n)) then
  41
  42          chem_inc = re%chemic_surf(n)%chem(ichem)
  43
  44          call da_stats_calculate (n, 0, iv%chemic_surf(n)%chem(ichem)%qc,  &
  45             chem_inc, nchem, &
  46             stats%minimum%chem(ichem), stats%maximum%chem(ichem), &
  47             stats%average%chem(ichem), stats%rms_err%chem(ichem))
  48       end if
  49    end do
  50
  51    ! Do inter-processor communication to gather statistics.
  52    call da_proc_sum_int (nchem)
  53
  54    iv%nstats(chemic_surf) = iv%nstats(chemic_surf)+nchem
  55
  56    call da_proc_stats_combine(stats%average%chem(ichem), stats%rms_err%chem(ichem), &
  57       stats%minimum%chem(ichem)%value, stats%maximum%chem(ichem)%value, &
  58       stats%minimum%chem(ichem)%n, stats%maximum%chem(ichem)%n, &
  59       stats%minimum%chem(ichem)%l, stats%maximum%chem(ichem)%l)
  60
  61    if (rootproc) then
  62       if (nchem /= 0 ) then
  63          write(unit=stats_unit, fmt='(/a/)') ' Diagnostics of AO for chem'
  64          if (chemicda_opt == 1 .or. chemicda_opt == 2) then
  65            call da_print_stats_chem_sfc(stats_unit, nchem, stats)
  66          else if (chemicda_opt >= 3) then
  67               call da_print_stats_chem_sfc(stats_unit, nchem, stats, ichem)
  68          end if
  69       end if
  70    end if
  71    end do
  72
  73     deallocate(stats%maximum%chem)
  74     deallocate(stats%minimum%chem)
  75     deallocate(stats%average%chem)
  76     deallocate(stats%rms_err%chem)
  77
  78    if (trace_use_dull) call da_trace_exit("da_ao_stats_chem_sfc")
  79
  80 end subroutine da_ao_stats_chem_sfc
  81
  82