Changeset 2757 for LMDZ5/trunk/libf/phylmd/wake.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/wake.F90 (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/wake.F90

@@ -171 +171 @@
   !$OMP THREADPRIVATE(stark, wdens_ref, coefgw, alpk, crep_upper, crep_sol)
 
+  LOGICAL, SAVE                                         :: flag_wk_check_trgl
+  !$OMP THREADPRIVATE(flag_wk_check_trgl)
+
   REAL                                                  :: delta_t_min
   INTEGER                                               :: nsub
@@ -200 +203 @@
   REAL, DIMENSION (klon)                                :: z, dz, hw0
   INTEGER, DIMENSION (klon)                             :: ktop, kupper
+
+  ! Variables liées au test de la forme triangulaire du profil de Delta_theta
+  REAL, DIMENSION (klon)                                :: sum_half_dth
+  REAL, DIMENSION (klon)                                :: dz_half
 
   ! Sub-timestep tendencies and related variables
@@ -321 +328 @@
   WRITE(*,*) 'wdens_ref=', wdens_ref
   WRITE(*,*) 'coefgw=', coefgw 
+
+  flag_wk_check_trgl=.false.
+  CALL getin_p('flag_wk_check_trgl ', flag_wk_check_trgl)
+  WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
 
   first=.false.
@@ -824 +835 @@
 
   IF (prt_level>=10) THEN
-    PRINT *, 'wake-4, sigmaw(igout), cstar(igout), wape(igout) ', &
-                      sigmaw(igout), cstar(igout), wape(igout)
+    PRINT *, 'wake-4, sigmaw(igout), cstar(igout), wape(igout), ktop(igout) ', &
+                      sigmaw(igout), cstar(igout), wape(igout), ktop(igout)
   ENDIF
 
@@ -846 +857 @@
     END DO
     IF (prt_level>=10) THEN
-      PRINT *, 'wake-4.1, isubstep,wk_adv(igout),cstar(igout),wape(igout) ', &
-                          isubstep,wk_adv(igout),cstar(igout),wape(igout)
+      PRINT *, 'wake-4.1, isubstep,wk_adv(igout),cstar(igout),wape(igout), ptop(igout) ', &
+                          isubstep,wk_adv(igout),cstar(igout),wape(igout), ptop(igout)
     ENDIF
 
@@ -962 +973 @@
     IF (prt_level>=10) THEN
       PRINT *, 'wake-4.2, omg(igout,k) ', (k,omg(igout,k), k=1,klev)
-      PRINT *, 'wake-4.2, omgtop(igout) ', omgtop(igout)
+      PRINT *, 'wake-4.2, omgtop(igout), ptop(igout), ktop(igout) ', &
+                          omgtop(igout), ptop(igout), ktop(igout)
     ENDIF
 
@@ -1609 +1621 @@
 
   IF (prt_level>=10) THEN
-    PRINT *, 'wake-5, sigmaw(igout), cstar(igout), wape(igout) ', &
-                      sigmaw(igout), cstar(igout), wape(igout)
+    PRINT *, 'wake-5, sigmaw(igout), cstar(igout), wape(igout), ptop(igout) ', &
+                      sigmaw(igout), cstar(igout), wape(igout), ptop(igout)
   ENDIF
 
@@ -1632 +1644 @@
       sum_thvu(i) = 0.
       sum_dth(i) = 0.
+      sum_half_dth(i) = 0.
       sum_dq(i) = 0.
       sum_rho(i) = 0.
@@ -1646 +1659 @@
       av_dtdwn(i) = 0.
       av_dqdwn(i) = 0.
+
+      dthmin(i) = -delta_t_min
     END IF
   END DO
@@ -1685 +1700 @@
       z(i) = 1.
       dz(i) = 1.
+      dz_half(i) = 1.
       sum_thvu(i) = thu(i, 1)*(1.+epsim1*qu(i,1))*dz(i)
       sum_dth(i) = 0.
@@ -1696 +1712 @@
         ! cc
         dz(i) = -(amax1(ph(i,k+1),ptop(i))-ph(i,k))/(rho(i,k)*rg)
+        dz_half(i) = -(amax1(ph(i,k+1),0.5*(ptop(i)+ph(i,1)))-ph(i,k))/(rho(i,k)*rg)
         IF (dz(i)>0) THEN
           z(i) = z(i) + dz(i)
@@ -1707 +1724 @@
           sum_dtdwn(i) = sum_dtdwn(i) + dtdwn(i, k)*dz(i)
           sum_dqdwn(i) = sum_dqdwn(i) + dqdwn(i, k)*dz(i)
+!
+          dthmin(i) = min(dthmin(i), dth(i,k))
+        END IF
+        IF (dz_half(i)>0) THEN
+          sum_half_dth(i) = sum_half_dth(i) + dth(i, k)*dz_half(i)
         END IF
       END IF
@@ -1750 +1772 @@
 
   ! Filter out bad wakes
+
+  IF (flag_wk_check_trgl) THEN
+    ! Check triangular shape of dth profile
+    DO i = 1, klon
+      IF (ok_qx_qw(i)) THEN
+        !! print *,'wake, hw0(i), dthmin(i) ', hw0(i), dthmin(i)
+        !! print *,'wake, 2.*sum_dth(i)/(hw0(i)*dthmin(i)) ', &
+        !!                2.*sum_dth(i)/(hw0(i)*dthmin(i))
+        !! print *,'wake, sum_half_dth(i), sum_dth(i) ', &
+        !!                sum_half_dth(i), sum_dth(i)
+        IF ((hw0(i) < 1.) .or. (dthmin(i) >= -delta_t_min) ) THEN
+          wape2(i) = -1.
+          !! print *,'wake, rej 1'
+        ELSE IF (abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN
+          wape2(i) = -1.
+          !! print *,'wake, rej 2'
+        ELSE IF (abs(sum_half_dth(i)) < 0.5*abs(sum_dth(i)) ) THEN
+          wape2(i) = -1.
+          !! print *,'wake, rej 3'
+        END IF
+      END IF
+    END DO
+  END IF
+
 
   DO k = 1, klev