Changeset 12116

Timestamp:

04/24/12 21:42:08 (13 years ago)

Author:

Eric.Larour

Message:

Indexing bug

File:

• issm/trunk-jpl/src/m/model/mesh/triangle.py (modified) (3 diffs)

issm/trunk-jpl/src/m/model/mesh/triangle.py

@@ -36 +36 @@
         [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 +42 @@
         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 +51 @@
         [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.