Index: trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 276)
+++ trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 277)
@@ -95,4 +95,5 @@
 
 #ifdef MESOSCALE
+         !! valid for timesteps < 200s 
          nsplit_thermals=2
          r_aspect_thermals=0.7
Index: trunk/LMDZ.MARS/libf/phymars/physiq.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 276)
+++ trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 277)
@@ -788,7 +788,10 @@
         endif
 
+        lmax_th_out(:)=real(lmax_th(:))
+
         else   !of if calltherm
         lmax_th(:)=0
         wmax_th(:)=0.
+        lmax_th_out(:)=0.
         end if
 
@@ -1631,6 +1634,4 @@
 !        endif
 
-        lmax_th_out(:)=real(lmax_th(:))
-
         call WRITEDIAGFI(ngridmx,'lmax_th',
      &              'hauteur du thermique','K',
@@ -1725,6 +1726,4 @@
          if(calltherm) then
 
-        lmax_th_out(:)=real(lmax_th(:))
-
         call WRITEDIAGFI(ngridmx,'lmax_th',
      &              'hauteur du thermique','point',
@@ -1755,6 +1754,6 @@
          call WRITEDIAGFI(ngrid,"tsurf","Surface temperature","K",0,
      &                  tsurf)
-!         call WRITEDIAGFI(ngrid,"u","u wind","m/s",1,zu)
-!         call WRITEDIAGFI(ngrid,"v","v wind","m/s",1,zv)
+         call WRITEDIAGFI(ngrid,"u","u wind","m/s",1,zu)
+         call WRITEDIAGFI(ngrid,"v","v wind","m/s",1,zv)
 
          call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",1,zplay)