Index: /issm/trunk/src/m/solutions/cielo/loadresults.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/loadresults.m	(revision 768)
+++ /issm/trunk/src/m/solutions/cielo/loadresults.m	(revision 769)
@@ -14,4 +14,5 @@
 isstokes=m_ds.parameters.isstokes;
 
+if isnan(md.results), md.results=struct(); end 
 if dim==2,
 	results.step=1;
Index: /issm/trunk/src/m/solutions/cielo/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 768)
+++ /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 769)
@@ -28,4 +28,5 @@
 
 	displaystring(md.debug,'\n%s',['load results...']);
+	if isnan(md.results), md.results=struct(); end
 	md.results.prognostic.step=1;
 	md.results.prognostic.time=0;
Index: /issm/trunk/src/m/solutions/cielo/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 768)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 769)
@@ -29,4 +29,5 @@
 
 	%plug onto the model
+	if isnan(md.results), md.results=struct(); end
 	md.results.thermal=results;
 
Index: /issm/trunk/src/m/solutions/cielo/transient2d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 768)
+++ /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 769)
@@ -109,3 +109,4 @@
 	results(i).thickness=solution(i).h_g;
 end
+if isnan(md.results), md.results=struct(); end
 md.results.transient=results;