Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 2026)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 2027)
@@ -15,5 +15,5 @@
 %recover parameters common to all solutions
 debug=m_dhu.parameters.debug;
-dim=m_dhu.parameters.dim;
+dim=m_dh.parameters.dim;
 ishutter=m_dhu.parameters.ishutter;
 ismacayealpattyn=m_dh.parameters.ismacayealpattyn;
Index: /issm/trunk/src/m/solutions/cielo/gradjcompute_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/gradjcompute_core.m	(revision 2026)
+++ /issm/trunk/src/m/solutions/cielo/gradjcompute_core.m	(revision 2027)
@@ -42,7 +42,7 @@
 end
 
-if analysis_type==SteadystateAnalysisEnum(),
+if m.parameters.control_steady,
 	%compute initial velocity from diagnostic_core (horiz+vertical)
-	displaystring(debug,'%s','          compute initial velocity...');
+	displaystring(debug,'%s','          compute 3d initial velocity...');
 	results_diag=diagnostic_core(models,inputs);
 	u_g=results_diag.u_g;
Index: /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m	(revision 2026)
+++ /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m	(revision 2027)
@@ -17,5 +17,5 @@
 
 %Run diagnostic with updated parameters.
-if (analysis_type==DiagnosticAnalysisEnum()), 
+if ~m.parameters.control_steady;
 	%do a simple diagnostic, with the current temperature profile, do not look for steady state.
 	u_g=diagnostic_core_nonlinear(m,inputs,analysis_type,sub_analysis_type);
Index: /issm/trunk/src/m/solutions/cielo/steadystate.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/steadystate.m	(revision 2026)
+++ /issm/trunk/src/m/solutions/cielo/steadystate.m	(revision 2027)
@@ -50,4 +50,7 @@
 	if ~models.dh.parameters.qmu_analysis,
 		if md.control_analysis,
+			%change control_steady to 1
+			models.dh.parameters.control_steady=1;
+			models.ds.parameters.control_steady=1;
 			%launch core of control solution.
 			results=control_core(models,inputs);