Changeset 1774 for LMDZ5/trunk/libf/phylmd/cv3_routines.F

Timestamp:

Jun 18, 2013, 2:02:30 PM (11 years ago)

Author:

Laurent Fairhead

Message:

Initialisation de variables dans la convection

JYG

---

Missing initialisations in convection

JYG

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F

@@ -3561 +3561 @@
       SUBROUTINE cv3_tracer(nloc,len,ncum,nd,na,
      &                       ment,sigij,da,phi,phi2,d1a,dam,
-     &                       ep,Vprecip,elij,clw,icb,inb)
+     &                       ep,Vprecip,elij,clw,epmlmMm,eplaMm,
+     &                       icb,inb)
         implicit none
 
@@ -3577 +3578 @@
         real phi2(nloc,na,na)
         real d1a(nloc,na),dam(nloc,na)
+        real epmlmMm(nloc,na,na),eplaMm(nloc,na)
 ! variables pour tracer dans precip de l'AA et des mel
 c local variables:
@@ -3595 +3597 @@
        dam(:,:)=0.
        epm(:,:,:)=0.
-c
-!  fraction deau condensee dans les melanges convertie en precip
+       eplaMm(:,:)=0.
+       epmlmMm(:,:,:)=0.
+       phi(:,:,:)=0.
+       phi2(:,:,:)=0.
+c
+!  fraction deau condensee dans les melanges convertie en precip : epm
+! et eau condensée précipitée dans masse d'air saturé : l_m*dM_m/dzdz.dzdz
         do j=1,na
          do k=1,na
            do i=1,ncum
             if(k.ge.icb(i).and.k.le.inb(i).and.
-     &         j.ge.k.and.j.le.inb(i)) then
+!!jyg     &         j.ge.k.and.j.le.inb(i)) then
 !!jyg             epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/elij(i,k,j)
+     &         j.gt.k.and.j.le.inb(i)) then
              epm(i,j,k)=1.-(1.-ep(i,j))*clw(i,j)/
      &                     max(elij(i,k,j),1.e-16)
-
+!!
              epm(i,j,k)=max(epm(i,j,k),0.0)
             endif
@@ -3611 +3619 @@
          end do 
         end do 
+
+!
+        do j=1,na
+         do k=1,na
+           do i=1,ncum
+            if(k.ge.icb(i).and.k.le.inb(i)) then
+             eplaMm(i,j)=eplaMm(i,j) + ep(i,j)*clw(i,j)
+     &                  *ment(i,j,k)*(1.-sigij(i,j,k))
+            endif
+           end do 
+         end do 
+        end do
+!
+        do j=1,na
+         do k=1,j-1
+           do i=1,ncum
+            if(k.ge.icb(i).and.k.le.inb(i).and.
+     &         j.le.inb(i)) then
+             epmlmMm(i,j,k)=epm(i,j,k)*elij(i,k,j)*ment(i,k,j)
+            endif
+           end do 
+         end do 
+        end do
 
 !  matrices pour calculer la tendance des concentrations dans cvltr.F90
@@ -3625 +3656 @@
 
              phi2(i,j,k)=phi(i,j,k)*epm(i,j,k)   
-            else
-             dam(i,j)=0.
-             phi2(i,j,k)=0.
             endif
             end do