Index: trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90	(revision 2980)
+++ trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90	(revision 2982)
@@ -106,20 +106,20 @@
    print*, 'Do we compute the ice table with the dynamic method?', icetable_dynamic
 
-  if ((not(soil_pem)).and.((icetable_equilibrium).or.(icetable_dynamic))) then
+  if ((.not.soil_pem).and.((icetable_equilibrium).or.(icetable_dynamic))) then
        print*,'Ice table  must be used when soil_pem = T'
        call abort_physic(modname,"Ice table  must be used when soil_pem = T",1)
   endif
 
-  if ((not(soil_pem)).and.adsorption_pem) then
+  if ((.not.soil_pem).and.adsorption_pem) then
        print*,'Adsorption must be used when soil_pem = T'
        call abort_physic(modname,"Adsorption must be used when soil_pem = T",1)
   endif
   
-  if ((not(soil_pem)).and.(fluxgeo.gt.0.)) then
+  if ((.not.soil_pem).and.(fluxgeo.gt.0.)) then
        print*,'Soil is not activated but Flux Geo > 0.'
        call abort_physic(modname,"Soil is not activated but Flux Geo > 0.",1)
   endif
  
-  if ((not(soil_pem)).and.reg_thprop_dependp) then
+  if ((.not.soil_pem).and.reg_thprop_dependp) then
      print*,'Regolith properties vary with Ps only when soil is set to true'
      call abort_physic(modname,'Regolith properties vary with Ps only when soil is set to true',1)
Index: trunk/LMDZ.COMMON/libf/evolution/orbit_param_criterion_mod.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/orbit_param_criterion_mod.F90	(revision 2980)
+++ trunk/LMDZ.COMMON/libf/evolution/orbit_param_criterion_mod.F90	(revision 2982)
@@ -14,4 +14,10 @@
 
       SUBROUTINE orbit_param_criterion(year_iter_max)
+#ifdef CPP_IOIPSL
+  use IOIPSL, only: getin
+#else
+  ! if not using IOIPSL, we still need to use (a local version of) getin
+  use ioipsl_getincom, only: getin
+#endif
 
       USE temps_mod_evol, ONLY: year_bp_ini, year_PEM, var_obl, var_ex, var_lsp
@@ -109,7 +115,13 @@
 !Constant max change case 
 
-        max_change_obl=0.1
+        max_change_obl=0.5
         max_change_ex=0.1
         max_change_lsp=40.
+
+        CALL getin('max_change_obl', max_change_obl)
+
+        CALL getin('max_change_ex', max_change_ex)
+
+        CALL getin('max_change_lsp', max_change_lsp)
 
 	max_obl=obliquit+max_change_obl
@@ -147,7 +159,7 @@
         endif
 
-        max_obl_iter=999999999999
-        max_ex_iter =999999999999
-        max_lsp_iter=999999999999
+        max_obl_iter=999999
+        max_ex_iter =999999
+        max_lsp_iter=999999
 
 !--------------------------------
Index: trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 2980)
+++ trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 2982)
@@ -921,5 +921,5 @@
 
       print *, "Updating the new Tsurf"
-      bool_sublim=0
+      bool_sublim=.false.
       Tsurfave_before_saved(:,:) = tsurf_ave(:,:)
       DO ig = 1,ngrid
@@ -931,9 +931,9 @@
                   if(initial_co2_ice(ig_loop,islope_loop).lt.0.5 .and. qsurf(ig_loop,igcm_co2,islope_loop).LT. 1E-10) then
                     tsurf_ave(ig,islope)=tsurf_ave(ig_loop,islope_loop)
-                    bool_sublim=1
+                    bool_sublim=.true.
                     exit
                   endif
                 enddo
-                if (bool_sublim.eq.1) then
+                if (bool_sublim.eqv. .true.) then
                   exit
                 endif
@@ -944,9 +944,9 @@
                   if(initial_co2_ice(ig_loop,islope_loop).lt.0.5 .and. qsurf(ig_loop,igcm_co2,islope_loop).LT. 1E-10) then
                     tsurf_ave(ig,islope)=tsurf_ave(ig_loop,islope_loop)
-                    bool_sublim=1
+                    bool_sublim=.true.
                     exit
                   endif
                 enddo
-                if (bool_sublim.eq.1) then
+                if (bool_sublim.eqv. .true.) then
                   exit
                 endif