Changeset 2294 for LMDZ5/trunk/libf

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.

Location:

LMDZ5/trunk/libf

Files:

• dynlonlat_phylonlat/grid_atob_m.F90 (modified) (2 diffs)
• phylmd/physiq.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/dynlonlat_phylonlat/grid_atob_m.F90

@@ -83 +83 @@
     RETURN
   END IF
-  nn=COUNT(found/=0); IF(nn==0) RETURN
+  nn=COUNT(found); IF(nn==0) RETURN
 
 !--- MISSING POINTS ; USE DISTANCE ON THE SPHERE TO FIND NEAREST POINT nr(2)
@@ -136 +136 @@
   DOUBLE PRECISION  :: out(SIZE(x),SIZE(y))
 !-------------------------------------------------------------------------------
-  CALL fine2coarse(xdata,ydata,x,y,out,DBLE(entree))
+!  CALL fine2coarse(xdata,ydata,x,y,out,DBLE(entree))
+  CALL fine2coarse(xdata,ydata,x,y,out,entree)
   sortie=REAL(out)
 

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