Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17527)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17528)
@@ -838,5 +838,5 @@
 	femmodel->parameters->FindParam(&newton,StressbalanceIsnewtonEnum);
 
-	if(isFS){
+	if(isFS && !(isSSA || isHO || isL1L2)){
 		if(VerboseSolution()) _printf0_("   computing velocities\n");
 
@@ -849,5 +849,5 @@
 		 solutionsequence_nonlinear(femmodel,conserve_loads); 
 	}
-	else if(isSSA || isHO || isL1L2){ 
+	else if(!isFS && (isSSA || isHO || isL1L2)){ 
 		if(VerboseSolution()) _printf0_("   computing velocities\n");
 
@@ -865,11 +865,10 @@
 		}
 	}
+	else if ((isSSA || isL1L2 || isHO) && isFS){
+		if(VerboseSolution()) _printf0_("   computing coupling between lower order models and FS\n");
+		solutionsequence_FScoupling_nonlinear(femmodel,conserve_loads);
+	}
 	else{
 		_error_("not supported");
-	}
-
-	if ((isSSA || isL1L2 || isHO) && isFS){
-		if(VerboseSolution()) _printf0_("   computing coupling between lower order models and FS\n");
-		solutionsequence_FScoupling_nonlinear(femmodel,conserve_loads);
 	}