Index: /issm/trunk/src/m/solutions/ice/Loadresults.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/Loadresults.m	(revision 769)
+++ /issm/trunk/src/m/solutions/ice/Loadresults.m	(revision 770)
@@ -8,4 +8,5 @@
 %      md=Loadresults(md,fem,u_g)
 
+if isnan(md.results), md.results=struct(); end
 if strcmpi(md.type,'2d')
 	if fem.ismacayealpattyn,
Index: /issm/trunk/src/m/solutions/ice/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/prognostic.m	(revision 769)
+++ /issm/trunk/src/m/solutions/ice/prognostic.m	(revision 770)
@@ -58,4 +58,5 @@
 
 %project collapsed thickness onto 3d mesh
+if isnan(md.results), md.results=struct(); end
 indx=1:6:m.gridset.gsize; indx=indx(m.tpart);
 md.results.prognostic.step=1;
Index: /issm/trunk/src/m/solutions/ice/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/thermal.m	(revision 769)
+++ /issm/trunk/src/m/solutions/ice/thermal.m	(revision 770)
@@ -45,4 +45,5 @@
 	indx=1:6:gridset.gsize;
 	indx=indx(m_t.tpart);
+	if isnan(md.results), md.results=struct(); end
 	md.results.thermal.step=1;
 	md.results.thermal.time=0;
@@ -92,4 +93,5 @@
 		solution(n).time=soln(n).time/md.yts; %in year
 	end
+	if isnan(md.results), md.results=struct(); end
 	md.results.thermal=solution;
 end
Index: /issm/trunk/src/m/solutions/ice/transient2d.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/transient2d.m	(revision 769)
+++ /issm/trunk/src/m/solutions/ice/transient2d.m	(revision 770)
@@ -114,3 +114,4 @@
 	solution2(i).step=i;
 end
+if isnan(md.results), md.results=struct(); end
 md.results.transient=solution2;
Index: /issm/trunk/src/m/solutions/ice/transient3d.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/transient3d.m	(revision 769)
+++ /issm/trunk/src/m/solutions/ice/transient3d.m	(revision 770)
@@ -164,3 +164,4 @@
 	solution2(i).step=i;
 end
+if isnan(md.results), md.results=struct(); end
 md.results.transient=solution2;
Index: /issm/trunk/src/m/solutions/macayeal/control.m
===================================================================
--- /issm/trunk/src/m/solutions/macayeal/control.m	(revision 769)
+++ /issm/trunk/src/m/solutions/macayeal/control.m	(revision 770)
@@ -239,4 +239,5 @@
 
 	 %Load results onto md:
+	 if isnan(md.results), md.results=struct(); end
 	 md.results.control.step=1;
 	 md.results.control.time=0;
Index: /issm/trunk/src/m/solutions/macayeal/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/macayeal/diagnostic.m	(revision 769)
+++ /issm/trunk/src/m/solutions/macayeal/diagnostic.m	(revision 770)
@@ -468,4 +468,5 @@
 
 %Load results onto md:
+if isnan(md.results), md.results=struct(); end
 md.results.diagnostic.step=1;
 md.results.diagnostic.time=0;