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

Last change on this file since 15749 was 15749, checked in by Mathieu Morlighem, 12 years ago

CHG: removing elements and vertices on water, which is now included in icelevelset

File size: 1.4 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 x=md.mesh.x;
15 y=md.mesh.y;
16 index=md.mesh.elements;
17else
18 x=varargin{1};
19 y=varargin{2};
20 index=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.dimension=2;
39md.mesh.numberofelements=size(md.mesh.elements,1);
40md.mesh.numberofvertices=length(md.mesh.x);
41md.mesh.numberofedges=size(md.mesh.edges,1);
42md.mesh.z=zeros(md.mesh.numberofvertices,1);
43md.mesh.vertexonbed=ones(md.mesh.numberofvertices,1);
44md.mesh.vertexonsurface=ones(md.mesh.numberofvertices,1);
45md.mesh.elementonbed=ones(md.mesh.numberofelements,1);
46md.mesh.elementonsurface=ones(md.mesh.numberofelements,1);
47md.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.