Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 7334)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 7335)
@@ -27,5 +27,7 @@
 else
 	%Guess where the ice front is
-	gridonicefront=double(md.gridonboundary & md.gridoniceshelf);
+	gridoniceshelf=zeros(md.numberofgrids,1);
+	gridoniceshelf(md.elements(find(md.elementoniceshelf),:))=1;
+	gridonicefront=double(md.gridonboundary & gridoniceshelf);
 end
 pos=find(md.gridonboundary & ~gridonicefront);