Changeset 2922 for LMDZ5/trunk/libf

Timestamp:

Jun 29, 2017, 5:45:27 PM (8 years ago)

Author:

fhourdin

Message:

Nouvelle option iflag_wk_check_trgl pour tester la bonne sante des poches.
L'ancien flag logique flag_wk_check_trgl est toujours utilisable.
Il vaut y/n si iflag_wk_check_trgl=1/0.
L'option 1 avait tendance à faire disparaitre les poches dans des cas 1D.
L'option 2 est donc précaunisée pour les réglages à venir.
Fredho pour Jean-Yves

File:

• LMDZ5/trunk/libf/phylmd/wake.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/wake.F90

@@ -178 +178 @@
   LOGICAL, SAVE                                         :: flag_wk_check_trgl
   !$OMP THREADPRIVATE(flag_wk_check_trgl)
+  INTEGER, SAVE                                         :: iflag_wk_check_trgl
+  !$OMP THREADPRIVATE(iflag_wk_check_trgl)
 
   REAL                                                  :: delta_t_min
@@ -347 +349 @@
   CALL getin_p('flag_wk_check_trgl ', flag_wk_check_trgl)
   WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
+  WRITE(*,*) 'flag_wk_check_trgl OBSOLETE. Utilisr iflag_wk_check_trgl plutot'
+  iflag_wk_check_trgl=0 ; IF (flag_wk_check_trgl) iflag_wk_check_trgl=1
+  CALL getin_p('iflag_wk_check_trgl ', iflag_wk_check_trgl)
+  WRITE(*,*) 'iflag_wk_check_trgl=', iflag_wk_check_trgl
 
   first=.false.
@@ -1793 +1799 @@
   ! Filter out bad wakes
 
-  IF (flag_wk_check_trgl) THEN
+  IF (iflag_wk_check_trgl>=1) THEN
     ! Check triangular shape of dth profile
     DO i = 1, klon
@@ -1805 +1811 @@
           wape2(i) = -1.
           !! print *,'wake, rej 1'
-        ELSE IF (abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN
+        ELSE IF (iflag_wk_check_trgl==1.AND.abs(2.*sum_dth(i)/(hw0(i)*dthmin(i)) - 1.) > 0.5) THEN
           wape2(i) = -1.
           !! print *,'wake, rej 2'