Index: /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m	(revision 4557)
+++ /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m	(revision 4558)
@@ -19,4 +19,7 @@
 	ug=GetSolutionFromInputs(femmodel.elements, femmodel.nodes, femmodel.vertices, loads, femmodel.materials, femmodel.parameters);
 	uf=Reducevectorgtof( ug, femmodel.nodesets);
+
+	%Update the solution to make sure that vx and vxold are similar
+	femmodel.elements=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,loads,femmodel.materials,femmodel.parameters,ug);
 
 	while(~converged),
@@ -78,7 +81,4 @@
 		
 	end
-	
-	%Now that the solution has converged, update once again the solution to make sure that vx and vxold are similar
-	femmodel.elements=InputUpdateFromSolution(femmodel.elements,femmodel.nodes,femmodel.vertices,loads,femmodel.materials,femmodel.parameters,ug);
 
 	%extrude if we are in 3D: