Changeset 2479

Timestamp:

10/22/09 10:45:28 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added time stepping in 2d

File:

• issm/trunk/src/m/utils/Mesh/timesteps.m (modified) (1 diff)

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

@@ -15 +15 @@
 end
 
-areas=GetAreas(md.elements,md.x,md.y);
-Dx=averaging(md,sqrt(2*areas),0);
-Du=sqrt(vx.^2+vy.^2);
-maxtime=min(1/2*Dx./Du);
+edgex=max(md.x(md.elements),[],2)-min(md.x(md.elements),[],2);
+edgey=max(md.y(md.elements),[],2)-min(md.y(md.elements),[],2);
+vx=max(vx(md.elements),[],2);
+vy=max(vy(md.elements),[],2);
 
-
+maxtime=1/2*(edgex/vx+edgey/vy);