Index: trunk/LMDZ.MARS/libf/phymars/dimradmars.h
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dimradmars.h	(revision 329)
+++ trunk/LMDZ.MARS/libf/phymars/dimradmars.h	(revision 330)
@@ -49,5 +49,5 @@
 
 !  Warning : the "nir" thermal IR bands are not ordered by wavelength:
-!      iir=1 : central 15um CO2 bands \    
+!      iir=1 : central 15um CO2 bands     
 !      iir=2 : CO2 band wings    [long1co2-long2co2] MINUS central band
 !      iir=3 : 9 um band [long1ir - long1co2]
Index: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F	(revision 329)
+++ trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F	(revision 330)
@@ -1,3 +1,3 @@
-         if (callrichsl .eq. .false.) then
+         if (.not.callrichsl) then
        
              DO ig=1,ngrid
Index: trunk/LMDZ.MARS/libf/phymars/physiq.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 329)
+++ trunk/LMDZ.MARS/libf/phymars/physiq.F	(revision 330)
@@ -89,8 +89,8 @@
 c    pt(ngrid,nlayer)      Temperature (K)
 c    pq(ngrid,nlayer,nq)   Advected fields
-c    pudyn(ngrid,nlayer)    \ 
-c    pvdyn(ngrid,nlayer)     \ Dynamical temporal derivative for the
-c    ptdyn(ngrid,nlayer)     / corresponding variables
-c    pqdyn(ngrid,nlayer,nq) /
+c    pudyn(ngrid,nlayer)    | 
+c    pvdyn(ngrid,nlayer)    | Dynamical temporal derivative for the
+c    ptdyn(ngrid,nlayer)    | corresponding variables
+c    pqdyn(ngrid,nlayer,nq) |
 c    pw(ngrid,?)           vertical velocity
 c
@@ -98,9 +98,9 @@
 c   -------
 c
-c    pdu(ngrid,nlayermx)        \
-c    pdv(ngrid,nlayermx)         \  Temporal derivative of the corresponding
-c    pdt(ngrid,nlayermx)         /  variables due to physical processes.
-c    pdq(ngrid,nlayermx,nqmx)   /
-c    pdpsrf(ngrid)             /
+c    pdu(ngrid,nlayermx)       |
+c    pdv(ngrid,nlayermx)       |  Temporal derivative of the corresponding
+c    pdt(ngrid,nlayermx)       |  variables due to physical processes.
+c    pdq(ngrid,nlayermx,nqmx)  |
+c    pdpsrf(ngrid)             |
 c    tracerdyn                 call tracer in dynamical part of GCM ?
 
@@ -1773,5 +1773,5 @@
          if (calltherm .and. (z_out .gt. 0.)) then
          call surflayer_interpol(ngrid,nlayer,z0,g,zzlay,zu,zv,wmax_th
-     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar)
+     &              ,tsurf,zh,z_out,Teta_out,u_out,ustar,tstar,L_mo)
 
          zu2(:)=sqrt(zu(:,1)*zu(:,1)+zv(:,1)*zv(:,1))
Index: trunk/LMDZ.MARS/libf/phymars/thermcell_dqupdown.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/thermcell_dqupdown.F90	(revision 329)
+++ trunk/LMDZ.MARS/libf/phymars/thermcell_dqupdown.F90	(revision 330)
@@ -212,5 +212,7 @@
                  if ((q(ig,k) .gt. 0.) .and. ( q(ig,k) .gt. 1.e-15 )) then
                    !! only print if it is the thermal scheme which makes qa<0
-                   print*,'qa<0 created by thermals!!!',charvar,ig,k,fm0(ig,k),qa(ig,k-1),entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
+                   print*,'qa<0 created by thermals!!!',charvar,ig,k, &
+     &                     fm0(ig,k),qa(ig,k-1),     &
+     &                     entr0(ig,k),q(ig,k),fm0(ig,k+1),detr0(ig,k)
                    print*,'---------> Cancelling qa'
                  endif