Changeset 2420 for LMDZ5/trunk/libf/phylmd/cv3_routines.F90

Timestamp:

Jan 5, 2016, 4:37:49 PM (8 years ago)

Author:

crio

Message:

Nouvelle option d'epluchage de l'ascendance adiabatique dans le schema d'Emanuel: epluchage fonction de B/w2 au lieu de w. S'active avec iflag_mix_adiab=1 (valeur par defaut iflag_mix_adiab=0). Fonctionne avec iflag_mix=0 et iflag_mix=1.
Correction de bugs dans le schema de convection pour le calcul de inb, cape et buoy (sous le meme flag pour l'instant).
New option for the erosion of the adiabatic ascent in the Emanuel scheme: erosion function of B/w2 instead of w. Activated by iflag_mix_adiab=1 (standard value iflag_mix_adiab=0). Should work with iflag_mix=0 and iflag_mix=1.
Various bug corrections in the convection scheme for the computation of inb, cape, buoy (protected by the same flag for now).

File:

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

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -1087 +1087 @@
 SUBROUTINE cv3_undilute2(nloc, ncum, nd, icb, icbs, nk, &
                          tnk, qnk, gznk, hnk, t, q, qs, gz, &
-                         p, h, tv, lv, lf, pbase, buoybase, plcl, &
+                         p, ph, h, tv, lv, lf, pbase, buoybase, plcl, &
                          inb, tp, tvp, clw, hp, ep, sigp, buoy, frac)
   IMPLICIT NONE
@@ -1113 +1113 @@
   include "conema3.h"
   include "cvflag.h"
+  include "YOMCST2.h"
 
 !inputs:
@@ -1119 +1120 @@
   REAL, DIMENSION (nloc, nd), INTENT (IN)            :: t, q, qs, gz
   REAL, DIMENSION (nloc, nd), INTENT (IN)            :: p
+  REAL, DIMENSION (nloc, nd+1), INTENT (IN)          :: ph
   REAL, DIMENSION (nloc), INTENT (IN)                :: tnk, qnk, gznk
   REAL, DIMENSION (nloc), INTENT (IN)                :: hnk
@@ -1144 +1146 @@
   INTEGER iposit(nloc)
   REAL fracg
+  REAL deltap
 
 ! =====================================================================
@@ -1476 +1479 @@
     END DO
   END DO
+
+!CR fix computation of inb
+!keep flag or modify in all cases?
+  IF (iflag_mix_adiab.eq.1) THEN
+  DO i = 1, ncum
+     cape(i)=0.
+     inb(i)=icb(i)+1
+  ENDDO
+  
+  DO k = 2, nl 
+    DO i = 1, ncum
+       IF ((k>=iposit(i))) THEN
+       deltap = min(plcl(i), ph(i,k-1)) - min(plcl(i), ph(i,k))
+       cape(i) = cape(i) + rrd*buoy(i, k-1)*deltap/p(i, k-1)
+       IF (cape(i).gt.0.) THEN
+        inb(i) = max(inb(i), k)
+       END IF
+       ENDIF
+    ENDDO
+  ENDDO
+
+!  DO i = 1, ncum
+!     print*,"inb",inb(i)
+!  ENDDO
+
+  endif
 
 ! -- end convect3