Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15815)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15816)
@@ -53,30 +53,11 @@
 		segments=numpy.nonzero(numpy.sum(segmentsfront,axis=1)!=numbernodesfront)[0]
 		#Find all nodes for these segments and spc them
-		pos=md.mesh.segments[segments,0:-1]
+		pos=md.mesh.segments[segments,0:-1]-1
 	else:
 		pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-	md.stressbalance.spcvx[pos[:]-1]=0
-	md.stressbalance.spcvy[pos[:]-1]=0
-	md.stressbalance.spcvz[pos[:]-1]=0
+	md.stressbalance.spcvx[pos[:]]=0
+	md.stressbalance.spcvy[pos[:]]=0
+	md.stressbalance.spcvz[pos[:]]=0
 																													   
-	if md.mesh.dimension==2:
-		numbernodesfront=2
-	else:
-		numbernodesfront=4
-	if any(md.mask.icelevelset<=0):
-		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]-1]
-		segmentsfront=1-values
-		numpy.sum(segmentsfront,axis=1)!=numbernodesfront
-		segments=numpy.nonzero(numpy.sum(segmentsfront,axis=1)!=numbernodesfront)[0]
-		print segments
-		#Find all nodes for these segments and spc them
-		pos=md.mesh.segments[segments,0:-1]
-		print pos
-	else:
-		pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-	md.stressbalance.spcvx[pos[:]-1]=0
-	md.stressbalance.spcvy[pos[:]-1]=0
-	md.stressbalance.spcvz[pos[:]-1]=0
-
 	#Dirichlet Values
 	if isinstance(md.inversion.vx_obs,numpy.ndarray) and numpy.size(md.inversion.vx_obs,axis=0)==md.mesh.numberofvertices and isinstance(md.inversion.vy_obs,numpy.ndarray) and numpy.size(md.inversion.vy_obs,axis=0)==md.mesh.numberofvertices:
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15815)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15816)
@@ -61,10 +61,10 @@
 		segments=numpy.nonzero(numpy.sum(segmentsfront,axis=1)!=numbernodesfront)[0]
 		#Find all nodes for these segments and spc them
-		pos=md.mesh.segments[segments,0:-1]
+		pos=md.mesh.segments[segments,0:-1]-1
 	else:
 		pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-	md.stressbalance.spcvx[pos[:]-1]=0
-	md.stressbalance.spcvy[pos[:]-1]=0
-	md.stressbalance.spcvz[pos[:]-1]=0
+	md.stressbalance.spcvx[pos[:]]=0
+	md.stressbalance.spcvy[pos[:]]=0
+	md.stressbalance.spcvz[pos[:]]=0
 
 	#Dirichlet Values