Changeset 2484

Timestamp:

10/22/09 11:06:18 (15 years ago)

Author:

seroussi

Message:

fixed timesteps for 2d

File:

• issm/trunk/src/m/utils/Mesh/timesteps.m (modified) (2 diffs)

issm/trunk/src/m/utils/Mesh/timesteps.m

@@ -2 +2 @@
 %TIMESTEPS - return the maximum timesteps for the model in years
 %
-%   Dt < c Dx/u
+%   Dt < c / ( u/Dx +v/Dy )
 %
 %   Usage:
@@ -21 +21 @@
 vy=max(abs(vy(index)),[],2);
 
-maxtime=min(1/2*(edgex./vx+edgey./vy));
+maxtime=1/2*min(1./(vx./edgex+vy./edgey));