phylmd/CV30_routines/cv30_tracer.f90

module cv30_tracer_m

  implicit none

contains

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

    ! Passive tracers.

    integer, intent(in):: ncum, na, nloc
    real, intent(in):: ment(nloc, na, na), sij(nloc, na, na)

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

    ! Local:
    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 cv30_tracer

end module cv30_tracer_m
1	module cv30_tracer_m
2
3	implicit none
4
5	contains
6
7	SUBROUTINE cv30_tracer(nloc, ncum, na, ment, sij, da, phi)
8
9	! Passive tracers.
10
11	integer, intent(in):: ncum, na, nloc
12	real, intent(in):: ment(nloc, na, na), sij(nloc, na, na)
13
14	! Ouputs:
15	real, intent(out):: da(nloc, na)
16	real phi(nloc, na, na)
17
18	! Local:
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 cv30_tracer
35
36	end module cv30_tracer_m