source: issm/oecreview/Archive/24307-24683/ISSM-24484-24485.diff@ 24684

Index: ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m
===================================================================
--- ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m     (revision 24484)
+++ ../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m     (revision 24485)
@@ -11,8 +11,13 @@
    time=[];
 
        %Inserting time 0. md0 must be last experiment (e.g., Ice1r for Ice1ra)
-       x                                                       = md0.results.TransientSolution(end).MeshX;
-       y                                                       = md0.results.TransientSolution(end).MeshY;
+       if(isfield(md0.results.TransientSolution,'MeshElements'))
+      x     = md0.results.TransientSolution(end).MeshX;
+      y     = md0.results.TransientSolution(end).MeshY;
+   else
+      x     = md0.mesh.x;
+      y     = md0.mesh.y;
+   end 
        time(1)                                 = 0;
        [xgl_step ygl_step]     = gl_position(md0,length(md0.results.TransientSolution),0);
        xGL{1}                                  = xgl_step;
@@ -25,8 +30,13 @@
        vBaseGL{1}                              = griddata(x,y,md0.results.TransientSolution(end).Vy,xgl_step,ygl_step);
 
    for i=2:length(step),
-               x = md.results.TransientSolution(step(i)).MeshX;
-               y = md.results.TransientSolution(step(i)).MeshY;
+               if(isfield(md.results.TransientSolution,'MeshElements'))
+                       x     = md.results.TransientSolution(step(i)).MeshX;
+                       y     = md.results.TransientSolution(step(i)).MeshY;
+               else
+                       x     = md.mesh.x;
+                       y     = md.mesh.y;
+               end     
                time(i)=md.results.TransientSolution(step(i)).time;     
                [xgl_step ygl_step]=gl_position(md,step(i),0);
       xGL{i}=xgl_step;