Changeset 1341

Timestamp:

07/17/09 07:47:44 (16 years ago)

Author:

Mathieu Morlighem

Message:

Added routine to compute mesh quality

Location:

issm/trunk/src/m/classes/public/mesh

Files:

• meshquality.m (added)
• meshyams.m (modified) (1 diff)

issm/trunk/src/m/classes/public/mesh/meshyams.m

@@ -111 +111 @@
 md.elementonbed=ones(md.numberofelements,1);
 md.elementonsurface=ones(md.numberofelements,1);
+if ~strcmp(groundeddomain,'N/A'),
+        gridground=ContourToMesh(md.elements,md.x,md.y,expread(groundeddomain,1),'node',2);
+        md.gridonwater=ones(md.numberofgrids,1);
+        md.gridonwater(find(gridground))=0;
+else
+        md.gridonwater=zeros(md.numberofgrids,1);
+end
+md.vx_obs=InterpFromGrid(x_m,y_m,vx,md.x,md.y,0);
+md.vy_obs=InterpFromGrid(x_m,y_m,vy,md.x,md.y,0);
+md.vel_obs=sqrt(md.vx_obs.^2+md.vy_obs.^2);