Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15811)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15812)
@@ -48,5 +48,5 @@
 		numbernodesfront=4
 	if any(md.mask.icelevelset<=0):
-		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]]
+		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]-1]
 		segmentsfront=1-values
 		segments=(numpy.sum(segmentsfront,axis=1)!=numbernodesfront).astype(int)
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15811)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15812)
@@ -56,5 +56,5 @@
 		numbernodesfront=4
 	if any(md.mask.icelevelset<=0):
-		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]]
+		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]-1]
 		segmentsfront=1-values
 		segments=(numpy.sum(segmentsfront,axis=1)!=numbernodesfront).astype(int)