Index: trunk/LMDZ.GENERIC/libf/phystd/comm_wrf.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/comm_wrf.F90	(revision 2010)
+++ trunk/LMDZ.GENERIC/libf/phystd/comm_wrf.F90	(revision 2019)
@@ -32,4 +32,6 @@
   REAL,SAVE,ALLOCATABLE :: comm_H2OICE_REFF(:,:)
   REAL,SAVE,ALLOCATABLE :: comm_LATENT_HF(:)
+  REAL,SAVE,ALLOCATABLE :: comm_REEVAP(:)
+  REAL,SAVE,ALLOCATABLE :: comm_SURFRAIN(:)
 
 contains
@@ -66,4 +68,6 @@
   allocate(comm_H2OICE_REFF(ngrid,nlayer))
   allocate(comm_LATENT_HF(ngrid))
+  allocate(comm_REEVAP(ngrid))
+  allocate(comm_SURFRAIN(ngrid)
 
   end subroutine allocate_comm_wrf
@@ -98,4 +102,6 @@
   deallocate(comm_H2OICE_REFF)
   deallocate(comm_LATENT_HF)
+  deallocate(comm_REEVAP)
+  deallocate(comm_SURFRAIN)
 
   end subroutine deallocate_comm_wrf
Index: trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90	(revision 2010)
+++ trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90	(revision 2019)
@@ -392,5 +392,8 @@
       real :: flux_sens_lat(ngrid)
       real :: qsurfint(ngrid,nq)
-
+#ifdef MESOSCALE 
+      REAL :: lsf_dt(nlayer)
+      REAL :: lsf_dq(nlayer)
+#endif
 
 !==================================================================================================
@@ -2157,9 +2160,14 @@
       comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer)
       comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer)
-      comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer)
-      comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid)
-      comm_RAIN(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap)
-      comm_SNOW(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_ice)
-      comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid)
+      if ((tracer).and.(water)) then
+        comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer)
+        comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid)
+        comm_RAIN(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap)
+        comm_SURFRAIN(1:ngrid)=zdqsrain(1:ngrid)
+        comm_SNOW(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_ice)
+        comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid)
+        comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o)
+        comm_REEVAP(1:ngrid)=reevap_precip(1:ngrid)
+      endif
       comm_FLUXTOP_DN(1:ngrid)=fluxtop_dn(1:ngrid)
       comm_FLUXABS_SW(1:ngrid)=fluxabs_sw(1:ngrid)
@@ -2169,6 +2177,21 @@
       comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid)
       comm_LSCEZ(1:ngrid,1:nlayer)=lscaledEz(1:ngrid,1:nlayer)
-      comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o)
       sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) !!! ????
+
+      if (turb_resolved) then
+        open(17,file='lsf.txt',form='formatted',status='old')
+        rewind(17)
+        DO l=1,nlayer
+          read(17,*) lsf_dt(l),lsf_dq(l)
+        ENDDO
+        close(17)
+        do ig=1,ngrid
+          if ((tracer).and.(water)) then
+           pdq(ig,:,igcm_h2o_vap) = pdq(ig,:,igcm_h2o_vap) + lsf_dq(:)
+          endif
+          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)
+          comm_HR_DYN(ig,:) = lsf_dt(:)
+        enddo
+      endif
 #endif