Index: /issm/trunk/src/c/solvers/solver_thermal_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/solvers/solver_thermal_nonlinear.cpp	(revision 4521)
+++ /issm/trunk/src/c/solvers/solver_thermal_nonlinear.cpp	(revision 4522)
@@ -83,6 +83,6 @@
 		MatFree(&Kgg);
 	
+		/*!Reduce load from g to f size: */
 		if (verbose) _printf_("   reducing load from g to f set\n");
-		/*!Reduce load from g to f size: */
 		Reduceloadfromgtofx(&pf, pg, fem->Gmn, Kfs, fem->ys, fem->nodesets);
 
@@ -100,15 +100,15 @@
 		MatFree(&Kff);VecFree(&pf);VecFree(&tg);
 
+		//Merge back to g set
 		if (verbose) _printf_("   merging solution from f to g set\n");
-		//Merge back to g set
 		Mergesolutionfromftogx(&tg, tf,fem->Gmn,fem->ys,fem->nodesets);
+
+		//Update inputs using new solution:
+		if (verbose) _printf_("   updating inputs\n");
+		InputUpdateFromSolutionx(fem->elements,fem->nodes, fem->vertices, fem->loads, fem->materials, fem->parameters,tg);
 
 		//Deal with penalty loads
 		if (verbose) _printf_("   penalty constraints\n");
 		PenaltyConstraintsx(&constraints_converged, &num_unstable_constraints, fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters);
-		
-		//Update inputs using new solution:
-		InputUpdateFromVectorx( fem->elements,fem->nodes, fem->vertices, fem->loads, fem->materials, fem->parameters,tg,TemperatureEnum,VertexEnum);
-		InputUpdateFromSolutionx(fem->elements,fem->nodes, fem->vertices, fem->loads, fem->materials, fem->parameters,tg);
 
 		if (!converged){
Index: /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m	(revision 4521)
+++ /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m	(revision 4522)
@@ -39,5 +39,5 @@
 		t_g= Mergesolutionfromftog( t_f, femmodel.Gmn, femmodel.ys, femmodel.nodesets ); 
 
-		displaystring(femmodel.parameters.Verbose,'%s',['   update inputs']);
+		displaystring(femmodel.parameters.Verbose,'%s',['   updating inputs']);
 		femmodel.elements=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,t_g);