phylmd/CV30_routines/cv30_tracer.f

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	guez	47	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