Index: trunk/LMDZ.GENERIC/libf/phystd/turb_mod.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/turb_mod.F90	(revision 1834)
+++ trunk/LMDZ.GENERIC/libf/phystd/turb_mod.F90	(revision 1834)
@@ -0,0 +1,35 @@
+module turb_mod
+
+  !! variables
+  REAL,SAVE,ALLOCATABLE :: q2(:,:)    ! Turbulent Kinetic Energy
+  REAL,allocatable,SAVE :: l0(:)
+  REAL,SAVE,ALLOCATABLE :: ustar(:)
+  REAL,SAVE,ALLOCATABLE :: wstar(:)
+  REAL,SAVE,ALLOCATABLE :: tstar(:)
+  REAL,SAVE,ALLOCATABLE :: hfmax_th(:)
+  REAL,SAVE,ALLOCATABLE :: zmax_th(:)
+  REAL,SAVE,ALLOCATABLE :: sensibFlux(:)
+  LOGICAL :: turb_resolved = .false.
+      ! this is a flag to say 'turbulence is resolved'
+      ! mostly for LES use. default is FALSE (for GCM and mesoscale)
+
+!contains
+
+!  subroutine ini_turb_mod(ngrid,nlayer)
+
+!  implicit none
+!  integer,intent(in) :: ngrid ! number of atmospheric columns
+!  integer,intent(in) :: nlayer ! number of atmospheric layers
+
+!    allocate(q2(ngrid,nlayer+1))
+!    allocate(l0(ngrid))
+!    allocate(wstar(ngrid))
+!   allocate(ustar(ngrid))
+!   allocate(tstar(ngrid))
+!   allocate(hfmax_th(ngrid))
+!   allocate(zmax_th(ngrid))
+!   allocate(sensibFlux(ngrid))
+
+!  end subroutine ini_turb_mod
+
+end module turb_mod
Index: trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90	(revision 1833)
+++ trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90	(revision 1834)
@@ -13,5 +13,8 @@
       use comcstfi_mod, only: rcp, g, r, cpp
       use callkeys_mod, only: water,tracer,nosurf
-
+#ifdef MESOSCALE
+      use turb_mod, only : ustar
+      use comm_wrf, only : comm_LATENT_HF
+#endif
       implicit none
 
@@ -259,5 +262,4 @@
       ENDDO
 !we treat only winds, energy and tracers coefficients will be computed with upadted winds
- 
 !JL12 calculate the flux coefficients (tables multiplied elements by elements)
       zfluxv(1:ngrid,1:nlay)=zkv(1:ngrid,1:nlay)*zb0(1:ngrid,1:nlay)
@@ -359,5 +361,7 @@
          zcdh(ig) = zcdh_true(ig)*sqrt(zu2)
          zkh(ig,1)= zcdh(ig)
-      ENDDO
+         ustar(ig)=sqrt(zcdv_true(ig))*sqrt(zu2)
+      ENDDO
+
 
 !     JL12 calculate the flux coefficients (tables multiplied elements by elements)
@@ -719,7 +723,14 @@
 
       if(water)then
+#ifndef MESOSCALE
          call writediagfi(ngrid,'beta','Dryness coefficient',' ',2,dryness)
+#endif
          if (tracer) then
+#ifndef MESOSCALE
             call writediagfi(ngrid,'dqevap','evaporated water vapor specific concentration','s-1',3,pdqevap)
+#else
+            comm_LATENT_HF(:)=0.0
+            comm_LATENT_HF(1:ngrid)=RLVTT*dqsdif_total(1:ngrid)/ptimestep
+#endif
 	 endif
       endif