Changeset 4632

Timestamp:

Jul 18, 2023, 3:27:56 PM (11 months ago)

Author:

evignon

Message:

petits ajustements suite a la derniere commission

Location:

LMDZ6/trunk/libf/phylmd

Files:

• atke_exchange_coeff_mod.F90 (modified) (2 diffs)
• atke_turbulence_ini_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90

@@ -176 +176 @@
    ENDDO
 
+ELSE IF (iflag_atke_lmix .EQ. 2 ) THEN
+! add effect of wind shear on lstrat following grisogono and belusic 2008, qjrms, eq 2
+DO ilay=2,nlay
+      DO igrid=1,ngrid
+          l_exchange(igrid,ilay) = kappa*l0*z_interf(igrid,ilay) / (kappa*z_interf(igrid,ilay) + l0)
+          IF (N2(igrid,ilay) .GT. 0.) THEN
+             lstrat=clmix*sqrt(tke(igrid,ilay))/(2.*max(sqrt(shear2(igrid,ilay)),1E-10)*(1.+sqrt(Ri(igrid,ilay))/2.))
+             IF (lstrat .LT. l_exchange(igrid,ilay)) THEN
+                l_exchange(igrid,ilay)=max(lstrat,lmin)
+             ENDIF
+          ENDIF
+      ENDDO
+   ENDDO
+
+
+
 ELSE
-! default: neglect effect of local stratification
+! default: neglect effect of local stratification and shear
 
    DO ilay=2,nlay+1
@@ -217 +233 @@
 
 ! full implicit scheme resolved with a second order polynomial equation
-! + exact resolution for very stable cases
+! + exact resolution for very stable cases (iflag_atke_lmix=1)
     DO ilay=2,nlay
         DO igrid=1,ngrid

LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90

@@ -97 +97 @@
 
   ! coefficient for mixing length depending on local stratification
-  clmix=0.2
+  clmix=0.5
   CALL getin_p('atke_clmix',clmix)