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
|
Line | |
---|
1 | function flag=outflo(md);
|
---|
2 |
|
---|
3 | A=md.segments(:,1);
|
---|
4 | B=md.segments(:,2);
|
---|
5 | Nx=-(md.y(A)-md.y(B));
|
---|
6 | Ny= md.x(A)-md.x(B);
|
---|
7 | Vx=(md.vx(A)+md.vy(B))/2;
|
---|
8 | Vy=(md.vy(A)+md.vy(B))/2;
|
---|
9 |
|
---|
10 | %dot product
|
---|
11 | VdotN=Vx.*Nx+Vy.*Ny;
|
---|
12 |
|
---|
13 | flag=zeros(md.numberofgrids,1);
|
---|
14 | flag(A(find(VdotN>0)))=1;
|
---|
Note:
See
TracBrowser
for help on using the repository browser.