Index: trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90	(revision 2948)
+++ trunk/LMDZ.MARS/libf/phymars/comsoil_h.F90	(revision 2951)
@@ -15,5 +15,5 @@
        !                 soil_settings.F)
 
-!$OMP THREADPRIVATE(layer,mlayer,inertiedat,volcapa)
+!$OMP THREADPRIVATE(layer,mlayer,inertiedat,inertiesoil,volcapa)
 
   ! variables (FC: built in firstcall in soil.F)
Index: trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F	(revision 2948)
+++ trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F	(revision 2951)
@@ -992,5 +992,5 @@
       write_prof=.false.
       call getin("write_prof",write_prof)
-      IF (write_prof.eq..true.) then
+      IF (write_prof) then
               DO iq = 1, nq
         	call writeprofile(nlayer,q(:,iq),noms(iq),iq)
Index: trunk/LMDZ.MARS/libf/phymars/soil_settings.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/soil_settings.F	(revision 2948)
+++ trunk/LMDZ.MARS/libf/phymars/soil_settings.F	(revision 2951)
@@ -320,10 +320,9 @@
       if (.not.ok) then
         write(*,*)'soil_settings: Field <inertiesoil> not found!'
-	write(*,*)' => Building <inertiesoil> from surface values <inertiedat>'
-	do iloop=1,nsoil
-      do islope=1,nslope
+	write(*,*)' => Building <inertiesoil> from surface values'//
+     & ' <inertiedat>'
+        do islope=1,nslope
 	        inertiesoil(:,:,islope)=inertiedat(:,:)
-       enddo
-	enddo
+        enddo
       else ! <inertiesoil> found
        if (interpol) then ! put values in oldtsoil