Index: /issm/trunk/src/m/solutions/balancedthickness2_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness2_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/balancedthickness2_core.m	(revision 4335)
@@ -8,4 +8,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
@@ -25,6 +26,8 @@
 	%femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	if solution_type==Balancedthickness2SolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	end
 	
 end %end function
Index: /issm/trunk/src/m/solutions/balancedthickness_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 4335)
@@ -8,4 +8,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
@@ -25,6 +26,8 @@
 	femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	if solution_type==BalancethicknessSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	end
 	
 end %end function
Index: /issm/trunk/src/m/solutions/balancedvelocities_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 4335)
@@ -8,4 +8,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
@@ -26,7 +27,9 @@
 	femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
+	if solution_type==BalancedvelocitiesSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
+	end
 
 end %end function
Index: /issm/trunk/src/m/solutions/bedslope_core.m
===================================================================
--- /issm/trunk/src/m/solutions/bedslope_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/bedslope_core.m	(revision 4335)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	displaystring(verbose,'\n%s',['computing bed slope...']);
@@ -25,5 +26,7 @@
 	}
 	
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeXEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeYEnum);
+	if solution_type==BesSlopeSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeXEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeYEnum);
+	end
Index: /issm/trunk/src/m/solutions/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 4335)
@@ -18,4 +18,5 @@
 	if(isstokes) stokesreconditioning=femmodel.parameters.stokesreconditioning; end
 	qmu_analysis=femmodel.parameters.QmuAnalysis;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%for qmu analysis, be sure the velocity input we are starting from  is the one in the parameters: 
@@ -70,8 +71,10 @@
 	end
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
-	if(dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum); end
+	if solution_type==DiagnosticSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
+		if(dim==3) femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum); end
+	end
Index: /issm/trunk/src/m/solutions/prognostic2_core.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic2_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/prognostic2_core.m	(revision 4335)
@@ -7,4 +7,5 @@
 	%recover parameters common to all solutions
 	verbose=femmodel.parameters.Verbose;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
@@ -21,6 +22,8 @@
 %	femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	if solution_type==Prognostic2SolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	end
 	
 end %end function
Index: /issm/trunk/src/m/solutions/prognostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/prognostic_core.m	(revision 4335)
@@ -7,4 +7,5 @@
 	%recover parameters common to all solutions
 	verbose=femmodel.parameters.Verbose;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
@@ -21,6 +22,8 @@
 	femmodel.elements=InputExtrude(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	if solution_type==PrognosticSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
+	end
 
 end %end function
Index: /issm/trunk/src/m/solutions/steadystate_core.m
===================================================================
--- /issm/trunk/src/m/solutions/steadystate_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/steadystate_core.m	(revision 4335)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Initialize counter
@@ -43,11 +44,13 @@
 	end
 
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
+	if solution_type==SteadystateSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum);
 	end
+end
 
 end %end of function
Index: /issm/trunk/src/m/solutions/surfaceslope_core.m
===================================================================
--- /issm/trunk/src/m/solutions/surfaceslope_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/surfaceslope_core.m	(revision 4335)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	solution_type=femmodel.parameters.SolutionType;
 
 	displaystring(verbose,'\n%s',['computing surface slope...']);
@@ -25,7 +26,9 @@
 	}
 	
-	displaystring(verbose,'\n%s',['saving results...']);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeXEnum);
-	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeYEnum);
+	if solution_type==SurfaceSlopeSolutionEnum,
+		displaystring(verbose,'\n%s',['saving results...']);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeXEnum);
+		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceSlopeYEnum);
+	end
 
 end %end function
Index: /issm/trunk/src/m/solutions/thermal_core.m
===================================================================
--- /issm/trunk/src/m/solutions/thermal_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/thermal_core.m	(revision 4335)
@@ -10,4 +10,5 @@
 	ndt=femmodel.parameters.Ndt;
 	dt=femmodel.parameters.Dt;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Compute number of timesteps
@@ -27,7 +28,9 @@
 		femmodel=thermal_core_step(femmodel,i,time); 
 
-		displaystring(verbose,'\n%s',['saving results...']);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum,i,time);
+		if solution_type==ThermalSolutionEnum,
+			displaystring(verbose,'\n%s',['saving results...']);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingRateEnum,i,time);
+		end
 	end
 
Index: /issm/trunk/src/m/solutions/transient2d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient2d_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/transient2d_core.m	(revision 4335)
@@ -9,4 +9,5 @@
 	ndt=femmodel.parameters.Ndt;
 	dt=femmodel.parameters.Dt;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Initialize
@@ -29,12 +30,13 @@
 		[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateGeometry(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 
-		displaystring(verbose,'\n%s',['saving results...']);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
+		if solution_type==Transient2DSolutionEnum,
+			displaystring(verbose,'\n%s',['saving results...']);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
+		end
 	end
-
 end %end of function
Index: /issm/trunk/src/m/solutions/transient3d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4334)
+++ /issm/trunk/src/m/solutions/transient3d_core.m	(revision 4335)
@@ -9,4 +9,5 @@
 	ndt=femmodel.parameters.Ndt;
 	dt=femmodel.parameters.Dt;
+	solution_type=femmodel.parameters.SolutionType;
 
 	%Initialize
@@ -38,14 +39,16 @@
 		[femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters]=UpdateVertexPositions(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 
-		displaystring(verbose,'\n%s',['saving results...']);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,step,time);
-		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingEnum,step,time);
+		if solution_type==Transient3DSolutionEnum,
+			displaystring(verbose,'\n%s',['saving results...']);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VyEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VzEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,PressureEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,SurfaceEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,step,time);
+			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,MeltingEnum,step,time);
+		end
 	end
 end %end of function