Index: /issm/trunk-jpl/src/m/parameterization/setflowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 17872)
+++ /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 17873)
@@ -89,5 +89,5 @@
 		                              numpy.logical_not(numpy.isnan(md.stressbalance.spcvy)).astype(int)+ \
 		                              numpy.logical_not(numpy.isnan(md.stressbalance.spcvz)).astype(int)==3, \
-		                              numpy.logical_and(nodeonHO,nodeonFS).reshape(-1,1)).astype(int)    #find all the nodes on the boundary of the domain without icefront
+		                              numpy.logical_and(nodeonHO,nodeonFS)).astype(int)    #find all the nodes on the boundary of the domain without icefront
 #		fullspcelems=double(sum(fullspcnodes(md.mesh.elements),2)==6);         %find all the nodes on the boundary of the domain without icefront
 		fullspcelems=(numpy.sum(fullspcnodes[md.mesh.elements-1],axis=1)==6).astype(int)    #find all the nodes on the boundary of the domain without icefront