Changeset 4777 for LMDZ6/trunk/libf/phylmd

Timestamp:

Dec 18, 2023, 11:03:57 PM (12 months ago)

Author:

evignon

Message:

petites corrections dans les routines atke

Location:

LMDZ6/trunk/libf/phylmd

Files:

• cdrag_mod.F90 (modified) (2 diffs)
• lmdz_atke_exchange_coeff.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/cdrag_mod.F90

@@ -26 +26 @@
   USE print_control_mod, ONLY: lunout, prt_level
   USE ioipsl_getin_p_mod, ONLY : getin_p
-  USE lmdz_atke_turbulence_ini, ONLY : ric, cinf, cepsilon, pr_slope, pr_asym, pr_neut
+  USE lmdz_atke_turbulence_ini, ONLY : smmin, ric, cinf, cepsilon, pr_slope, pr_asym, pr_neut
 
   IMPLICIT NONE
@@ -472 +472 @@
                   endif
              CASE (6) ! Consistent with turbulence scheme (in stationary case) derived in atke (2023)
-                  sm(i) = MAX(0.,cn*(1.-zri(i)/ric))
-                  prandtl(i) = pr_neut + zri(i) * pr_slope
+                  sm(i) = MAX(smmin,cn*(1.-zri(i)/ric))
+                           ! prandlt expression from venayagamoorthy and stretch 2010, Li et al 2019
+                  prandtl(i) = pr_neut*exp(-pr_slope/pr_neut*zri(i)+zri(i)/pr_neut) &
+                                + zri(i) * pr_slope
                   FM(i) = MAX((sm(i)**(3./2.) * sqrt(cepsilon) * (1 - zri(i) / prandtl(i))**(1./2.)),f_ri_cd_min)
                   FH(i) = MAX((FM(i) / prandtl(i)), f_ri_cd_min)

LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90

@@ -191 +191 @@
 ELSE IF (iflag_atke_lmix .EQ. 2 ) THEN
 ! add effect of wind shear on lstrat following grisogono and belusic 2008, qjrms
+! implies 2 tuning coefficients clmix and clmixshear
 DO ilay=2,nlay
       DO igrid=1,ngrid
@@ -206 +207 @@
 ELSE IF (iflag_atke_lmix .EQ. 3 ) THEN
 ! add effect of wind shear on lstrat following grisogono 2010, qjrms
+! keeping a single tuning coefficient clmix
 DO ilay=2,nlay
       DO igrid=1,ngrid