meshprocessrifts

PURPOSE ^

MESHPROCESSRIFTS - process mesh when rifts are present

SYNOPSIS ^

function md=meshprocessrifts(md)

DESCRIPTION ^

MESHPROCESSRIFTS - process mesh when rifts are present

   split rifts inside mesh (rifts are defined by presence of
   segments inside the domain outline)

   Usage:
      md=meshprocessrifts(md)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SUBFUNCTIONS ^

SOURCE CODE ^

0001 function md=meshprocessrifts(md)
0002 %MESHPROCESSRIFTS - process mesh when rifts are present
0003 %
0004 %   split rifts inside mesh (rifts are defined by presence of
0005 %   segments inside the domain outline)
0006 %
0007 %   Usage:
0008 %      md=meshprocessrifts(md)
0009 
0010 %some checks on list of arguments
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 %Fill in rest of fields:
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)');

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003