Index: /issm/trunk-jpl/src/m/parameterization/setflowequation.m
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setflowequation.m	(revision 24057)
+++ /issm/trunk-jpl/src/m/parameterization/setflowequation.m	(revision 24058)
@@ -116,5 +116,7 @@
 
 %Now take care of the coupling between SSA and HO
-md.stressbalance.vertex_pairing=[];
+if strcmpi(coupling_method,'penalties'),
+	md.stressbalance.vertex_pairing=[];
+end
 nodeonSSAHO=zeros(md.mesh.numberofvertices,1);
 nodeonHOFS=zeros(md.mesh.numberofvertices,1);
Index: /issm/trunk-jpl/src/m/parameterization/setflowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 24057)
+++ /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 24058)
@@ -107,5 +107,6 @@
 
 	#Now take care of the coupling between SSA and HO
-	md.stressbalance.vertex_pairing=np.array([])
+	if not coupling_method in ['penalties']:
+            md.stressbalance.vertex_pairing=np.array([])
 	nodeonSSAHO=np.zeros(md.mesh.numberofvertices,bool)
 	nodeonHOFS=np.zeros(md.mesh.numberofvertices,bool)