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	module cv3_tracer_m
2
3	implicit none
4
5	contains
6
7	SUBROUTINE cv3_tracer(nloc, ncum, na, ment, sij, da, phi)
8
9	implicit none
10
11	! Inputs:
12	integer, intent(in):: ncum, na, nloc
13	real ment(nloc, na, na), sij(nloc, na, na)
14
15	! Ouputs:
16	real da(nloc, na), phi(nloc, na, na)
17
18	! Local variables:
19	integer i, j, k
20
21	!------------------------------------------------------------
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