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

Timestamp:

Apr 28, 2010, 4:27:15 PM (14 years ago)

Author:

musat

Message:

Some changes for new physics :

• new formulation for wakes' alp_offset (physiq)
• bug fix for lalim in thermals plume
• add iflag_pbl=9 option (yamada4)

FH/IM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/yamada4.F (modified) (2 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