da_set_boundary_3d.inc

1 subroutine da_set_boundary_3d(var)
2    !------------------------------------------------------------------------
3    !  Purpose: 
4    !
5    !  Merge East-West boundary values for input 3d-array (var)
6    !------------------------------------------------------------------------
7 
8    implicit none
9 
10    real, intent(inout) :: var(ims:ime, jms:jme, kms:kme)
11 
12    integer :: n, j, k
13 
14    if ((its /= ids) .or. (ite /= ide)) return
15 
16    if (trace_use) call da_trace_entry("da_set_boundary_3d")
17 
18    do k=kts, kte
19       do j=jts, jte
20          do n=1,bdyzone
21             var(ids-n,j,k) = var(ide+1-n,j,k)
22             var(ide+n,j,k) = var(ids-1+n,j,k)
23          end do
24       end do
25    end do
26 
27    if (trace_use) call da_trace_exit("da_set_boundary_3d")
28 
29 end subroutine da_set_boundary_3d
30 
31