Changeset 26930

Timestamp:

03/11/22 06:01:41 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: handle both 2d and 3d models

File:

• issm/trunk-jpl/src/m/plot/plot_transient_movie.m (modified) (2 diffs)

issm/trunk-jpl/src/m/plot/plot_transient_movie.m

@@ -93 +93 @@
 
                         %Add grounding line
-                        if exist(options,'groundingline')
-                                contours=isoline(md, results(i).MaskOceanLevelset,'output','matrix');
-                                hold on
+                        if exist(options,'groundingline') 
+                                if dimension(md.mesh)==2
+                                        contours=isoline(md, results(i).MaskOceanLevelset,'output','matrix');
+                                else
+                                        ocean = project2d(md, results(i).MaskOceanLevelset, 1);
+                                        contours=isoline(md, ocean,'output','matrix');
+                                end
                                 plot(contours(:,1),contours(:,2),getfieldvalue(options,'groundingline'));
                         end
@@ -101 +105 @@
                         %Add ice front
                         if exist(options,'icefront')
-                                contours=isoline(md, results(i).MaskIceLevelset,'output','matrix');
+                                if dimension(md.mesh)==2
+                                        contours=isoline(md, results(i).MaskIceLevelset,'output','matrix');
+                                else
+                                        ice = project2d(md, results(i).MaskIceLevelset, 1);
+                                        contours=isoline(md, ice,'output','matrix');
+                                end
                                 hold on
                                 plot(contours(:,1),contours(:,2),getfieldvalue(options,'icefront'));