Index: trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 658)
+++ trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90	(revision 659)
@@ -381,5 +381,6 @@
       a1=1.4716 ; b1=0.0005698 ; ae=0.03683 ; be = 0.57421
       ad = 0.00048088  ; bd = -0.6697
-      fdfu = -1.3
+!      fdfu = -1.3
+      fdfu=-0.8
       a1inv=a1
       b1inv=b1
@@ -404,6 +405,9 @@
 
 ! Best config for norad 222 levels:
-
-       omega=0.06
+! and more specifically to C.large :
+
+!       omega=0.06
+!       omega=0.
+       omega=-0.03
 !       omega=0.
        a1=1.
@@ -413,8 +417,9 @@
        be=1.1*be
        ad = 0.0004
-       adalim=0.
-
-       b1inv=0.00025
-!       b1inv=b1
+!       ad=0.0003
+!       adalim=0.
+
+!       b1inv=0.00025
+       b1inv=b1
 
 !       b1inv = 0.0003
@@ -478,5 +483,5 @@
 !        do l=2,4
          do ig=1,ngridmx
-           if (ztv(ig,l)>(ztv(ig,l+1)+0.) .and. ztv(ig,1)>=ztv(ig,l) .and. (alim_star(ig,l-1) .ne. 0.)) then ! .and. (zlev(ig,l+1) .lt. 1000.)) then
+           if (ztv(ig,l)>(ztv(ig,l+1 + 0.)) .and. ztv(ig,1)>=ztv(ig,l) .and. (alim_star(ig,l-1).ne. 0.)) then ! .and. (zlev(ig,l+1) .lt. 1000.)) then
                alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.)  &
 !     &                       *log(1.+zlev(ig,l+1))
@@ -1279,5 +1284,5 @@
               if(zlay(ig,l) .le. zmax(ig)) then
               fm_down(ig,l) =fm(ig,l)* &
-     &      max(fdfu,-3*max(0.,(zlay(ig,l)/zmax(ig)))-0.9)
+     &      max(fdfu,-4*max(0.,(zlay(ig,l)/zmax(ig)))-0.6)
               endif