Index: /issm/trunk-jpl/src/c/solutions/adjointbalancethickness_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/adjointbalancethickness_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/adjointbalancethickness_core.cpp	(revision 11832)
@@ -17,9 +17,7 @@
 	/*parameters: */
 	bool save_results;
-	int  solution_type;
 
 	/*retrieve parameters:*/
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	/*compute thickness */
Index: /issm/trunk-jpl/src/c/solutions/adjointdiagnostic_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/adjointdiagnostic_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/adjointdiagnostic_core.cpp	(revision 11832)
@@ -19,10 +19,8 @@
 	bool save_results;
 	bool conserve_loads   = true;
-	int  solution_type;
 
 	/*retrieve parameters:*/
 	femmodel->parameters->FindParam(&isstokes,FlowequationIsstokesEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	/*Compute velocities*/
Index: /issm/trunk-jpl/src/c/solutions/balancethickness_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/balancethickness_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/balancethickness_core.cpp	(revision 11832)
@@ -16,6 +16,4 @@
 
 	/*parameters: */
-	int  dim;
-	int  solution_type;
 	bool save_results;
 
@@ -24,7 +22,5 @@
 	
 	/*recover parameters: */
-	femmodel->parameters->FindParam(&dim,MeshDimensionEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	_printf_(VerboseSolution(),"call computational core:\n");
Index: /issm/trunk-jpl/src/c/solutions/bedslope_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/bedslope_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/bedslope_core.cpp	(revision 11832)
@@ -14,14 +14,8 @@
 
 	/*parameters: */
-	int  dim;
-	bool isstokes;
-	bool ishutter;
 	bool save_results;
-	int  solution_type;
 
 	/*Recover some parameters: */
-	femmodel->parameters->FindParam(&dim,MeshDimensionEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	_printf_(VerboseSolution(),"%s\n","   computing slope");
Index: /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/hydrology_core.cpp	(revision 11832)
@@ -22,5 +22,4 @@
 	double dt;
 	bool   save_results;
-	int    solution_type;
 	int    output_frequency;
 
@@ -30,5 +29,4 @@
 	femmodel->parameters->FindParam(&dt,TimesteppingTimeStepEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 	femmodel->parameters->FindParam(&output_frequency,SettingsOutputFrequencyEnum);
 	/*first compute slopes: */
Index: /issm/trunk-jpl/src/c/solutions/prognostic_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/prognostic_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/prognostic_core.cpp	(revision 11832)
@@ -15,5 +15,4 @@
 
 	/*parameters: */
-	int solution_type;
 	bool save_results;
 
@@ -23,5 +22,4 @@
 	/*recover parameters: */
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	_printf_(VerboseSolution(),"   call computational core\n");
Index: /issm/trunk-jpl/src/c/solutions/steadystate_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/steadystate_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/steadystate_core.cpp	(revision 11832)
@@ -26,6 +26,4 @@
 	/*parameters: */
 	bool save_results,isenthalpy;
-	int  dim;
-	int  solution_type;
 	int  maxiter;
 	int  numoutputs         = 0;
@@ -33,10 +31,9 @@
 	
 	/* recover parameters:*/
-	femmodel->parameters->FindParam(&dim,MeshDimensionEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 	femmodel->parameters->FindParam(&maxiter,SteadystateMaxiterEnum);
 	femmodel->parameters->FindParam(&numoutputs,SteadystateNumRequestedOutputsEnum);
 	femmodel->parameters->FindParam(&isenthalpy,ThermalIsenthalpyEnum);
+	femmodel->parameters->SetParam(false,SaveResultsEnum);
 	if(numoutputs) femmodel->parameters->FindParam(&requested_outputs,&numoutputs,SteadystateRequestedOutputsEnum);
 
Index: /issm/trunk-jpl/src/c/solutions/surfaceslope_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/surfaceslope_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/surfaceslope_core.cpp	(revision 11832)
@@ -14,14 +14,8 @@
 
 	/*parameters: */
-	int dim;
-	bool isstokes;
-	bool ishutter;
 	bool save_results;
-	int solution_type;
 
 	/*Recover some parameters: */
-	femmodel->parameters->FindParam(&dim,MeshDimensionEnum);
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 
 	_printf_(VerboseSolution(),"%s\n","computing slope...");
Index: /issm/trunk-jpl/src/c/solutions/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutions/transient_core.cpp	(revision 11831)
+++ /issm/trunk-jpl/src/c/solutions/transient_core.cpp	(revision 11832)
@@ -27,5 +27,4 @@
 	bool   save_results,dakota_analysis;
 	bool   time_adapt=false;
-	int    solution_type;
 	int    output_frequency;
 	int    dim,groundingline_migration;
@@ -44,5 +43,4 @@
 	femmodel->parameters->FindParam(&yts,ConstantsYtsEnum);
 	femmodel->parameters->FindParam(&dakota_analysis,QmuIsdakotaEnum);
-	femmodel->parameters->FindParam(&solution_type,SolutionTypeEnum);
 	femmodel->parameters->FindParam(&output_frequency,SettingsOutputFrequencyEnum);
 	femmodel->parameters->FindParam(&time_adapt,TimesteppingTimeAdaptEnum);
Index: /issm/trunk-jpl/src/m/solutions/adjointbalancethickness_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/adjointbalancethickness_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/adjointbalancethickness_core.m	(revision 11832)
@@ -7,5 +7,4 @@
 	%recover parameters common to all solutions
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	%set analysis type to compute velocity:
Index: /issm/trunk-jpl/src/m/solutions/adjointdiagnostic_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/adjointdiagnostic_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/adjointdiagnostic_core.m	(revision 11832)
@@ -7,7 +7,5 @@
 	%recover parameters common to all solutions
 	isstokes=femmodel.parameters.FlowequationIsstokes;
-	dim=femmodel.parameters.MeshDimension;
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 	conserve_loads=true;
 
Index: /issm/trunk-jpl/src/m/solutions/balancethickness_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/balancethickness_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/balancethickness_core.m	(revision 11832)
@@ -6,7 +6,5 @@
 
 	%recover parameters common to all solutions
-	dim=femmodel.parameters.MeshDimension;
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
Index: /issm/trunk-jpl/src/m/solutions/bedslope_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/bedslope_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/bedslope_core.m	(revision 11832)
@@ -7,7 +7,5 @@
 
 	%Recover some parameters:
-	dim=femmodel.parameters.MeshDimension;
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	issmprintf(VerboseSolution,'   computing bed slope');
Index: /issm/trunk-jpl/src/m/solutions/control_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/control_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/control_core.m	(revision 11832)
@@ -17,5 +17,4 @@
 	tol_cm=femmodel.parameters.InversionCostFunctionThreshold;
 	gradient_only=femmodel.parameters.InversionGradientOnly;
-	dim=femmodel.parameters.MeshDimension;
 	isstokes=femmodel.parameters.FlowequationIsstokes;
 	femmodel.parameters.SaveResults=false;
Index: /issm/trunk-jpl/src/m/solutions/diagnostic_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/diagnostic_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/diagnostic_core.m	(revision 11832)
@@ -11,5 +11,4 @@
 
 	%recover parameters common to all solutions
-	dim=femmodel.parameters.MeshDimension;
 	ishutter=femmodel.parameters.FlowequationIshutter;
 	ismacayealpattyn=femmodel.parameters.FlowequationIsmacayealpattyn;
Index: /issm/trunk-jpl/src/m/solutions/enthalpy_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/enthalpy_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/enthalpy_core.m	(revision 11832)
@@ -7,5 +7,4 @@
 	%recover parameters common to all solutions
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	issmprintf(VerboseSolution,'\n%s',['   computing enthalpy']);
Index: /issm/trunk-jpl/src/m/solutions/prognostic_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/prognostic_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/prognostic_core.m	(revision 11832)
@@ -7,5 +7,4 @@
 	%recover parameters common to all solutions
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	%Activate formulation
Index: /issm/trunk-jpl/src/m/solutions/steadystate_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/steadystate_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/steadystate_core.m	(revision 11832)
@@ -7,8 +7,7 @@
 
 	%recover parameters common to all solutions
-	dim=femmodel.parameters.MeshDimension;
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 	isenthalpy=femmodel.parameters.ThermalIsenthalpy;
+	femmodel.parameters.SaveResults=false;
 
 	%Initialize counter
Index: /issm/trunk-jpl/src/m/solutions/surfaceslope_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/surfaceslope_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/surfaceslope_core.m	(revision 11832)
@@ -7,7 +7,5 @@
 
 	%Recover some parameters:
-	dim=femmodel.parameters.MeshDimension;
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	%Call on core computations: 
Index: /issm/trunk-jpl/src/m/solutions/thermal_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/thermal_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/thermal_core.m	(revision 11832)
@@ -8,5 +8,4 @@
 	%recover parameters common to all solutions
 	save_results=femmodel.parameters.SaveResults;
-	solution_type=femmodel.parameters.SolutionType;
 
 	issmprintf(VerboseSolution,'\n%s',['   computing temperature']);
Index: /issm/trunk-jpl/src/m/solutions/transient_core.m
===================================================================
--- /issm/trunk-jpl/src/m/solutions/transient_core.m	(revision 11831)
+++ /issm/trunk-jpl/src/m/solutions/transient_core.m	(revision 11832)
@@ -6,10 +6,8 @@
 
 	%recover parameters common to all solutions
-	dim=femmodel.parameters.MeshDimension;
 	starttime=femmodel.parameters.TimesteppingStartTime;
 	finaltime=femmodel.parameters.TimesteppingFinalTime;
 	dt=femmodel.parameters.TimesteppingTimeStep;
 	yts=femmodel.parameters.ConstantsYts;
-	solution_type=femmodel.parameters.SolutionType;
 	output_frequency=femmodel.parameters.SettingsOutputFrequency;
 	time_adapt=femmodel.parameters.TimesteppingTimeAdapt;