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 |