Index: trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 334)
+++ trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 336)
@@ -100,6 +100,6 @@
          r_aspect_thermals=0.7
 #else
-         nsplit_thermals=25
-         r_aspect_thermals=1.8
+         nsplit_thermals=35
+         r_aspect_thermals=1.5
 #endif
 
Index: trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 334)
+++ trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 336)
@@ -182,4 +182,7 @@
 !   ---------------
 
+      entr(:,:)=0.
+      detr(:,:)=0.
+      fm(:,:)=0.
       zu(:,:)=pu(:,:)
       zv(:,:)=pv(:,:)
@@ -927,4 +930,5 @@
                detr(ig,l)=detr(ig,l)+fm(ig,l+1)
                fm(ig,l+1)=0.
+               entr(ig,l+1)=0.
                ncorecfm2=ncorecfm2+1
                else
@@ -988,6 +992,6 @@
 !            endif
 
-!            if(l.gt.lmax(ig)) then
-            if(l.gt.lalim(ig)) then
+            if(l.gt.lmax(ig)) then
+!            if(l.gt.lalim(ig)) then
                detr(ig,l)=0.
                fm(ig,l+1)=0.