Changeset 5172 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_atke_turbulence_ini.F90

Timestamp:

Sep 9, 2024, 10:58:36 AM (2 months ago)

Author:

abarral

Message:

Merge r5164, r5166, r5167, r5168

Location:

LMDZ6/branches/Amaury_dev

Files:

• . (modified) (1 prop)
• libf/phylmd/lmdz_atke_turbulence_ini.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/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
     !!----------------------------------------------------------------------
 
@@ -74 +75 @@
 
     ! 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 @@
     endif
 
-    ! activate vertical diffusion of TKE or not
-    atke_ok_vdiff = .FALSE.
-    CALL getin_p('atke_ok_vdiff', atke_ok_vdiff)
+      ! activate vertical diffusion of TKE or not
+      atke_ok_vdiff = .TRUE.
+      CALL getin_p('atke_ok_vdiff', atke_ok_vdiff)
 
     ! account for vapor for flottability
@@ -102 +103 @@
     ! 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)
 
@@ -132 +133 @@
 
     ! 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 <= 1) THEN
@@ -140 +141 @@
 
     ! 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)
 
@@ -152 +153 @@
 
     ! coefficient for mixing length depending on local stratification
-    clmix = 0.5
+    clmix = 0.648055235325291
     CALL getin_p('atke_clmix', clmix)
 
@@ -161 +162 @@
     ! 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)