Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp	(revision 15924)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/CreateConstraintsStressbalance.cpp	(revision 15925)
@@ -74,5 +74,15 @@
 		else if(isL1L2) finiteelement = P1Enum;
 		else if(isHO)   iomodel->Constant(&finiteelement,FlowequationFeHOEnum);
-		else if(isFS)   iomodel->Constant(&finiteelement,FlowequationFeFSEnum);
+		else if(isFS){  iomodel->Constant(&finiteelement,FlowequationFeFSEnum);
+			/*Deduce velocity interpolation from finite element*/
+			switch(finiteelement){
+				case P1P1Enum          : finiteelement = P1Enum;       break;
+				case P1P1GLSEnum       : finiteelement = P1Enum;       break;
+				case MINIcondensedEnum : finiteelement = P1bubbleEnum; break;
+				case MINIEnum          : finiteelement = P1bubbleEnum; break;
+				case TaylorHoodEnum    : finiteelement = P2Enum;       break;
+				default: _error_("finite element "<<finiteelement<<" not supported");
+			}
+		}
 
 		IoModelToConstraintsx(constraints,iomodel,StressbalanceSpcvxEnum,StressbalanceAnalysisEnum,finiteelement,1);