Index: /issm/trunk/src/m/solutions/balancedthickness.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/balancedthickness.m	(revision 4135)
@@ -22,4 +22,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=balancedthickness_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/balancedthickness2.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness2.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/balancedthickness2.m	(revision 4135)
@@ -19,4 +19,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=balancedthickness2_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/balancedvelocities.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedvelocities.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/balancedvelocities.m	(revision 4135)
@@ -22,4 +22,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=balancedvelocities_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/prognostic.m	(revision 4135)
@@ -23,4 +23,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=prognostic_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/prognostic2.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic2.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/prognostic2.m	(revision 4135)
@@ -22,4 +22,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=prognostic2_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/steadystate.m
===================================================================
--- /issm/trunk/src/m/solutions/steadystate.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/steadystate.m	(revision 4135)
@@ -33,4 +33,8 @@
 
 		end
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
+
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/thermal.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/thermal.m	(revision 4135)
@@ -22,4 +22,7 @@
 		displaystring(verbose,'%s',['call computational core']);
 		femmodel=thermal_core(femmodel);
+
+		displaystring(verbose,'%s',['write results']);
+		OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/transient3d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4134)
+++ /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4135)
@@ -49,4 +49,10 @@
 		InputToResult(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,TemperatureEnum,step,time);
 		InputToResult(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,MeltingEnum,step,time);
+
+		%Save temporary results
+		if mod(step,5)==0,
+			displaystring(verbose,'\n%s',['saving temporary results...']);
+			OutputResults(femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices, femmodel.materials, femmodel.parameters);
+		end
 	end