mlosch/interp_llc/remove_inlets.m

function [h,hx,hy]=remove_inlets(h,hx,hy)
%
% Remove simple inlets
%
% e.g.
%      [H2,Hx2,Hy2]=remove_inlets(H2,Hx2,Hy2);

if size(h)~=size(hx) | size(h)~=size(hy)
 error('Arguments must all be the same size!\n');
end

oh=h+1;

while ( sum(oh(:)-h(:)) )
oh=h;

% Find second deepest side
hh(1,:,:,:)=hx;
hh(2,:,:,:)=hx([2:end 1],:,:);
hh(3,:,:,:)=hy;
hh(4,:,:,:)=hy(:,[2:end end],:);
hh=sort(hh);

% Remove inlets
h=squeeze(hh(2,:,:,:));

% Regenerate side walls
hX=max(h,h([end 1:end-1],:,:)); hx=max(hx,hX);
hY=max(h,h(:,[1 1:end-1],:)); hy=max(hy,hY);

end
1	function [h,hx,hy]=remove_inlets(h,hx,hy)
2	%
3	% Remove simple inlets
4	%
5	% e.g.
6	% [H2,Hx2,Hy2]=remove_inlets(H2,Hx2,Hy2);
7
8	if size(h)~=size(hx) \| size(h)~=size(hy)
9	error('Arguments must all be the same size!\n');
10	end
11
12	oh=h+1;
13
14	while ( sum(oh(:)-h(:)) )
15	oh=h;
16
17	% Find second deepest side
18	hh(1,:,:,:)=hx;
19	hh(2,:,:,:)=hx([2:end 1],:,:);
20	hh(3,:,:,:)=hy;
21	hh(4,:,:,:)=hy(:,[2:end end],:);
22	hh=sort(hh);
23
24	% Remove inlets
25	h=squeeze(hh(2,:,:,:));
26
27	% Regenerate side walls
28	hX=max(h,h([end 1:end-1],:,:)); hx=max(hx,hX);
29	hY=max(h,h(:,[1 1:end-1],:)); hy=max(hy,hY);
30
31	end