Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 434)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 435)
@@ -25,10 +25,7 @@
 	gridinsideicefront=ContourToMesh(md.elements,md.x,md.y,expread(icefrontfile,1),'node',2);
 	gridonicefront=double(md.gridonboundary & gridinsideicefront);
-
 else
 	%Guess where the ice front is
-	pos=find(md.gridonboundary & md.gridoniceshelf);
-	md.gridondirichlet_diag=zeros(md.numberofgrids,1);
-	md.gridondirichlet_diag(pos)=1;
+	gridonicefront=double(md.gridonboundary & md.gridoniceshelf);
 end
 pos=find(md.gridonboundary & ~gridonicefront);