Index: trunk/LMDZ.COMMON/libf/evolution/surf_ice.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/surf_ice.F90	(revision 4174)
+++ trunk/LMDZ.COMMON/libf/evolution/surf_ice.F90	(revision 4180)
@@ -16,5 +16,5 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, qp, di, k4, minieps
+use numerics, only: dp, qp, di, k4, eps, tol
 
 ! DECLARATION
@@ -30,5 +30,5 @@
 integer(di),                              parameter :: WEIGHT_RESERVOIR = 5_di ! Weight by locally available H2O ice mass
 integer(di),                              parameter :: WEIGHT_DIRECTION = 6_di ! Weight by directional capacity (up/down)
-integer(di),                              private   :: weight_method = WEIGHT_DIRECTION ! Default method for balancing
+integer(di),                              private   :: weight_method = WEIGHT_COMBINED ! Default method for balancing
 real(dp),                                 parameter :: rho_co2ice = 1650._dp ! Density of CO2 ice [kg.m-3]
 real(dp),                                 parameter :: rho_h2oice = 920._dp  ! Density of H2O ice [kg.m-3]
@@ -420,6 +420,4 @@
 ! ---------------
 integer(di), parameter :: max_iter = 50_di ! Maximum number of iterations for the balancing procedure
-real(dp),    parameter :: eps = 1.e-12_dp ! Small number to prevent division by zero in weights normalization
-real(dp),    parameter :: tiny_corr = minieps ! Minimum correction to consider that progress is made in the balancing procedure
 integer(di)                              :: i, islope, iter
 integer(di)                              :: method_used
@@ -556,5 +554,6 @@
     Delta = Delta - Delta_used
 
-    if (abs(Delta_used) < tiny_corr) exit
+    ! Test if the balancing procedure progresses
+    if (abs(Delta_used) < tol) exit
 end do