phylmd/CV30_routines/cv3_tracer.f

module cv3_tracer_m

  implicit none

contains

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

    implicit none

    ! Inputs:
    integer, intent(in):: ncum, na, nloc
    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

end module cv3_tracer_m
1	guez	97	module cv3_tracer_m
2	guez	47
3			implicit none
4
5	guez	97	contains
6	guez	47
7	guez	139	SUBROUTINE cv3_tracer(nloc, ncum, na, ment, sij, da, phi)
8	guez	47
9	guez	97	implicit none
10	guez	47
11	guez	97	! Inputs:
12	guez	139	integer, intent(in):: ncum, na, nloc
13	guez	97	real ment(nloc, na, na), sij(nloc, na, na)
14	guez	47
15	guez	97	! Ouputs:
16			real da(nloc, na), phi(nloc, na, na)
17	guez	47
18	guez	97	! Local variables:
19			integer i, j, k
20	guez	47
21	guez	97	!------------------------------------------------------------
22
23			da = 0.
24
25			do j = 1, na
26			do k = 1, na
27			do i = 1, ncum
28			da(i, j) = da(i, j)+(1.-sij(i, k, j))*ment(i, k, j)
29			phi(i, j, k) = sij(i, k, j) * ment(i, k, j)
30			end do
31			end do
32			end do
33
34			end SUBROUTINE cv3_tracer
35
36			end module cv3_tracer_m