Index: /issm/trunk-jpl/src/m/classes/mesh3dprisms.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/mesh3dprisms.m	(revision 27647)
+++ /issm/trunk-jpl/src/m/classes/mesh3dprisms.m	(revision 27648)
@@ -35,4 +35,5 @@
 		y2d                         = NaN;
 		elements2d                  = NaN;
+		segments2d                  = NaN;
 		numberofvertices2d          = 0;
 		numberofelements2d          = 0;
@@ -136,4 +137,5 @@
 			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,4 +192,5 @@
 			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,4 +238,5 @@
 			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);
Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 27647)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 27648)
@@ -1112,4 +1112,6 @@
 			md.mesh.extractedvertices           = mesh2d.extractedvertices;
 			md.mesh.extractedelements           = mesh2d.extractedelements;
+
+			md.mesh.segments2d                  = mesh2d.segments;
 
 			x3d=[];