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, dimension(ims:ime, jms:jme, kms:kme), intent(in) :: var
10    character(len=*), intent(in) :: name
11 
12    real :: partial, total
13 
14    partial = sum(var(its:ite,jts:jte,kts:kte)*var(its:ite,jts:jte,kts:kte))
15 
16    total = wrf_dm_sum_real( partial)
17 
18    write(unit=stdout, fmt='(3a, e24.14)') &
19         'Square sum of <', trim(name), '>=', total
20 
21 end subroutine da_get_3d_sum
22 
23