1 |
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 |