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