0001 function md=rifttipsrefine(md,resolution,circleradius);
0002
0003
0004
0005
0006
0007 numberofgrids=50;
0008
0009
0010 if strcmpi(md.riftoutline,''),
0011 error('No rifts found in the outlines!');
0012 end
0013
0014
0015 writefile('Rifts.exp.temp',md.riftoutline);
0016
0017
0018 rifts=expread('Rifts.exp.temp',1);
0019 !echo -n "" > Circles.exp
0020 for i=1:length(rifts),
0021 tip1=[rifts(i).x(1) rifts(i).y(1)];
0022 tip2=[rifts(i).x(end) rifts(i).y(end)];
0023
0024 expcreatecircle('Circle1.exp',tip1(1),tip1(2),circleradius,numberofgrids);
0025 expcreatecircle('Circle2.exp',tip2(1),tip2(2),circleradius,numberofgrids);
0026 !cat Circles.exp Circle1.exp Circle2.exp > Circles2.exp
0027 !mv Circles2.exp Circles.exp
0028 !rm -rf Circle1.exp Circle2.exp
0029 end
0030
0031 md=meshexprefine(md,'Circles.exp',resolution);
0032
0033 system('rm -rf Circles.exp');