Context Navigation

outflow.m@ 9714

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

Rev	Line
[7512]	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);
[6959]	9	Nx=-(md.y(A)-md.y(B));
	10	Ny= md.x(A)-md.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
[8298]	17	flag=zeros(md.numberofnodes,1);
[6959]	18	flag(A(find(VdotN>0)))=1;

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