Index: /issm/trunk/src/m/solutions/cielo/loadcontrolfinalsol.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/loadcontrolfinalsol.m	(revision 914)
+++ /issm/trunk/src/m/solutions/cielo/loadcontrolfinalsol.m	(revision 915)
@@ -1,5 +1,5 @@
 function md=loadfinalcontrolsol(md,solution);
 
-if isnan(md.results), md.results=struct(); end
+if ~isstruct(md.results), md.results=struct(); end
 md.results.control.vx=solution.vx*md.yts;
 md.results.control.vy=solution.vy*md.yts;
Index: /issm/trunk/src/m/solutions/cielo/loadresults.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/loadresults.m	(revision 914)
+++ /issm/trunk/src/m/solutions/cielo/loadresults.m	(revision 915)
@@ -14,5 +14,5 @@
 isstokes=m_ds.parameters.isstokes;
 
-if isnan(md.results), md.results=struct(); end 
+if ~isstruct(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 914)
+++ /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 915)
@@ -28,5 +28,5 @@
 
 	displaystring(md.debug,'\n%s',['load results...']);
-	if isnan(md.results), md.results=struct(); end
+	if ~isstruct(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 914)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 915)
@@ -29,5 +29,5 @@
 
 	%plug onto the model
-	if isnan(md.results), md.results=struct(); end
+	if ~isstruct(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 914)
+++ /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 915)
@@ -109,4 +109,4 @@
 	results(i).thickness=solution(i).h_g;
 end
-if isnan(md.results), md.results=struct(); end
+if ~isstruct(md.results), md.results=struct(); end
 md.results.transient=results;
Index: /issm/trunk/src/m/solutions/cielo/transient3d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 914)
+++ /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 915)
@@ -149,4 +149,4 @@
 	results(i).melting=solution(i).m_g;
 end
-if isnan(md.results), md.results=struct(); end
+if ~isstruct(md.results), md.results=struct(); end
 md.results.transient=results;