Changeset 24029

Timestamp:

06/19/19 05:39:24 (6 years ago)

Author:

hakesson

Message:

CHG: added support for 3d meshes

File:

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

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

@@ -11 +11 @@
 %process data 
 if dimension(md.mesh)==3,
-        error('contourlevelzero error message: routine not supported for 3d meshes, project on a layer');
-end
-
-x=md.mesh.x;
-y=md.mesh.y;
+%       error('contourlevelzero error message: routine not supported for 3d meshes, project on a layer');
+        x = md.mesh.x2d;
+        y = md.mesh.y2d;
+        index=md.mesh.elements2d;
+else
+        x=md.mesh.x;
+        y=md.mesh.y;
+        index=md.mesh.elements;
+end
+
 if isfield(md.mesh,'z'),
         z=md.mesh.z;
@@ -21 +26 @@
         z=zeros(md.mesh.numberofvertices,1);
 end
-index=md.mesh.elements;
 
 if isempty(mask), error('mask provided is empty'); end
-if length(mask)~=md.mesh.numberofvertices, error('mask provided should be specified at the vertices of the mesh'); end
+if dimension(md.mesh)==3,
+        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
+end
 
 %initialization of some variables