<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE> <A NAME='DA_COUPLE_UV'><A href='../../html_code/update_bc/da_couple_uv.inc.html#DA_COUPLE_UV' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>subroutine da_couple_uv (u, v, mu, mub, msfu, msfv, ids, ide, jds, jde, kds, kde) 3,3 implicit none integer, intent(in) :: ids, ide, jds, jde, kds, kde real, intent(inout) :: u(ids:ide+1,jds:jde,kds:kde) real, intent(inout) :: v(ids:ide,jds:jde+1,kds:kde) real, intent(in) :: msfu(ids:ide+1,jds:jde) real, intent(in) :: msfv(ids:ide,jds:jde+1) real, intent(in) :: mu(ids:ide,jds:jde) real, intent(in) :: mub(ids:ide,jds:jde) real, allocatable :: muu(:,:), muv(:,:) ! update_bc ! if (trace_use) call da_trace_entry("da_couple_uv") allocate(muu(ids:ide+1, jds:jde )) allocate(muv(ids:ide , jds:jde+1)) ! couple variables u, v call da_calc_mu_uv
(mu, mub, muu, muv, ids, ide, jds, jde) call da_couple
(muu, u, msfu, ids, ide+1, jds, jde, kds, kde) call da_couple
(muv, v, msfv, ids, ide, jds, jde+1, kds, kde) deallocate(muu) deallocate(muv) ! if (trace_use) call da_trace_exit("da_couple_uv") end subroutine da_couple_uv