Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_thermal_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_thermal_nonlinear.cpp	(revision 15456)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_thermal_nonlinear.cpp	(revision 15457)
@@ -47,9 +47,9 @@
 	for(;;){
 
+		delete tf_old; tf_old=tf;
 		femmodel->SystemMatricesx(&Kff, &Kfs, &pf,&df, &melting_offset);
 		CreateNodalConstraintsx(&ys,femmodel->nodes,configuration_type);
-		Reduceloadx(pf, Kfs, ys); delete Kfs; delete tf;
+		Reduceloadx(pf, Kfs, ys); delete Kfs;
 		Solverx(&tf, Kff, pf,tf_old, df, femmodel->parameters);
-		delete tf_old; tf_old=tf->Duplicate();
 		delete Kff;delete pf;delete tg; delete df;
 		Mergesolutionfromftogx(&tg, tf,ys,femmodel->nodes,femmodel->parameters); delete ys;