Index: trunk/LMDZ.COMMON/libf/evolution/stopping_crit.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/stopping_crit.F90	(revision 4110)
+++ trunk/LMDZ.COMMON/libf/evolution/stopping_crit.F90	(revision 4152)
@@ -450,10 +450,10 @@
         S_atm_2_h2o = S_atm_2_h2o + delta_h2o_ads(i)*cell_area(i)
     else
-        S_h2o_2_atm = S_h2o_2_atm + delta_h2o_ads(i)*cell_area(i)
+        S_h2o_2_atm = S_h2o_2_atm - delta_h2o_ads(i)*cell_area(i)
     end if
     if (delta_icetable(i) > 0._dp) then
         S_atm_2_h2o = S_atm_2_h2o + delta_icetable(i)*cell_area(i)
     else
-        S_h2o_2_atm = S_h2o_2_atm + delta_icetable(i)*cell_area(i)
+        S_h2o_2_atm = S_h2o_2_atm - delta_icetable(i)*cell_area(i)
     end if
     do islope = 1,nslope