Index: trunk/LMDZ.COMMON/libf/evolution/sorption.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/sorption.F90	(revision 4174)
+++ trunk/LMDZ.COMMON/libf/evolution/sorption.F90	(revision 4180)
@@ -18,5 +18,5 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, qp, di, k4, minieps
+use numerics, only: dp, qp, di, k4, eps
 
 ! DECLARATION
@@ -230,5 +230,5 @@
         deltam_reg_slope(ig,islope) = 0._dp
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 if (iloop == 1) then
                     deltam_reg_complete(ig,iloop,islope) = (dm_h2o_regolith_slope(ig,iloop,islope) - h2o_ads_reg(ig,iloop,islope))*(layer(iloop))
@@ -347,9 +347,9 @@
     do islope = 1,nslope
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1._dp - theta_h2o_ads(ig,iloop,islope))*alpha*pco2_avg(ig)/ &
                                                          (alpha*pco2_avg(ig) + sqrt(tsoil(ig,iloop,islope))*exp(beta/tsoil(ig,iloop,islope)))
             else
-                if (abs(co2_ads_reg(ig,iloop,islope)) < minieps) then !!! we are at first call
+                if (abs(co2_ads_reg(ig,iloop,islope)) < eps) then !!! we are at first call
                     dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1._dp - theta_h2o_ads(ig,iloop,islope))*alpha*pco2_avg(ig) &
                                                              /(alpha*pco2_avg(ig)+sqrt(tsoil(ig,iloop,islope))*exp(beta/tsoil(ig,iloop,islope)))
@@ -368,5 +368,5 @@
         deltam_reg_slope(ig,islope) = 0._dp
         do iloop = 1,index_breccia
-            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < minieps .and. abs(co2_ice(ig,islope)) < minieps) then
+            if (TI(ig,iloop,islope) < inertie_thresold .and. abs(h2o_ice(ig,islope)) < eps .and. abs(co2_ice(ig,islope)) < eps) then
                 if (iloop == 1) then
                     deltam_reg_complete(ig,iloop,islope) = (dm_co2_regolith_slope(ig,iloop,islope) - co2_ads_reg(ig,iloop,islope))*(layer(iloop))