Last change
on this file since 9714 was 9714, checked in by Mathieu Morlighem, 14 years ago |
moved some stuff to mesh
|
File size:
409 bytes
|
Line | |
---|
1 | function flag=outflow(md);
|
---|
2 | %OUTFLOW - flag nodes on outflux boundary
|
---|
3 | %
|
---|
4 | % Usage:
|
---|
5 | % flag=outflow(md);
|
---|
6 |
|
---|
7 | A=md.mesh.segments(:,1);
|
---|
8 | B=md.mesh.segments(:,2);
|
---|
9 | Nx=-(md.y(A)-md.y(B));
|
---|
10 | Ny= md.x(A)-md.x(B);
|
---|
11 | Vx=(md.initialization.vx(A)+md.initialization.vx(B))/2;
|
---|
12 | Vy=(md.initialization.vy(A)+md.initialization.vy(B))/2;
|
---|
13 |
|
---|
14 | %dot product
|
---|
15 | VdotN=Vx.*Nx+Vy.*Ny;
|
---|
16 |
|
---|
17 | flag=zeros(md.numberofnodes,1);
|
---|
18 | flag(A(find(VdotN>0)))=1;
|
---|
Note:
See
TracBrowser
for help on using the repository browser.