Index: /issm/trunk/src/c/parallel/OutputResults.cpp
===================================================================
--- /issm/trunk/src/c/parallel/OutputResults.cpp	(revision 692)
+++ /issm/trunk/src/c/parallel/OutputResults.cpp	(revision 693)
@@ -36,4 +36,3 @@
 		pfclose(fid,filename);
 	}
-
 }
Index: /issm/trunk/src/c/parallel/ProcessResults.cpp
===================================================================
--- /issm/trunk/src/c/parallel/ProcessResults.cpp	(revision 692)
+++ /issm/trunk/src/c/parallel/ProcessResults.cpp	(revision 693)
@@ -286,3 +286,2 @@
 	*presults=newresults;
 }
-
Index: /issm/trunk/src/c/parallel/thermal.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal.cpp	(revision 692)
+++ /issm/trunk/src/c/parallel/thermal.cpp	(revision 693)
@@ -86,4 +86,7 @@
 	inputs->Add("dt",dt);
 
+	_printf_("initialize results:\n");
+	results=new DataSet(ResultsEnum());
+
 	//erase velocities and pressure embedded in parameters
 	param=(Param*)femmodels[0].parameters->FindParamObject("u_g");
@@ -115,8 +118,7 @@
 		#endif
 	}
-
 	
 	_printf_("process results:\n");
-	ProcessResults(&results,&femmodels[0],DiagnosticAnalysisEnum());
+	ProcessResults(&results,&femmodels[0],ThermalAnalysisEnum());
 	
 	_printf_("write results to disk:\n");