Changeset 1244

Timestamp:

07/07/09 10:56:02 (15 years ago)

Author:

Mathieu Morlighem

Message:

fixed dirichlet icefront

File:

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

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

@@ -54 +54 @@
 
 %now, make sure all elements on water have grids that are spc'd, otherwise, we'll get a singular problem.
-pos=find(md.elementonwater);
-grids=unique(md.elements(pos,:));
+pos=find(~md.elementonwater);
+numpos=unique(md.elements(pos,:));
+grids=setdiff(1:1:md.numberofgrids,numpos);
 
 md.gridondirichlet_diag(grids)=1;