source: issm/trunk/src/m/model/outflow.m@ 6959

Last change on this file since 6959 was 6959, checked in by Mathieu Morlighem, 14 years ago

Added outflow.m which flag all node on boundary that are on the outflow
File size: 258 bytes
RevLine 
[6959]1function flag=outflo(md);
2
3A=md.segments(:,1);
4B=md.segments(:,2);
5Nx=-(md.y(A)-md.y(B));
6Ny= md.x(A)-md.x(B);
7Vx=(md.vx(A)+md.vy(B))/2;
8Vy=(md.vy(A)+md.vy(B))/2;
9
10%dot product
11VdotN=Vx.*Nx+Vy.*Ny;
12
13flag=zeros(md.numberofgrids,1);
14flag(A(find(VdotN>0)))=1;
Note: See TracBrowser for help on using the repository browser.