var/da/da_tools/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