Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 13464)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 13465)
@@ -25,4 +25,5 @@
 		if not os.path.exists(icefrontfile):
 			raise IOError("SetIceShelfBC error message: ice front file '%s' not found." % icefrontfile)
+		#nodeonicefront= numpy.logical_and( md.mesh.vertexonboundary, nodeinsideicefront.astype(float)) )
 		nodeinsideicefront=ContourToMesh(md.mesh.elements,md.mesh.x.reshape(-1,1),md.mesh.y.reshape(-1,1),icefrontfile,'node',2)
 		nodeonicefront= numpy.bitwise_and( map(int,md.mesh.vertexonboundary), map(int,nodeinsideicefront[0].ravel()) )
Index: /issm/trunk-jpl/test/Par/SquareShelf.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelf.py	(revision 13464)
+++ /issm/trunk-jpl/test/Par/SquareShelf.py	(revision 13465)
@@ -81,7 +81,5 @@
 #Boundary conditions:
 # #md=SetIceShelfBC(md)
-print '...calling seticeshelfbc:'
 md=SetIceShelfBC(md,'../Exp/SquareFront.exp')
-print '...called seticeshelfbc:'
 
 #Change name so that no test have the same name