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

Last change on this file since 16279 was 16279, checked in by Eric.Larour, 11 years ago

CHG: issue with order of varargin

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 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.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.