Changeset 5164 for LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90

Timestamp:

Aug 26, 2024, 2:57:30 PM (3 months ago)

Author:

evignon

Message:

mise par defaut des parametres de la param ATKE suivant la BEST du papier Vignon et al. 2024

File:

• LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90

@@ -50 +50 @@
       !!
       !! ** Purpose :   Initialization of the atke module and choice of some constants
-      !!               
+      !!                Default values correspond to the  'best' configuration
+      !!                from tuning on GABLS1 in Vignon et al. 2024, JAMES
       !!----------------------------------------------------------------------
 
@@ -73 +74 @@
 
       ! flag that controls options in atke_compute_km_kh
-      iflag_atke=0
+      iflag_atke=1
       CALL getin_p('iflag_atke',iflag_atke)
 
       ! flag that controls the calculation of mixing length in atke
-      iflag_atke_lmix=0
+      iflag_atke_lmix=3
       CALL getin_p('iflag_atke_lmix',iflag_atke_lmix)
 
@@ -86 +87 @@
 
       ! activate vertical diffusion of TKE or not
-      atke_ok_vdiff=.false.
+      atke_ok_vdiff=.true.
       CALL getin_p('atke_ok_vdiff',atke_ok_vdiff)
 
@@ -101 +102 @@
       ! Sun et al 2011, JAMC 
       ! between 10 and 40
-      l0=15.0 
+      l0=42.5279652116005 
       CALL getin_p('atke_l0',l0)
 
       ! critical Richardson number 
-      ric=0.25
+      ric=0.190537327781655
       CALL getin_p('atke_ric',ric)
 
       ! constant for tke dissipation calculation
-      cepsilon=5.87 ! default value as in yamada4
+      cepsilon=8.89273387537601
       CALL getin_p('atke_cepsilon',cepsilon)
 
@@ -131 +132 @@
 
       ! slope of Pr=f(Ri) for stable conditions
-      pr_slope=5.0 ! default value from Zilitinkevich et al. 2005
+      pr_slope=4.67885738180385
       CALL getin_p('atke_pr_slope',pr_slope)
       if (pr_slope .le. 1) then
@@ -139 +140 @@
 
       ! value of turbulent prandtl number in neutral conditions (Ri=0)
-      pr_neut=0.8
+      pr_neut=0.837372701768868
       CALL getin_p('atke_pr_neut',pr_neut)
 
@@ -151 +152 @@
 
       ! coefficient for mixing length depending on local stratification
-      clmix=0.5
+      clmix=0.648055235325291
       CALL getin_p('atke_clmix',clmix)
 
@@ -160 +161 @@
       ! minimum anisotropy coefficient (defined here as minsqrt(Ez/Ek)) at large Ri. 
       ! From Zilitinkevich et al. 2013, it equals sqrt(0.03)~0.17  
-      smmin=0.17
+      smmin=0.0960838631869678
       CALL getin_p('atke_smmin',smmin)
 
       ! ratio between the eddy diffusivity coeff for tke wrt that for momentum
       ! default value from Lenderink et al. 2004
-      cke=2.
+      cke=2.47069655134662
       CALL getin_p('atke_cke',cke)