Changeset 1371 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/yamada4.F

Timestamp:

Apr 23, 2010, 4:09:31 PM (14 years ago)

Author:

idelkadi

Message:

• Corrections d'erreurs
• Declarations de variables additionnelles pour les sorties 1D

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/yamada4.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/yamada4.F

@@ -366 +366 @@
 c     print*,'0L=',k,l(ig,k),delta(ig,k),km(ig,k)
          qpre=sqrt(q2(ig,k))
-         if (iflag_pbl.eq.8 ) then
+!        if (iflag_pbl.eq.8 ) then
             if (aa(ig,k).gt.0.) then
                q2(ig,k)=(qpre+aa(ig,k)*qpre*qpre)**2
@@ -372 +372 @@
                q2(ig,k)=(qpre/(1.-aa(ig,k)*qpre))**2
             endif
-         else ! iflag_pbl=9
-            if (aa(ig,k)*qpre.gt.0.9) then
-               q2(ig,k)=(qpre*10.)**2
-            else
-               q2(ig,k)=(qpre/(1.-aa(ig,k)*qpre))**2
-            endif
-         endif
+!        else ! iflag_pbl=9
+!           if (aa(ig,k)*qpre.gt.0.9) then
+!              q2(ig,k)=(qpre*10.)**2
+!           else
+!              q2(ig,k)=(qpre/(1.-aa(ig,k)*qpre))**2
+!           endif
+!        endif
          q2(ig,k)=min(max(q2(ig,k),1.e-10),1.e4)
 c     print*,'Q2 L=',k,q2(ig,k),qpre*qpre
@@ -403 +403 @@
       enddo
 
+! Transport diffusif vertical de la TKE.
       if (iflag_pbl.ge.9) then
 !       print*,'YAMADA VDIF'