Changeset 26221

Timestamp:

04/29/21 11:52:01 (4 years ago)

Author:

Eric.Larour

Message:

CHG: extended to a 3D plot type.

File:

• issm/trunk-jpl/src/m/exp/contourlevelzero.m (modified) (4 diffs)

issm/trunk-jpl/src/m/exp/contourlevelzero.m

@@ -1 +1 @@
 function contours=contourlevelzero(md,mask,level)
-%CONTOURLEVELZERO - figure out the zero level (or offset thereof, specified by 
-%the level value) of a vectorial mask, and vectorialize it into an exp or shp 
-%compatible structure.
+%CONTOURLEVELZERO - figure out the zero level (or offset thereof, specified by the level value) 
+%                   of a vectorial mask, and vectorialize it into an exp or shp compatible 
+%structure.
 %
 %   Usage:
@@ -9 +9 @@
 %   See also: PLOT_CONTOUR
 
-if isempty(mask)
-        error('mask provided is empty');
-end
-
 %process data 
 if dimension(md.mesh)==3,
+%       error('contourlevelzero error message: routine not supported for 3d meshes, project on a layer');
         x = md.mesh.x2d;
         y = md.mesh.y2d;
@@ -24 +21 @@
 end
 
-if isfield(md.mesh,'z'),
+if isprop(md.mesh,'z'),
         z=md.mesh.z;
 else
@@ -30 +27 @@
 end
 
+if isempty(mask), error('mask provided is empty'); end
 if dimension(md.mesh)==3,
-        if length(mask)~=md.mesh.numberofvertices2d, error('Since this is a 3d model, you need to project the mask on a layer (e.g. project2d(md,mask,0))'); end
+        if length(mask)~=md.mesh.numberofvertices2d, error('mask provided should be specified at the vertices of the mesh'); end
 else
         if length(mask)~=md.mesh.numberofvertices, error('mask provided should be specified at the vertices of the mesh'); end