Changeset 24485

Timestamp:

12/19/19 15:08:08 (5 years ago)

Author:

tsantos

Message:

CHG: minor

File:

• issm/trunk-jpl/src/m/contrib/tsantos/mismip/writeNetCDF.m (modified) (2 diffs)

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

@@ -12 +12 @@
 
         %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);
@@ -26 +31 @@
 
    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);