Index: /issm/trunk/src/m/solutions/cielo/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 1848)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 1849)
@@ -26,11 +26,17 @@
 	inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.numberofnodes);
 	inputs=add(inputs,'dt',models.t.parameters.dt*models.t.parameters.yts,'double');
+
+	%compute solution
+	if ~models.t.parameters.qmu_analysis,
+		%launch core of diagnostic solution.
+		results=thermal_core(models,inputs);
 	
-	%call core
-	results=thermal_core(models,inputs);
-
-	%plug onto the model
-	if ~isstruct(md.results), md.results=struct(); end
-	md.results.thermal=processresults(models,results);
+		%process results
+		if ~isstruct(md.results), md.results=struct(); end
+		md.results.thermal=processresults(models,results);
+	else
+		%launch dakota driver for diagnostic core solution
+		Qmu(models,inputs,models.t.parameters);
+	end
 
 	%stop timing