Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17869)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17870)
@@ -42,7 +42,7 @@
 		print "SetMarineIceSheetBC warning: ice front all around the glacier, no dirichlet found. Dirichlet must be added manually."
 
-	md.stressbalance.spcvx=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
-	md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
-	md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
+	md.stressbalance.spcvx=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+	md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+	md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,))
 	md.stressbalance.referential=float('nan')*numpy.ones((md.mesh.numberofvertices,6))
 	md.stressbalance.loadingforce=0*numpy.ones((md.mesh.numberofvertices,3))