Index: /issm/trunk/src/m/solutions/cielo/control.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/control.m	(revision 1038)
+++ /issm/trunk/src/m/solutions/cielo/control.m	(revision 1039)
@@ -12,6 +12,18 @@
 	models.analysis_type='control'; %needed for processresults
 
-	displaystring(md.debug,'%s',['reading diagnostic horiz model data']);
-	md.analysis_type='control'; md.sub_analysis_type=''; models.dh=CreateFemModel(md); 
+	displaystring(md.debug,'%s',['reading control horiz model data']);
+	md.analysis_type='control'; md.sub_analysis_type='horiz'; models.dh=CreateFemModel(md);
+
+	displaystring(md.debug,'\n%s',['reading control vert model data']);
+	md.analysis_type='control'; md.sub_analysis_type='vert'; models.dv=CreateFemModel(md);
+
+	displaystring(md.debug,'\n%s',['reading control stokes model data']);
+	md.analysis_type='control'; md.sub_analysis_type='stokes'; models.ds=CreateFemModel(md);
+
+	displaystring(md.debug,'\n%s',['reading control hutter model data']);
+	md.analysis_type='control'; md.sub_analysis_type='hutter'; models.dhu=CreateFemModel(md);
+
+	displaystring(md.debug,'\n%s',['reading surface and bed slope computation model data']);
+	md.analysis_type='slope_compute'; md.sub_analysis_type=''; models.sl=CreateFemModel(md);
 
 	% figure out number of dof: just for information purposes.
@@ -31,5 +43,5 @@
 		md.results.control=processresults(models,results);
 	else
-		%launch dakota driver for diagnostic core solution
+		%launch dakota driver for control core solution
 		Qmu(models,inputs,models.dh.parameters);
 	end