Changeset 1518 for LMDZ5/trunk/libf/phylmd/cv3p1_closure.F

Timestamp:

May 12, 2011, 3:32:22 PM (13 years ago)

Author:

idelkadi

Message:

Modifications des routines de convection :

• correction de bug : le champ sigd n'etait pas decompresse avant d'etre sorti de la convection.
• sortir de nouveaux champs de convection (plcl, plfc, wb)

File:

• LMDZ5/trunk/libf/phylmd/cv3p1_closure.F (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/cv3p1_closure.F

@@ -7 +7 @@
      o                      ,sig,w0,ptop2,cape,cin,m,iflag,coef
      :                      ,Plim1,Plim2,asupmax,supmax0
-     :                      ,asupmaxmin,cbmf)
+     :                      ,asupmaxmin,cbmf,plfc,wbeff)
 
 *
@@ -49 +49 @@
       real asupmax(nloc,nd),supmax0(nloc)
       real asupmaxmin(nloc)
+      real cbmf(nloc),plfc(nloc)
+      real wbeff(nloc)
       integer iflag(nloc)
 c
@@ -63 +65 @@
       integer nsupmax(nloc)
       real supcrit,temp(nloc,nd)
-      real P1(nloc),Pmin(nloc),plfc(nloc)
+      real P1(nloc),Pmin(nloc)
       real asupmax0(nloc)
       logical ok(nloc)
       real siglim(nloc,nd),wlim(nloc,nd),mlim(nloc,nd)
       real wb2(nloc)
-      real cbmflim(nloc),cbmf1(nloc),cbmfmax(nloc),cbmf(nloc)
+      real cbmflim(nloc),cbmf1(nloc),cbmfmax(nloc)
       real cbmflast(nloc)
-      real wbeff(nloc)
-
       real coef(nloc)
       real xp(nloc),xq(nloc),xr(nloc),discr(nloc),b3(nloc),b4(nloc)
@@ -203 +203 @@
            nsupmax(il)=i
            ok(il)=.false.
-        ENDIF    ! end IF (P(i) ...
+        ENDIF    ! end IF (P(i) ...  )
         ENDIF    ! end IF (icb+1 le i le inb)
         ENDDO
@@ -308 +308 @@
      $             /(P(il,i)-P(il,i-1))
          GO TO 425
-        ENDIF    ! end IF (P(i) ...
+        ENDIF    ! end IF (P(i) ... )
         ENDIF    ! end IF (icb+1 le i le inb)
         ENDDO
@@ -517 +517 @@
       ENDDO
 c
-      if (flag_wb==0) then
-         wbeff(:)=wbmax
-      else if (flag_wb==1) then
-         wbeff(1:ncum)=wbmax/(1.+500./(ph(1:ncum,1)-plfc(1:ncum)))
-      else
-         stop'cas non prevu iflag_wb'
-! ligne a suprimer un jour ...
-      endif
-!     if (ncum==1) print*,'WB ',plfc,wbeff,alp2,cin
-
+      IF (flag_wb==0) THEN
+        wbeff(:) = wbmax
+      ELSE IF (flag_wb==1) THEN
+        wbeff(1:ncum) = wbmax/(1.+500./(ph(1:ncum,1)-plfc(1:ncum)))
+      ELSE IF (flag_wb==2) THEN
+        wbeff(1:ncum) = wbmax*(0.01*(ph(1:ncum,1)-plfc(1:ncum)))**2
+      ENDIF
+c
       DO il = 1,ncum
-
-!jyg    Modification du coef de wb*wb pour conformite avec papier Wake
-!       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
-! A partir de mai 2011 on revient au 2 wb^2, conformement au papier.
-! les versions d'avant (NPv2 et al) tournaient avec 0.5 wb^2 et wb=1m/s
-! equivalent  a wb=1 m/s pour la nouvelle formulation
-! Puis passage de wb constant a wbeff variable suivant flag_wb
-
+cjyg    Modification du coef de wb*wb pour conformite avec papier Wake
+cc       cbmf1(il) = alp2(il)/(0.5*wb*wb-Cin(il))
        cbmf1(il) = alp2(il)/(2.*wbeff(il)*wbeff(il)-Cin(il))
-
        if(cbmf1(il).EQ.0.AND.alp2(il).NE.0.) THEN
         write(lunout,*)