var/da/da_tools/da_get_3d_sum.inc

   1 subroutine da_get_3d_sum(var, name)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real,             intent(in) :: var(ims:ime, jms:jme, kms:kme)
  10    character(len=*), intent(in) :: name
  11
  12    real :: partial, total
  13
  14    if (trace_use) call da_trace_entry("da_get_3d_sum")
  15
  16    partial = sum(var(its:ite,jts:jte,kts:kte)*var(its:ite,jts:jte,kts:kte))
  17
  18    total = wrf_dm_sum_real (partial)
  19
  20    write(unit=stdout, fmt='(3a, e24.14)') 'Square sum of <', trim(name), '>=', total
  21
  22    if (trace_use) call da_trace_exit("da_get_3d_sum")
  23
  24 end subroutine da_get_3d_sum
  25
  26