Index: LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90	(revision 4713)
+++ LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90	(revision 4714)
@@ -20,6 +20,5 @@
 !
 ! Main assumptions of the model :
-! (1) dry atmosphere
-! (2) horizontal homogeneity (Dx=Dy=0.)
+! (1) horizontal homogeneity (Dx=Dy=0.)
 !=======================================================================
 
@@ -205,4 +204,18 @@
    ENDDO
 
+ELSE IF (iflag_atke_lmix .EQ. 3 ) THEN
+! add effect of wind shear on lstrat following grisogono 2010, qjrms
+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. .AND. shear2(igrid,ilay) .GT. 0.) THEN
+             lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay))*(1.0+Ri(igrid,ilay)/(2.*Prandtl(igrid,ilay)))
+             lstrat=max(lstrat,lmin)
+             !Inverse interpolation, Van de Wiel et al. 2010   
+             l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0)
+          ENDIF
+      ENDDO
+   ENDDO
+