Changeset 24514

Timestamp:

01/14/20 12:42:55 (5 years ago)

Author:

tsantos

Message:

CHG: verification

File:

• issm/trunk-jpl/src/m/contrib/tsantos/mismip/gl_position.m (modified) (1 diff)

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

@@ -2 +2 @@
 
                 %initialization of some variables
-                data                                    = md.results.TransientSolution(step).MaskGroundediceLevelset;
-                if(isfield(md.results.TransientSolution,'MeshElements'))
-                        index                                   = md.results.TransientSolution(step).MeshElements;
-                        x                                               = md.results.TransientSolution(step).MeshX;
-                        y                                               = md.results.TransientSolution(step).MeshY;
+                if isfield(md.results,'TransientSolution')
+                        data            = md.results.TransientSolution(step).MaskGroundediceLevelset;
+                        if isfield(md.results.TransientSolution,'MeshElements')
+                                index   = md.results.TransientSolution(step).MeshElements;
+                                x               = md.results.TransientSolution(step).MeshX;
+                                y               = md.results.TransientSolution(step).MeshY;
+                        else
+                                index   = md.mesh.elements;
+                                x               = md.mesh.x;
+                                y               = md.mesh.y;
+                        end
                 else
-                        index                                   = md.mesh.elements;
-                        x                                               = md.mesh.x;
-                        y                                               = md.mesh.y;
+                        data    = md.mask.groundedice_levelset;
+                        index   = md.mesh.elements;
+                        x               = md.mesh.x;
+                        y               = md.mesh.y;
                 end
                 numberofelements        = size(index,1);