Changeset 2420 for LMDZ5/trunk/libf/phylmd/cv3_buoy.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_buoy.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cv3_buoy.F90

@@ -14 +14 @@
   include "cvthermo.h"
   include "cv3param.h"
+  include "YOMCST2.h"
 
   ! input:
@@ -139 +140 @@
   END DO
 
-
+!CR:Correction of buoy for what comes next
+!keep flag or to modify in all cases?
+  IF (iflag_mix_adiab.eq.1) THEN
+  DO k = 1, nl 
+    DO il = 1, ncum
+       IF ((k>=kmx(il)) .AND. (k<=inb(il)) .AND. (buoy(il,k).lt.0.)) THEN
+          buoy(il,k)=buoy(il,k-1)
+       END IF
+    ENDDO
+  ENDDO
+  ENDIF
 
   RETURN