Changeset 770

Timestamp:

06/04/09 13:46:13 (16 years ago)

Author:

Mathieu Morlighem

Message:

initialize md.results if needed

Location:

issm/trunk/src/m/solutions

Files:

• ice/Loadresults.m (modified) (1 diff)
• ice/prognostic.m (modified) (1 diff)
• ice/thermal.m (modified) (2 diffs)
• ice/transient2d.m (modified) (1 diff)
• ice/transient3d.m (modified) (1 diff)
• macayeal/control.m (modified) (1 diff)
• macayeal/diagnostic.m (modified) (1 diff)

issm/trunk/src/m/solutions/ice/Loadresults.m

@@ -8 +8 @@
 %      md=Loadresults(md,fem,u_g)
 
+if isnan(md.results), md.results=struct(); end
 if strcmpi(md.type,'2d')
         if fem.ismacayealpattyn,

issm/trunk/src/m/solutions/ice/prognostic.m

@@ -58 +58 @@
 
 %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;

issm/trunk/src/m/solutions/ice/thermal.m

@@ -45 +45 @@
         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 +93 @@
                 solution(n).time=soln(n).time/md.yts; %in year
         end
+        if isnan(md.results), md.results=struct(); end
         md.results.thermal=solution;
 end

issm/trunk/src/m/solutions/ice/transient2d.m

@@ -114 +114 @@
         solution2(i).step=i;
 end
+if isnan(md.results), md.results=struct(); end
 md.results.transient=solution2;

issm/trunk/src/m/solutions/ice/transient3d.m

@@ -164 +164 @@
         solution2(i).step=i;
 end
+if isnan(md.results), md.results=struct(); end
 md.results.transient=solution2;

issm/trunk/src/m/solutions/macayeal/control.m

@@ -239 +239 @@
 
          %Load results onto md:
+         if isnan(md.results), md.results=struct(); end
          md.results.control.step=1;
          md.results.control.time=0;

issm/trunk/src/m/solutions/macayeal/diagnostic.m

@@ -468 +468 @@
 
 %Load results onto md:
+if isnan(md.results), md.results=struct(); end
 md.results.diagnostic.step=1;
 md.results.diagnostic.time=0;