0001 function md=meshprocessrifts(md)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 if ((nargin~=1) | (nargout~=1)),
0012 meshprocessriftsusage();
0013 error('meshprocessrifts error message');
0014 end
0015
0016 [md.elements,md.x,md.y,md.segments,md.segmentmarkers,md.rifts]=TriMeshProcessRifts(md.elements,md.x,md.y,md.segments,md.segmentmarkers);
0017
0018
0019 md.numberofelements=length(md.elements);
0020 md.numberofgrids=length(md.x);
0021 md.z=zeros(md.numberofgrids,1);
0022 md.gridonboundary=zeros(length(md.x),1); md.gridonboundary(md.segments(:,1:2))=1;
0023 md.numrifts=length(md.rifts);
0024 md.elements_type=3*ones(md.numberofelements,1);
0025 md.gridonbed=ones(md.numberofgrids,1);
0026 md.gridonsurface=ones(md.numberofgrids,1);
0027 md.elementonbed=ones(md.numberofelements,1);
0028 md.elementonsurface=ones(md.numberofelements,1);
0029
0030 function meshprocessriftsusage(),
0031 disp('usage: md=meshprocessrifts(md)');