da_couple.inc

1 subroutine da_couple (mut, field, msf, &
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, jds:jde),          intent(in  ) :: mut, msf
9   
10    real, dimension(ids:ide, jds:jde, kds:kde), intent(inout) :: field
11   
12   
13    integer :: i, j, k
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 end subroutine da_couple