Index: /issm/trunk-jpl/src/m/model/mesh/triangle.py
===================================================================
--- /issm/trunk-jpl/src/m/model/mesh/triangle.py	(revision 12115)
+++ /issm/trunk-jpl/src/m/model/mesh/triangle.py	(revision 12116)
@@ -36,4 +36,5 @@
 	[md.mesh.elements,md.mesh.x,md.mesh.y,md.mesh.segments,md.mesh.segmentmarkers]=tm.TriMesh(domainname,riftname,area)
 
+
 	#Fill in rest of fields:
 	md.mesh.numberofelements = size(md.mesh.elements)
@@ -41,5 +42,5 @@
 	md.mesh.z = zeros(md.mesh.numberofvertices)
 	md.mesh.vertexonboundary = zeros(md.mesh.numberofvertices)
-	md.mesh.vertexonboundary[md.mesh.segments[:,0:2].astype(int)] = 1.
+	md.mesh.vertexonboundary[md.mesh.segments[:,0:2].astype(int)-1] = 1.
 	md.mesh.vertexonbed = ones(md.mesh.numberofvertices)
 	md.mesh.vertexonsurface = ones(md.mesh.numberofvertices)
@@ -50,5 +51,5 @@
 	[md.mesh.vertexconnectivity]= nc.NodeConnectivity(md.mesh.elements, md.mesh.numberofvertices)
 	[md.mesh.elementconnectivity] = ec.ElementConnectivity(md.mesh.elements, md.mesh.vertexconnectivity)
-
+	
 	#type of model
 	md.mesh.dimension = 2.