Index: /issm/trunk/src/m/solutions/adjoint_core.m
===================================================================
--- /issm/trunk/src/m/solutions/adjoint_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/adjoint_core.m	(revision 5317)
@@ -9,4 +9,5 @@
 	isstokes=femmodel.parameters.IsStokes;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 	conserve_loads=true;
@@ -34,6 +35,6 @@
 	femmodel=solver_adjoint_linear(femmodel);
 
-	displaystring(verbose,'\n%s',['      saving results...']);
-	if(solution_type==AdjointSolutionEnum)
+	if(solution_type==AdjointSolutionEnum & ~control_analysis)
+		displaystring(verbose,'\n%s',['      saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,AdjointxEnum);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,AdjointyEnum);
Index: /issm/trunk/src/m/solutions/balancedthickness_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/balancedthickness_core.m	(revision 5317)
@@ -8,4 +8,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -16,5 +17,5 @@
 	femmodel=solver_linear(femmodel);
 
-	if solution_type==BalancedthicknessSolutionEnum,
+	if (solution_type==BalancedthicknessSolutionEnum & ~control_analysis),
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
Index: /issm/trunk/src/m/solutions/balancedvelocities_core.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/balancedvelocities_core.m	(revision 5317)
@@ -8,4 +8,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -16,5 +17,5 @@
 	femmodel=solver_linear(femmodel);
 	
-	if solution_type==BalancedvelocitiesSolutionEnum,
+	if (solution_type==BalancedvelocitiesSolutionEnum & ~control_analysis),
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VelEnum);
Index: /issm/trunk/src/m/solutions/bedslope_core.m
===================================================================
--- /issm/trunk/src/m/solutions/bedslope_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/bedslope_core.m	(revision 5317)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -19,5 +20,5 @@
 	femmodel=solver_linear(femmodel);
 	
-	if solution_type==BedSlopeSolutionEnum,
+	if (solution_type==BedSlopeSolutionEnum & ~control_analysis),
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,BedSlopeXEnum);
Index: /issm/trunk/src/m/solutions/control_core.m
===================================================================
--- /issm/trunk/src/m/solutions/control_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/control_core.m	(revision 5317)
@@ -72,4 +72,5 @@
 	%generate output
 	displaystring(verbose,'\n%s',['      preparing final velocity solution...']);
+	femmodel.parameters.ControlAnalysis=0;
 	eval(['femmodel=' solutioncore '(femmodel);']);
 
Index: /issm/trunk/src/m/solutions/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/diagnostic_core.m	(revision 5317)
@@ -18,4 +18,5 @@
 	if(isstokes) stokesreconditioning=femmodel.parameters.StokesReconditioning; end
 	qmu_analysis=femmodel.parameters.QmuAnalysis;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -71,5 +72,5 @@
 	end
 
-	if solution_type==DiagnosticSolutionEnum,
+	if (solution_type==DiagnosticSolutionEnum & ~control_analysis)
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
Index: /issm/trunk/src/m/solutions/prognostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/prognostic_core.m	(revision 5317)
@@ -7,4 +7,5 @@
 	%recover parameters common to all solutions
 	verbose=femmodel.parameters.Verbose;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -15,5 +16,5 @@
 	femmodel=solver_linear(femmodel);
 	
-	if solution_type==PrognosticSolutionEnum,
+	if (solution_type==PrognosticSolutionEnum & ~control_analysis)
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,ThicknessEnum);
Index: /issm/trunk/src/m/solutions/steadystate_core.m
===================================================================
--- /issm/trunk/src/m/solutions/steadystate_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/steadystate_core.m	(revision 5317)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -38,5 +39,5 @@
 	end
 
-	if solution_type==SteadystateSolutionEnum,
+	if (solution_type==SteadystateSolutionEnum & ~control_analysis),
 		displaystring(verbose,'\n%s',['saving results...']);
 		femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum);
Index: /issm/trunk/src/m/solutions/surfaceslope_core.m
===================================================================
--- /issm/trunk/src/m/solutions/surfaceslope_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/surfaceslope_core.m	(revision 5317)
@@ -9,4 +9,5 @@
 	verbose=femmodel.parameters.Verbose;
 	dim=femmodel.parameters.Dim;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -19,5 +20,5 @@
 	femmodel=solver_linear(femmodel);
 	
-	if solution_type==SurfaceSlopeSolutionEnum,
+	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);
Index: /issm/trunk/src/m/solutions/thermal_core.m
===================================================================
--- /issm/trunk/src/m/solutions/thermal_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/thermal_core.m	(revision 5317)
@@ -10,4 +10,5 @@
 	ndt=femmodel.parameters.Ndt;
 	dt=femmodel.parameters.Dt;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 
@@ -28,5 +29,5 @@
 		femmodel=thermal_core_step(femmodel,i,time); 
 
-		if solution_type==ThermalSolutionEnum,
+		if (solution_type==ThermalSolutionEnum & ~control_analysis),
 			displaystring(verbose,'\n%s',['saving results...']);
 			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,TemperatureEnum,i,time);
Index: /issm/trunk/src/m/solutions/transient2d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient2d_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/transient2d_core.m	(revision 5317)
@@ -10,4 +10,5 @@
 	dt=femmodel.parameters.Dt;
 	yts=femmodel.parameters.Yts;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 	output_frequency=femmodel.parameters.OutputFrequency;
@@ -35,5 +36,5 @@
 		[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);
 
-		if (solution_type==Transient2DSolutionEnum & mod(step,output_frequency)==0),
+		if (solution_type==Transient2DSolutionEnum & ~control_analysis & mod(step,output_frequency)==0),
 			displaystring(verbose,'\n%s',['saving results...']);
 			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);
Index: /issm/trunk/src/m/solutions/transient3d_core.m
===================================================================
--- /issm/trunk/src/m/solutions/transient3d_core.m	(revision 5316)
+++ /issm/trunk/src/m/solutions/transient3d_core.m	(revision 5317)
@@ -10,4 +10,5 @@
 	dt=femmodel.parameters.Dt;
 	yts=femmodel.parameters.Yts;
+	control_analysis=femmodel.parameters.ControlAnalysis;
 	solution_type=femmodel.parameters.SolutionType;
 	output_frequency=femmodel.parameters.OutputFrequency;
@@ -41,5 +42,5 @@
 		[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);
 
-		if (solution_type==Transient3DSolutionEnum & mod(step,output_frequency)==0),
+		if (solution_type==Transient3DSolutionEnum & ~control_analysis & mod(step,output_frequency)==0),
 			displaystring(verbose,'\n%s',['saving results...']);
 			femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,VxEnum,step,time);