Index: LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90	(revision 4631)
+++ LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90	(revision 4632)
@@ -176,6 +176,22 @@
    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,5 +233,5 @@
 
 ! 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
Index: LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90	(revision 4631)
+++ LMDZ6/trunk/libf/phylmd/atke_turbulence_ini_mod.F90	(revision 4632)
@@ -97,5 +97,5 @@
 
   ! coefficient for mixing length depending on local stratification
-  clmix=0.2
+  clmix=0.5
   CALL getin_p('atke_clmix',clmix)