da_couple.inc

1 subroutine da_couple (mut, field, msf, ids, ide, jds, jde, kds, kde)
2 
3    implicit none
4 
5    integer, intent(in)    :: ids, ide, jds, jde, kds, kde
6    real,    intent(in)    :: mut(ids:ide, jds:jde)
7    real,    intent(in)    :: msf(ids:ide, jds:jde)
8    real,    intent(inout) :: field(ids:ide, jds:jde, kds:kde)
9   
10    integer :: i, j, k
11 
12    ! update_bc
13    ! if (trace_use) call da_trace_entry("da_couple")
14   
15    do j=jds,jde
16       do k=kds,kde
17          do i=ids,ide
18             field(i,j,k)=field(i,j,k)*mut(i,j)/msf(i,j)
19          end do
20       end do
21    end do
22 
23    ! if (trace_use) call da_trace_exit("da_couple")
24 
25 end subroutine da_couple