Index: /issm/trunk/src/m/classes/public/BasinConstrain.m
===================================================================
--- /issm/trunk/src/m/classes/public/BasinConstrain.m	(revision 1243)
+++ /issm/trunk/src/m/classes/public/BasinConstrain.m	(revision 1244)
@@ -54,6 +54,7 @@
 
 %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;