1 |
SUBROUTINE cv3_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi) |
module cv3_tracer_m |
2 |
|
|
3 |
implicit none |
implicit none |
4 |
|
|
5 |
! Inputs: |
contains |
|
integer ncum, nd, na, nloc, len |
|
|
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 |
|
6 |
|
|
7 |
end SUBROUTINE cv3_tracer |
SUBROUTINE cv3_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi) |
8 |
|
|
9 |
|
implicit none |
10 |
|
|
11 |
|
! Inputs: |
12 |
|
integer, intent(in):: ncum, nd, na, nloc, len |
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 |