Changeset 9733 for issm/trunk/src/m/model/bamg.m

Timestamp:

09/09/11 14:58:09 (14 years ago)

Author:

seroussi

Message:

keep building mesh

File:

• issm/trunk/src/m/model/bamg.m (modified) (4 diffs)

issm/trunk/src/m/model/bamg.m

@@ -271 +271 @@
         else
                 bamg_mesh.Vertices=[md.x md.y ones(md.mesh.numberofvertices,1)];
-                bamg_mesh.Triangles=[md.elements ones(md.mesh.numberofelements,1)];
+                bamg_mesh.Triangles=[md.mesh.elements ones(md.mesh.numberofelements,1)];
         end
 
@@ -319 +319 @@
 md.x=bamgmesh_out.Vertices(:,1);
 md.y=bamgmesh_out.Vertices(:,2);
-md.elements=bamgmesh_out.Triangles(:,1:3);
-md.edges=bamgmesh_out.IssmEdges;
+md.mesh.elements=bamgmesh_out.Triangles(:,1:3);
+md.mesh.edges=bamgmesh_out.IssmEdges;
 md.mesh.segments=bamgmesh_out.IssmSegments(:,1:3);
 md.mesh.segmentmarkers=bamgmesh_out.IssmSegments(:,4);
@@ -326 +326 @@
 %Fill in rest of fields:
 md.mesh.dimension=2;
-md.mesh.numberofelements=size(md.elements,1);
+md.mesh.numberofelements=size(md.mesh.elements,1);
 md.mesh.numberofvertices=length(md.x);
-md.mesh.numberofedges=size(md.edges,1);
+md.mesh.numberofedges=size(md.mesh.edges,1);
 md.z=zeros(md.mesh.numberofvertices,1);
 md.mesh.vertexonbed=ones(md.mesh.numberofvertices,1);
@@ -340 +340 @@
 
 %Check for orphan
-if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.elements(:)))))
+if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))))
         error('Output mesh has orphans. Decrease MaxCornerAngle to prevent outside points (ex: 0.01)');
 end