phylmd/CV3_routines/cv3_tracer.f90

SUBROUTINE cv3_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi)

  implicit none

  ! Inputs:
  integer ncum, nd, na, nloc, len
  real ment(nloc, na, na), sij(nloc, na, na)

  ! Ouputs:
  real da(nloc, na), phi(nloc, na, na)

  ! Local variables:
  integer i, j, k

  !------------------------------------------------------------

  da = 0.

  do j = 1, na
     do k = 1, na
        do i = 1, ncum
           da(i, j) = da(i, j)+(1.-sij(i, k, j))*ment(i, k, j)
           phi(i, j, k) = sij(i, k, j) * ment(i, k, j)
        end do
     end do
  end do

end SUBROUTINE cv3_tracer
1	SUBROUTINE cv3_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi)
2
3	implicit none
4
5	! Inputs:
6	integer ncum, nd, na, nloc, len
7	real ment(nloc, na, na), sij(nloc, na, na)
8
9	! Ouputs:
10	real da(nloc, na), phi(nloc, na, na)
11
12	! Local variables:
13	integer i, j, k
14
15	!------------------------------------------------------------
16
17	da = 0.
18
19	do j = 1, na
20	do k = 1, na
21	do i = 1, ncum
22	da(i, j) = da(i, j)+(1.-sij(i, k, j))*ment(i, k, j)
23	phi(i, j, k) = sij(i, k, j) * ment(i, k, j)
24	end do
25	end do
26	end do
27
28	end SUBROUTINE cv3_tracer