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
|
Rev | Line | |
---|
[6959] | 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.