Index: /issm/trunk/src/m/classes/public/thicknessevolution.m
===================================================================
--- /issm/trunk/src/m/classes/public/thicknessevolution.m	(revision 1096)
+++ /issm/trunk/src/m/classes/public/thicknessevolution.m	(revision 1097)
@@ -4,5 +4,5 @@
 %   This routine compute the new thickness of a model after a time step
 %   according to the following formula:
-%   dh/dt=div(Hu)
+%   dh/dt=-div(Hu)
 %
 %   Usage:
@@ -40,3 +40,3 @@
 %compute dhdt=div(Hu)
 Hvx=H.*vx; Hvy=H.*vy;
-dhdt=sum(Hvx(index).*alpha,2)+sum(Hvy(index).*beta,2);
+dhdt=-sum(Hvx(index).*alpha,2)-sum(Hvy(index).*beta,2);