da_couple_uv.inc

1 subroutine da_couple_uv (u,  v,  mu, mub,  msfu, msfv, &
2                        ids, ide, jds, jde, kds, kde)
3 
4    implicit none
5 
6    integer, intent(in) :: ids, ide, jds, jde, kds, kde
7 
8    real, dimension(ids:ide+1, jds:jde  , kds:kde), intent(inout) :: u
9    real, dimension(ids:ide  , jds:jde+1, kds:kde), intent(inout) :: v
10 
11    real, dimension(ids:ide+1, jds:jde ),          intent(in  ) :: msfu
12    real, dimension(ids:ide  , jds:jde+1),          intent(in  ) :: msfv
13 
14    real, dimension(ids:ide  , jds:jde ),          intent(in  ) :: mu, mub
15 
16    real, allocatable, dimension(:, :) :: muu, muv
17 
18    allocate(muu(ids:ide+1, jds:jde ))
19    allocate(muv(ids:ide  , jds:jde+1))
20 
21    ! couple variables u, v
22 
23    call da_calc_mu_uv (mu, mub, muu, muv, &
24                      ids, ide, jds, jde)
25 
26    call da_couple (muu, u, msfu, &
27                  ids, ide+1, jds, jde, kds, kde)
28 
29    call da_couple (muv, v, msfv, &
30                  ids, ide, jds, jde+1, kds, kde)
31 
32    deallocate(muu)
33    deallocate(muv)
34 
35 end subroutine da_couple_uv
36 
37