Index: ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15797)
+++ ../trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15798)
@@ -55,7 +55,7 @@
 	else:
 		numbernodesfront=4
 	segmentsfront=numpy.nonzeros(md.mask.icelevelset(md.mesh.segments[:,0:-1]))
-	segments=numpy.sum(segmentsfront,axis=1)!=numbernodesfront).astype(int) 
+	segments=(numpy.sum(segmentsfront,axis=1)!=numbernodesfront).astype(int) 
 	#Find all nodes for these segments and spc them
 	pos=md.mesh.segments[segments,0:-1]
 	md.stressbalance.spcvx[pos[:]]=0