Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.py	(revision 17875)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceSheetBC.py	(revision 17876)
@@ -15,7 +15,7 @@
 	#node on Dirichlet
 	pos=numpy.nonzero(md.mesh.vertexonboundary)
-	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.spcvx[pos]=0
 	md.stressbalance.spcvy[pos]=0
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 17875)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 17876)
@@ -33,7 +33,7 @@
 #	pos=find(md.mesh.vertexonboundary & ~nodeonicefront);
 	pos=numpy.nonzero(numpy.logical_and(md.mesh.vertexonboundary,numpy.logical_not(nodeonicefront)))[0]
-	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))
@@ -59,7 +59,7 @@
 	else:
 		pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-	md.stressbalance.spcvx[pos[:]]=0
-	md.stressbalance.spcvy[pos[:]]=0
-	md.stressbalance.spcvz[pos[:]]=0
+	md.stressbalance.spcvx[pos]=0
+	md.stressbalance.spcvy[pos]=0
+	md.stressbalance.spcvz[pos]=0
 																													   
 	#Dirichlet Values
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17875)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 17876)
@@ -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,))
-	md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,))
-	md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+	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))
@@ -68,7 +68,7 @@
 	else:
 		pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-	md.stressbalance.spcvx[pos[:]]=0
-	md.stressbalance.spcvy[pos[:]]=0
-	md.stressbalance.spcvz[pos[:]]=0
+	md.stressbalance.spcvx[pos]=0
+	md.stressbalance.spcvy[pos]=0
+	md.stressbalance.spcvz[pos]=0
 
 	#Dirichlet Values