Context Navigation

outflow.m

Last change on this file was 20500, checked in by Mathieu Morlighem, 9 years ago
merged trunk-jpl and trunk for revision 20497
File size: 436 bytes

Rev	Line
[13730]	1	function flag=outflow(md)
[9075]	2	%OUTFLOW - flag nodes on outflux boundary
	3	%
	4	% Usage:
	5	% flag=outflow(md);
[6959]	6
[9714]	7	A=md.mesh.segments(:,1);
	8	B=md.mesh.segments(:,2);
[9734]	9	Nx=-(md.mesh.y(A)-md.mesh.y(B));
	10	Ny= md.mesh.x(A)-md.mesh.x(B);
[9684]	11	Vx=(md.initialization.vx(A)+md.initialization.vx(B))/2;
	12	Vy=(md.initialization.vy(A)+md.initialization.vy(B))/2;
[6959]	13
	14	%dot product
	15	VdotN=Vx.Nx+Vy.Ny;
	16
[9725]	17	flag=zeros(md.mesh.numberofvertices,1);
[6959]	18	flag(A(find(VdotN>0)))=1;

Note: See TracBrowser for help on using the repository browser.