var/da/da_update_bc/da_couple.inc

   1 subroutine da_couple (mut, field, msf, c1, c2, 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(in)    :: c1(kds:kde)
   9    real,    intent(in)    :: c2(kds:kde)
  10    real,    intent(inout) :: field(ids:ide, jds:jde, kds:kde)
  11
  12    integer :: i, j, k
  13
  14    ! update_bc
  15    ! if (trace_use) call da_trace_entry("da_couple")
  16
  17    do k=kds,kde
  18       do j=jds,jde
  19          do i=ids,ide
  20             field(i,j,k)=field(i,j,k)*(c1(k)*mut(i,j)+c2(k))/msf(i,j)
  21          end do
  22       end do
  23    end do
  24
  25    ! if (trace_use) call da_trace_exit("da_couple")
  26
  27 end subroutine da_couple