Index: trunk/LMDZ.MARS/libf/phymars/phys_state_var_init_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/phys_state_var_init_mod.F90	(revision 4008)
+++ trunk/LMDZ.MARS/libf/phymars/phys_state_var_init_mod.F90	(revision 4009)
@@ -88,5 +88,4 @@
 #ifndef MESOSCALE
       if(is_master) then
-#endif
         status = nf90_open(filename, nf90_nowrite, ncid)
         if (status /= nf90_noerr) then
@@ -105,7 +104,9 @@
           endif
         endif
-#ifndef MESOSCALE
       endif
       call bcast(nslope)
+#else
+      print *, "subgrid slopes temporarily switched off in mesoscale, nslope=1"
+      nslope=1
 #endif
 
Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 4008)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 4009)
@@ -3172,5 +3172,5 @@
       comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer)
       !state  real  SWDOWNZ    ij   misc  1  -  h  "SWDOWNZ"   "DOWNWARD SW FLUX AT SURFACE"     "W m-2"
-      comm_SWDOWNZ(1:ngrid) = fluxsurf_dn_sw_tot(1:ngrid)
+      comm_SWDOWNZ(1:ngrid) = fluxsurf_dn_sw_tot(1:ngrid,1)
       !state  real  TAU_DUST   ij   misc  1  -  h  "TAU_DUST"  "REFERENCE VISIBLE DUST OPACITY"  ""
       comm_TAU_DUST(1:ngrid) = tau_pref_gcm(1:ngrid)
@@ -3179,5 +3179,5 @@
       !state  real  QSURFDUST  ij   misc  1  -  h  "QSURFDUST" "DUST MASS AT SURFACE"            "kg m-2"
       IF (igcm_dust_mass .ne. 0) THEN
-        comm_QSURFDUST(1:ngrid) = qsurf(1:ngrid,igcm_dust_mass)
+        comm_QSURFDUST(1:ngrid) = qsurf(1:ngrid,igcm_dust_mass,1)
       ELSE
         comm_QSURFDUST(1:ngrid) = 0.
@@ -3208,6 +3208,6 @@
       !! -- vdifc must have been called
       if (.not.callrichsl) then
-        sensibFlux(1:ngrid) = zflubid(1:ngrid)
-     .         - capcal(1:ngrid)*zdtsdif(1:ngrid)
+        sensibFlux(1:ngrid) = zflubid(1:ngrid,1)
+     .         - capcal(1:ngrid,1)*zdtsdif(1:ngrid,1)
       else
         sensibFlux(1:ngrid) =
@@ -3215,5 +3215,5 @@
      &   *sqrt(pu(1:ngrid,1)*pu(1:ngrid,1)+pv(1:ngrid,1)*pv(1:ngrid,1)
      &         +(log(1.+0.7*wstar(1:ngrid) + 2.3*wstar(1:ngrid)**2))**2)
-     &   *zcdh(1:ngrid)*(tsurf(1:ngrid)-zh(1:ngrid,1))
+     &   *zcdh(1:ngrid,1)*(tsurf(1:ngrid,1)-zh(1:ngrid,1))
       endif