/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_tracer.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_tracer.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/CV3_routines/cv3_tracer.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/Sources/phylmd/CV30_routines/cv30_tracer.f revision 189 by guez, Tue Mar 29 15:20:23 2016 UTC
# Line 1  Line 1 
1  SUBROUTINE cv3_tracer(nloc, len, ncum, nd, na, ment, sij, da, phi)  module cv30_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 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

Legend:
Removed from v.82  
changed lines
  Added in v.189

  ViewVC Help
Powered by ViewVC 1.1.21