source: issm/trunk-jpl/src/m/mesh/meshconvert.m@ 17558

Last change on this file since 17558 was 17558, checked in by Mathieu Morlighem, 11 years ago

NEW: removed old mesh classes, use subclasses instead

File size: 1.2 KB
Line 
1function md=meshconvert(md,varargin)
2%CONVERTMESH - convert mesh to bamg mesh
3%
4% Usage:
5% md=meshconvert(md);
6% md=meshconvert(md,index,x,y);
7
8if nargin~=1 & nargin~=4,
9 help meshconvert
10 error('meshconvert error message: bad usage');
11end
12
13if nargin==1,
14 index = md.mesh.elements;
15 x = md.mesh.x;
16 y = md.mesh.y;
17else
18 index = varargin{1};
19 x = varargin{2};
20 y = varargin{3};
21end
22
23%call Bamg
24[bamgmesh_out bamggeom_out]=BamgConvertMesh(index,x,y);
25
26% plug results onto model
27md.private.bamg = struct();
28md.private.bamg.mesh = bamgmesh(bamgmesh_out);
29md.private.bamg.geometry = bamggeom(bamggeom_out);
30md.mesh.x = bamgmesh_out.Vertices(:,1);
31md.mesh.y = bamgmesh_out.Vertices(:,2);
32md.mesh.elements = bamgmesh_out.Triangles(:,1:3);
33md.mesh.edges = bamgmesh_out.IssmEdges;
34md.mesh.segments = bamgmesh_out.IssmSegments(:,1:3);
35md.mesh.segmentmarkers = bamgmesh_out.IssmSegments(:,4);
36
37%Fill in rest of fields:
38md.mesh.numberofelements = size(md.mesh.elements,1);
39md.mesh.numberofvertices = length(md.mesh.x);
40md.mesh.numberofedges = size(md.mesh.edges,1);
41md.mesh.vertexonboundary = zeros(md.mesh.numberofvertices,1); md.mesh.vertexonboundary(md.mesh.segments(:,1:2)) = 1;
Note: See TracBrowser for help on using the repository browser.