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

../trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m

@@ -11 +11 @@
    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 +30 @@
         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;