Index: trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 313)
+++ trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90	(revision 314)
@@ -100,6 +100,6 @@
          r_aspect_thermals=0.7
 #else
-         nsplit_thermals=20
-         r_aspect_thermals=2.
+         nsplit_thermals=25
+         r_aspect_thermals=1.8
 #endif
 
Index: trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 313)
+++ trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 314)
@@ -886,11 +886,11 @@
 ! autres corrections.
 
-      do l=1,nlayermx
-         do ig=1,ngridmx
-            if (detr(ig,l).gt.fm(ig,l)) then
-               ncorecfm8=ncorecfm8+1
-            endif
-         enddo
-      enddo
+!      do l=1,nlayermx
+!         do ig=1,ngridmx
+!            if (detr(ig,l).gt.fm(ig,l)) then
+!               ncorecfm8=ncorecfm8+1
+!            endif
+!         enddo
+!      enddo
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -986,11 +986,15 @@
 ! Dans le cas ou on est au dessus de la couche d'alimentation et que le
 ! detrainement est plus fort que le flux de masse, on stope le thermique.
-            endif
-
-            if(l.gt.lmax(ig)) then
+!            endif
+
+!            if(l.gt.lmax(ig)) then
+            if(l.gt.lalim(ig)) then
                detr(ig,l)=0.
                fm(ig,l+1)=0.
                entr(ig,l)=0.
             endif
+            
+            endif
+
          enddo