Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17102)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17103)
@@ -3303,4 +3303,5 @@
 void StressbalanceAnalysis::InputUpdateFromSolutionFS(IssmDouble* solution,Element* element){/*{{{*/
 
+	bool         results_on_nodes;
 	int          i,dim,meshtype;
 	int*         vdoflist=NULL;
@@ -3310,4 +3311,5 @@
 	element->FindParam(&meshtype,MeshTypeEnum);
 	element->FindParam(&FSreconditioning,StressbalanceFSreconditioningEnum);
+	element->FindParam(&results_on_nodes,SettingsResultsOnNodesEnum);
 	switch(meshtype){
 		case Mesh2DverticalEnum: dim = 2; break;
@@ -3377,14 +3379,18 @@
 
 	/*Add vx and vy as inputs to the tria element: */
-	//element->AddInput(VxEnum,      vx,      element->VelocityInterpolation());
-	//element->AddInput(VyEnum,      vy,      element->VelocityInterpolation());
-	//element->AddInput(VelEnum,     vel,     element->VelocityInterpolation());
-	//element->AddInput(PressureEnum,pressure,element->PressureInterpolation());
-	//if(dim==3) element->AddInput(VzEnum,vz, element->VelocityInterpolation());
-	element->AddInput(VxEnum,vx,P1Enum);
-	element->AddInput(VyEnum,vy,P1Enum);
-	element->AddInput(VelEnum,vel,P1Enum);
-	element->AddInput(PressureEnum,pressure,P1Enum);
-	if(dim==3) element->AddInput(VzEnum,vz,P1Enum);
+	if(results_on_nodes){
+		element->AddInput(VxEnum,      vx,      element->VelocityInterpolation());
+		element->AddInput(VyEnum,      vy,      element->VelocityInterpolation());
+		element->AddInput(VelEnum,     vel,     element->VelocityInterpolation());
+		element->AddInput(PressureEnum,pressure,element->PressureInterpolation());
+		if(dim==3) element->AddInput(VzEnum,vz, element->VelocityInterpolation());
+	}
+	else{
+		element->AddInput(VxEnum,vx,P1Enum);
+		element->AddInput(VyEnum,vy,P1Enum);
+		element->AddInput(VelEnum,vel,P1Enum);
+		element->AddInput(PressureEnum,pressure,P1Enum);
+		if(dim==3) element->AddInput(VzEnum,vz,P1Enum);
+	}
 
 	/*Free ressources:*/