Index: /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp	(revision 4338)
+++ /issm/trunk/src/c/solvers/solver_diagnostic_nonlinear.cpp	(revision 4339)
@@ -125,4 +125,7 @@
 	}
 	
+	//Now that the solution has converged, update once again the solution to make sure that vx and vxold are similar (for next step in transient or steadystate)
+	InputUpdateFromSolutionx( femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters,ug);
+
 	/*extrude if we are in 3D: */
 	if (dim==3){
Index: /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m	(revision 4338)
+++ /issm/trunk/src/m/solvers/solver_diagnostic_nonlinear.m	(revision 4339)
@@ -79,4 +79,7 @@
 	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: 
 	if dim==3,