Index: /issm/trunk/src/m/solutions/cielo/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 971)
@@ -10,4 +10,6 @@
 
 	%Build all models requested for diagnostic simulation
+	models.analysis_type='diagnostic'; %needed for processresults
+	
 	displaystring(md.debug,'%s',['reading diagnostic horiz model data']);
 	md.analysis_type='diagnostic'; md.sub_analysis_type='horiz'; models.dh=CreateFemModel(md);
@@ -39,5 +41,5 @@
 		%process results
 		if ~isstruct(md.results), md.results=struct(); end
-		md.results.diagnostic=processresults(results,models,'diagnostic');
+		md.results.diagnostic=processresults(results,models);
 	else
 		%launch dakota driver for diagnostic core solution
Index: /issm/trunk/src/m/solutions/cielo/processresults.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/processresults.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/processresults.m	(revision 971)
@@ -1,3 +1,3 @@
-function newresults=processresults(results,models,analysis_type);
+function newresults=processresults(results,models)
 %PROCESSRESULTS - process results from core solution
 %
@@ -10,4 +10,7 @@
 %initialize output
 newresults=struct();
+
+%get analysis_type:
+analysis_type=models.analysis_type;
 
 %recover models first
Index: /issm/trunk/src/m/solutions/cielo/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 971)
@@ -9,4 +9,6 @@
 
 	%Build all models requested for diagnostic simulation
+	models.analysis_type='prognostic'; %needed for processresults
+	
 	displaystring(md.debug,'%s',['reading prognostic model data']);
 	md.analysis_type='prognostic'; models.p=CreateFemModel(md);
@@ -29,5 +31,5 @@
 	displaystring(md.debug,'\n%s',['load results...']);
 	if ~isstruct(md.results), md.results=struct(); end
-	md.results.prognostic=processresults(results,models,'prognostic');
+	md.results.prognostic=processresults(results,models);
 
 	%stop timing
Index: /issm/trunk/src/m/solutions/cielo/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 971)
@@ -7,4 +7,6 @@
 	%timing
 	t1=clock;
+	
+	models.analysis_type='thermal'; %needed for processresults
 
 	%Build all models requested for diagnostic simulation
@@ -30,5 +32,5 @@
 	%plug onto the model
 	if ~isstruct(md.results), md.results=struct(); end
-	md.results.thermal=processresults(results,models,'thermal');
+	md.results.thermal=processresults(results,models);
 
 	%stop timing
Index: /issm/trunk/src/m/solutions/cielo/transient2d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 971)
@@ -8,4 +8,6 @@
 
 %Build all models requested
+models.analysis_type='transient'; %needed for processresults
+
 displaystring(md.debug,'%s',['reading diagnostic horiz model data']);
 md.analysis_type='diagnostic'; md.sub_analysis_type='horiz'; models.dh=CreateFemModel(md);
Index: /issm/trunk/src/m/solutions/cielo/transient3d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 970)
+++ /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 971)
@@ -8,4 +8,6 @@
 
 %Build all models related to diagnostic
+models.analysis_type='transient'; %needed for processresults
+
 displaystring(md.debug,'%s',['reading diagnostic horiz model data']);
 md.analysis_type='diagnostic'; md.sub_analysis_type='horiz'; models.dh=CreateFemModel(md);
@@ -137,3 +139,3 @@
 %process results
 if ~isstruct(md.results), md.results=struct(); end
-md.results.transient=processresults(results,models,'transient');
+md.results.transient=processresults(results,models);