Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15805)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15806)
@@ -38,4 +38,8 @@
 	md.stressbalance.loadingforce=0*numpy.ones((md.mesh.numberofvertices,3))
 
+	#Icefront position
+	pos=numpy.nonzero(nodeonicefront)[0]
+	md.mask.icelevelset[pos]=0
+
 	#First find segments that are not completely on the front
 	if md.mesh.dimension==2:
@@ -64,8 +68,4 @@
 	else:
 		print "      boundary conditions for stressbalance model: spc set as zero"
-
-	#Icefront position
-	pos=numpy.nonzero(nodeonicefront)[0]
-	md.mask.icelevelset[pos]=0
 
 	#Create zeros basalforcings and surfaceforcings