Changeset 1097

Timestamp:

06/25/09 11:17:34 (16 years ago)

Author:

Mathieu Morlighem

Message:

dh/dt= - div(Hu)

File:

• issm/trunk/src/m/classes/public/thicknessevolution.m (modified) (2 diffs)

issm/trunk/src/m/classes/public/thicknessevolution.m

@@ -4 +4 @@
 %   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 +40 @@
 %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);