Changeset 1411 for LMDZ4/trunk/libf/phylmd/fisrtilp.F

Timestamp:

Jul 9, 2010, 1:06:15 PM (14 years ago)

Author:

idelkadi

Message:

Nettoyage dans physiq.F des anciennes versions versions des options iflag_cldcon=5 ou 6.
Dans isrtilp.F, iflag_cldcon=5 : la version bi-gaussienne partout et iflag_cldcon=6 : bi-gaussiennes pour les thermiques et la lognormale ailleurs

File:

• LMDZ4/trunk/libf/phylmd/fisrtilp.F (modified) (8 diffs)

LMDZ4/trunk/libf/phylmd/fisrtilp.F

@@ -50 +50 @@
       REAL zthl(klon,klev)
 
+      logical lognormale(klon)
+
 cAA
 c Coeffients de fraction lessivee : pour OFF-LINE
@@ -140 +142 @@
       zdelq=0.0
       
-!      print*,'CLOUDTH4 A. JAM'
+      print*,'NUAGES4 A. JAM'
       IF (appel1er) THEN
 c
@@ -289 +291 @@
      .                            /RG/dtime
 
-c pour la glace, on r�vapore toute la pr�ip dans la couche du dessous
+c pour la glace, on r�vapore toute la pr�ip dans la couche du dessous
 c la glace venant de la couche du dessus est simplement dans la couche
 c du dessous.
@@ -343 +345 @@
 c   zqn   : eau totale dans le nuage
 c   zcond : eau condensee moyenne dans la maille.
-c           on prend en compte le r�hauffement qui diminue la partie condensee
+c           on prend en compte le r�hauffement qui diminue la partie condensee
 c
 c   Version avec les raqts
@@ -365 +367 @@
            enddo 
 
-              if (iflag_cldcon.eq.5) then
+              if (iflag_cldcon>=5) then
 
                  call cloudth(klon,klev,k,ztv,
@@ -377 +379 @@
                  enddo
 
-              else
-
+              endif
+
+! Pour iflag_cldcon<=4, on prend toujours la lognormale
+! Dans le cas iflag_cldcon=5, on prend systématiquement la bi-gaussienne
+! Dans le cas iflagÃ_cldcon=6, on prend la lognormale en absence des thermiques
+
+            lognormale(:)=
+     .      iflag_cldcon<=4.or.(iflag_cldcon==6.and.fraca(:,k)<1.e-10)
             do i=1,klon
+            if (lognormale(i)) then
             zpdf_sig(i)=ratqs(i,k)*zq(i)
             zpdf_k(i)=-sqrt(log(1.+(zpdf_sig(i)/zq(i))**2))
@@ -391 +400 @@
             zpdf_e2(i)=sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i))
             zpdf_e2(i)=1.-erf(zpdf_e2(i))
+            endif
+            enddo
+
+            do i=1,klon
+            if (lognormale(i)) then
             if (zpdf_e1(i).lt.1.e-10) then
                rneb(i,k)=0.
@@ -398 +412 @@
                zqn(i)=zq(i)*zpdf_e2(i)/zpdf_e1(i)
             endif
+            endif
             
            enddo 
 
-         endif ! iflag_cldcon
 
         endif ! iflag_pdf