var/da/da_update_bc/da_calc_mu_uv.inc

   1 subroutine da_calc_mu_uv (mu, mub, muu, muv, &
   2                         ids, ide, jds, jde)
   3
   4    implicit none
   5
   6    integer, intent(in) :: ids, ide, jds, jde
   7
   8    real, dimension(ids:ide, jds:jde),     intent(in  ) :: mu,  mub
   9
  10    real, dimension(ids:ide+1, jds:jde ), intent( out) :: muu
  11    real, dimension(ids:ide  , jds:jde+1), intent( out) :: muv
  12
  13    real, dimension(ids-1:ide+1, jds-1:jde+1) :: mut
  14
  15    integer :: i, j
  16
  17    do j=jds,jde
  18       do i=ids,ide
  19          mut(i,j) = mu(i,j)+mub(i,j)
  20       end do
  21
  22       mut(ids-1,j) = mut(ids,j)
  23       mut(ide+1,j) = mut(ide,j)
  24    end do
  25
  26    do i=ids-1,ide+1
  27       mut(i,jds-1)=mut(i,jds)
  28       mut(i,jde+1)=mut(i,jde)
  29    end do
  30
  31    do j=jds,jde
  32       do i=ids,ide+1
  33          muu(i,j) = 0.5*(mut(i,j)+mut(i-1,j))
  34       end do
  35    end do
  36
  37    do j=jds,jde+1
  38       do i=ids,ide
  39          muv(i,j) = 0.5*(mut(i,j)+mut(i,j-1))
  40       end do
  41    end do
  42
  43 end subroutine da_calc_mu_uv
  44
  45