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

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

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

trunk/Sources/phylmd/CV3_routines/cv3_mixing.f revision 145 by guez, Tue Jun 16 15:23:29 2015 UTC trunk/Sources/phylmd/CV30_routines/cv30_mixing.f revision 190 by guez, Thu Apr 14 15:15:56 2016 UTC
# Line 1  Line 1 
1  module cv3_mixing_m  module cv30_mixing_m
2    
3    implicit none    implicit none
4    
5  contains  contains
6    
7    SUBROUTINE cv3_mixing(nloc, ncum, nd, na, icb, nk, inb, t, rr, rs, u, v, h, &    SUBROUTINE cv30_mixing(nloc, ncum, nd, na, icb, nk, inb, t, rr, rs, u, v, h, &
8         lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, sij, elij, &         lv, hp, ep, clw, m, sig, ment, qent, uent, vent, nent, sij, elij, &
9         ments, qents)         ments, qents)
10      use cv3_param_m      use cv30_param_m
11      use cvthermo      use cv_thermo_m
12    
13      !---------------------------------------------------------------------      !
14      ! a faire:      ! a faire:
15      ! - changer rr(il, 1) -> qnk(il)      ! - changer rr(il, 1) -> qnk(il)
16      ! - vectorisation de la partie normalisation des flux (do 789)      ! - vectorisation de la partie normalisation des flux (do 789)
17      !---------------------------------------------------------------------      !
18    
19      ! inputs:      ! inputs:
20      integer, intent(in):: ncum, nd, na, nloc      integer, intent(in):: ncum, nd, na, nloc
21      integer icb(nloc), inb(nloc), nk(nloc)      integer, intent(in):: icb(nloc), inb(nloc), nk(nloc)
22      real sig(nloc, nd)      real sig(nloc, nd)
23      real t(nloc, nd), rr(nloc, nd), rs(nloc, nd)      real t(nloc, nd), rr(nloc, nd), rs(nloc, nd)
24      real u(nloc, nd), v(nloc, nd)      real u(nloc, nd), v(nloc, nd)
# Line 31  contains Line 31  contains
31      real uent(nloc, na, na), vent(nloc, na, na)      real uent(nloc, na, na), vent(nloc, na, na)
32      real sij(nloc, na, na), elij(nloc, na, na)      real sij(nloc, na, na), elij(nloc, na, na)
33      real ments(nloc, nd, nd), qents(nloc, nd, nd)      real ments(nloc, nd, nd), qents(nloc, nd, nd)
     real sigij(nloc, nd, nd)  
34      integer nent(nloc, nd)      integer nent(nloc, nd)
35    
36      ! local variables:      ! local variables:
# Line 45  contains Line 44  contains
44      real zm(nloc, na)      real zm(nloc, na)
45      logical lwork(nloc)      logical lwork(nloc)
46    
47      ! --- INITIALIZE VARIOUS ARRAYS USED IN THE COMPUTATIONS      ! INITIALIZE VARIOUS ARRAYS USED IN THE COMPUTATIONS
48    
49      do j=1, nl      do j=1, nl
50         do i=1, ncum         do i=1, ncum
51            nent(i, j)=0            nent(i, j)=0
52            ! in convect3, m is computed in cv3_closure            ! in convect3, m is computed in cv30_closure
53            ! ori m(i, 1)=0.0            ! ori m(i, 1)=0.0
54         end do         end do
55      end do      end do
# Line 74  contains Line 73  contains
73    
74      zm(:, :)=0.      zm(:, :)=0.
75    
76      ! --- CALCULATE ENTRAINED AIR MASS FLUX (ment), TOTAL WATER MIXING      ! CALCULATE ENTRAINED AIR MASS FLUX (ment), TOTAL WATER MIXING
77      ! --- RATIO (QENT), TOTAL CONDENSED WATER (elij), AND MIXING      ! RATIO (QENT), TOTAL CONDENSED WATER (elij), AND MIXING
78      ! --- FRACTION (sij)      ! FRACTION (sij)
79    
80      do i=minorig+1, nl      do i=minorig+1, nl
81    
# Line 121  contains Line 120  contains
120            end do            end do
121         end do         end do
122    
123         ! *** if no air can entrain at level i assume that updraft detrains ***         !  if no air can entrain at level i assume that updraft detrains
124         ! *** at that level and calculate detrained air flux and properties ***         !  at that level and calculate detrained air flux and properties
   
        !@ do 170 i=icb(il), inb(il)  
125    
126         do il=1, ncum         do il=1, ncum
127            if ((i >= icb(il)).and.(i <= inb(il)).and.(nent(il, i) == 0)) then            if ((i >= icb(il)).and.(i <= inb(il)).and.(nent(il, i) == 0)) then
# Line 140  contains Line 137  contains
137         end do         end do
138      end do      end do
139    
140      do j=minorig, nl      ! NORMALIZE ENTRAINED AIR MASS FLUXES
141         do i=minorig, nl      ! TO REPRESENT EQUAL PROBABILITIES OF MIXING
           do il=1, ncum  
              if ((j >= (icb(il)-1)).and.(j <= inb(il)) &  
                   .and.(i >= icb(il)).and.(i <= inb(il)))then  
                 sigij(il, i, j)=sij(il, i, j)  
              endif  
           end do  
        end do  
     end do  
     !@ enddo  
   
     !@170 continue  
   
     ! --- NORMALIZE ENTRAINED AIR MASS FLUXES  
     ! --- TO REPRESENT EQUAL PROBABILITIES OF MIXING  
142    
143      call zilch(asum, nloc*nd)      asum = 0.
144      call zilch(csum, nloc*nd)      csum = 0.
     call zilch(csum, nloc*nd)  
145    
146      do il=1, ncum      do il=1, ncum
147         lwork(il) = .FALSE.         lwork(il) = .FALSE.
# Line 331  contains Line 313  contains
313         enddo         enddo
314      enddo      enddo
315    
316    end SUBROUTINE cv3_mixing    end SUBROUTINE cv30_mixing
317    
318  end module cv3_mixing_m  end module cv30_mixing_m
Legend:
Removed from v.145  
changed lines
  Added in v.190
  ViewVC Help
Powered by ViewVC 1.1.21