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

Timestamp:

Jan 10, 2017, 3:30:46 PM (8 years ago)

Author:

jyg

Message:

Changes to convection and wakes:
+ convection is allowed only when

Plcl < 0.99*psurf

Controlled by cv_flag_feed=2 (D:1)

+ check on the triangular shape of
DT_wk.

Controlled by flag_wk_check_trgl=y
(D : n)

File:

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

LMDZ5/trunk/libf/phylmd/cv3_routines.F90

@@ -83 +83 @@
 
 ! Var interm pour le getin
+     cv_flag_feed=1
+     CALL getin_p('cv_flag_feed',cv_flag_feed)
      dpbase=-40.
      CALL getin_p('dpbase',dpbase)
@@ -453 +455 @@
     END DO
   END DO !  iter
+
   DO i = 1, len
     p2feed(i) = pfeed(i)
@@ -467 +470 @@
 ! --- are reasonable
 ! -------------------------------------------------------------------
-  DO i = 1, len
-    IF (((tnk(i)<250.0) .OR. (qnk(i)<=0.0)) .AND. (iflag(i)==0)) iflag(i) = 7
-  END DO
+  IF (cv_flag_feed == 1) THEN
+    DO i = 1, len
+      IF (((tnk(i)<250.0)                       .OR.  &
+           (qnk(i)<=0.0))                       .AND. &
+          (iflag(i)==0)) iflag(i) = 7
+    END DO
+  ELSEIF (cv_flag_feed >= 2) THEN
+! --- and demand that LCL be high enough
+    DO i = 1, len
+      IF (((tnk(i)<250.0)                       .OR.  &
+           (qnk(i)<=0.0)                        .OR.  &
+           (plcl(i)>min(0.99*ph(i,1),ph(i,3)))) .AND. &
+          (iflag(i)==0)) iflag(i) = 7
+    END DO
+  ENDIF
 
 ! -------------------------------------------------------------------