Index: /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 80)
+++ /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 81)
@@ -63,10 +63,12 @@
 		if (debug) _printf_("   Updating inputs\n");
 
+		/*Set input parameters: */
+		if(old_ug)ParameterInputsAddFromVec(inputs,old_ug,"old_velocity");
+		if(ug)ParameterInputsAddFromVec(inputs,ug,"velocity");
+
 		//save pointer to old velocity
 		VecFree(&old_ug);old_ug=ug;
 		VecFree(&old_uf);old_uf=uf;
 
-		/*Set input parameters: */
-		if(ug)ParameterInputsAddFromVec(inputs,ug,"velocity");
 
 		/*Update parameters: */
@@ -111,4 +113,5 @@
 		if (debug) _printf_("   penalty constraints\n");
 		ParameterInputsAddFromVec(inputs,ug,"velocity");
+		ParameterInputsAddFromVec(inputs,old_ug,"old_velocity");
 		
 		PenaltyConstraintsx(&constraints_converged, &num_unstable_constraints, fem->elements,fem->nodes,loads,fem->materials,inputs,analysis_type);