Index: /issm/trunk-jpl/src/c/cores/control_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/control_core.cpp	(revision 18137)
+++ /issm/trunk-jpl/src/c/cores/control_core.cpp	(revision 18138)
@@ -252,7 +252,4 @@
 	}
 
-	/*Needed for output results (FIXME: should be placed 6 lines below)*/
-	ControlInputSetGradientx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,G);
-
 	for(long i=0;i<nsize;i++){
 		if(X[i]>=XU[i]) G[i]=0.;
@@ -260,28 +257,6 @@
 	}
 
-	/*solve forward: (FIXME: not needed actually...)*/
-	switch(solution_type){
-		case SteadystateSolutionEnum:
-			femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum);
-			stressbalance_core(femmodel);	//We need a 3D velocity!! (vz is required for the next thermal run)
-			break;
-		case StressbalanceSolutionEnum:
-			femmodel->SetCurrentConfiguration(StressbalanceAnalysisEnum);
-			solutionsequence_nonlinear(femmodel,conserve_loads); 
-			break;
-		case BalancethicknessSolutionEnum:
-			femmodel->SetCurrentConfiguration(BalancethicknessAnalysisEnum);
-			solutionsequence_linear(femmodel); 
-			break;
-		case BalancethicknessSoftSolutionEnum:
-			/*NOTHING*/
-			break;
-		case Balancethickness2SolutionEnum:
-			femmodel->SetCurrentConfiguration(Balancethickness2AnalysisEnum);
-			solutionsequence_linear(femmodel); 
-			break;
-		default:
-			_error_("Solution " << EnumToStringx(solution_type) << " not implemented yet");
-	}
+	/*Needed for output results*/
+	ControlInputSetGradientx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,G);
 
 	/*Compute misfit for this velocity field.*/