Changeset 2294 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Jun 11, 2015, 8:06:49 PM (9 years ago)

Author:

fhourdin

Message:

Introduction d'un seuil sur la probabilité de non déclenchement
stochastique random_notrig_max
Si le tirage aleatoire uniforme entre 0 et 1 est > random_notrig_max
on ne déclenche pas.
random_notrig_max=1-epsilon avec epsilon petit.

Corrections dans grid_atob_m.F90 pour la compilation Linux.
A vérifier.

File:

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

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -60 +60 @@
   USE YOERAD   , ONLY : NRADLP
 #endif
+  USE ioipsl_getin_p_mod, ONLY : getin_p
+
 
   !IM stations CFMIP
@@ -440 +442 @@
   integer :: tau_trig(klon)
 
+  REAL,SAVE :: random_notrig_max=1.01
+  !$OMP THREADPRIVATE(random_notrig_max) 
+
   !--------Statistical Boundary Layer Closure: ALP_BL--------
   !---Profils de TKE dans et hors du thermique
@@ -993 +998 @@
   IF (debut) THEN
      CALL suphel ! initialiser constantes et parametres phys.
+     CALL getin_p('random_notrig_max',random_notrig_max)
   ENDIF
 
@@ -2734 +2740 @@
               proba_notrig(i)=1.
               random_notrig(i)=1e6*ale_bl_stat(i)-int(1e6*ale_bl_stat(i))
+              if ( random_notrig(i) > random_notrig_max ) random_notrig(i)=1.1
               if ( ale_bl_trig(i) .lt. abs(cin(i))+1.e-10 ) then 
                  tau_trig(i)=tau_trig_shallow