Changeset 27648

Timestamp:

03/18/23 13:04:05 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: keep 2d segments in 3d mesh for hydro models

Location:

issm/trunk-jpl/src/m/classes

Files:

• mesh3dprisms.m (modified) (4 diffs)
• model.m (modified) (1 diff)

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

@@ -35 +35 @@
                 y2d                         = NaN;
                 elements2d                  = NaN;
+                segments2d                  = NaN;
                 numberofvertices2d          = 0;
                 numberofelements2d          = 0;
@@ -136 +137 @@
                         fielddisplay(self,'numberofvertices2d','number of vertices');
                         fielddisplay(self,'elements2d','vertex indices of the mesh elements');
+                        fielddisplay(self,'segments2d','edges on 2d domain boundary (vertex1 vertex2 element)');
                         fielddisplay(self,'x2d','vertices x coordinate [m]');
                         fielddisplay(self,'y2d','vertices y coordinate [m]');
@@ -190 +192 @@
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','average_vertex_connectivity','format','Integer');
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','elements2d','format','DoubleMat','mattype',3);
+                        WriteData(fid,prefix,'object',self,'class','mesh','fieldname','segments2d','format','DoubleMat','mattype',3);
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofvertices2d','format','Integer');
                         WriteData(fid,prefix,'object',self,'class','mesh','fieldname','numberofelements2d','format','Integer');
@@ -235 +238 @@
                         writejs1Darray(fid,[modelname '.mesh.y2d'],self.y2d);
                         writejs2Darray(fid,[modelname '.mesh.elements2d'],self.elements2d);
+                        writejs2Darray(fid,[modelname '.mesh.segments2d'],self.segments2d);
                         writejsdouble(fid,[modelname '.mesh.numberofvertices2d'],self.numberofvertices2d);
                         writejsdouble(fid,[modelname '.mesh.numberofelements2d'],self.numberofelements2d);

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

@@ -1112 +1112 @@
                         md.mesh.extractedvertices           = mesh2d.extractedvertices;
                         md.mesh.extractedelements           = mesh2d.extractedelements;
+
+                        md.mesh.segments2d                  = mesh2d.segments;
 
                         x3d=[];