Changeset 2323

Timestamp:

09/25/09 08:58:08 (15 years ago)

Author:

Mathieu Morlighem

Message:

fixed Stokes

File:

• issm/trunk/src/m/classes/public/setelementstype.m (modified) (1 diff)

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

@@ -61 +61 @@
 gridonstokes=ones(md.numberofgrids,1);
 gridonstokes(md.elements(nonstokesflag,:))=0;      %non stokes grids
-pos=find(sum(md.spcvelocity(:,1:3),2)==3);         %find all the grids on the boundary of the domain without icefront
+pos=find(sum(md.spcvelocity(:,1:2),2)==2);         %find all the grids on the boundary of the domain without icefront
 gridonstokes(pos)=0;                               %we have to constrain all the boundary of the domain without icefront