Changeset 16395

Timestamp:

10/11/13 14:13:03 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added vertexonsurface for 2d flowline models

Location:

issm/trunk-jpl/src/m

Files:

• classes/mesh2dvertical.m (modified) (4 diffs)
• contrib/dassflow/importgmsh.m (modified) (1 diff)
• mesh/bamgflowband.m (modified) (1 diff)

issm/trunk-jpl/src/m/classes/mesh2dvertical.m

@@ -19 +19 @@
                 vertexonboundary            = NaN
                 vertexonbed                 = NaN
+                vertexonsurface             = NaN
 
                 edges                       = NaN
@@ -67 +68 @@
                         md = checkfield(md,'mesh.numberofvertices','>',0);
                         md = checkfield(md,'mesh.vertexonbed','size',[md.mesh.numberofvertices 1],'values',[0 1]);
+                        md = checkfield(md,'mesh.vertexonsurface','size',[md.mesh.numberofvertices 1],'values',[0 1]);
                         md = checkfield(md,'mesh.average_vertex_connectivity','>=',9,'message','''mesh.average_vertex_connectivity'' should be at least 9 in 2d');
 
@@ -89 +91 @@
                         fielddisplay(obj,'vertexonboundary','vertices on the boundary of the domain flag list');
                         fielddisplay(obj,'vertexonbed','vertices on the bed of the domain flag list');
+                        fielddisplay(obj,'vertexonsurface','vertices on the surface of the domain flag list');
                         fielddisplay(obj,'segments','edges on domain boundary (vertex1 vertex2 element)');
                         fielddisplay(obj,'segmentmarkers','number associated to each segment');
@@ -109 +112 @@
                         WriteData(fid,'object',obj,'class','mesh','fieldname','numberofvertices','format','Integer');
                         WriteData(fid,'object',obj,'class','mesh','fieldname','vertexonbed','format','BooleanMat','mattype',1);
+                        WriteData(fid,'object',obj,'class','mesh','fieldname','vertexonsurface','format','BooleanMat','mattype',1);
                         WriteData(fid,'object',obj,'class','mesh','fieldname','average_vertex_connectivity','format','Integer');
                 end % }}}

issm/trunk-jpl/src/m/contrib/dassflow/importgmsh.m

@@ -87 +87 @@
 md.mesh.vertexonbed=zeros(md.mesh.numberofvertices,1);
 md.mesh.vertexonbed(find(vertexflags(md.mesh,1)))=1;
+md.mesh.vertexonsurface=zeros(md.mesh.numberofvertices,1);
+md.mesh.vertexonsurface(find(vertexflags(md.mesh,3)))=1;

issm/trunk-jpl/src/m/mesh/bamgflowband.m

@@ -33 +33 @@
 md.mesh.vertexonbed=zeros(md.mesh.numberofvertices,1);
 md.mesh.vertexonbed(find(vertexflags(md.mesh,1)))=1;
+md.mesh.vertexonsurface=zeros(md.mesh.numberofvertices,1);
+md.mesh.vertexonsurface(find(vertexflags(md.mesh,3)))=1;