mlosch/interp_llc/drop_thinwalls.m

function [H]=drop_thinwalls(h,hx,hy,periodic)
%
% Returns the depths at T points based on T, U and V points.
%
% e.g.
%      H=drop_thinwalls(H2,Hx2,Hy2);

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

[nx,ny]=size(h);

% Fill-up cell-centers where all sides are higher
Hsides=min( hx, hx([2:end 1],:,:) );
Hsides=min( Hsides, hy );
Hsides=min( Hsides, hy(:,[2:end end],:) );
H=max( h, Hsides );

%H=h;
for jj=1:ny;
 jp=min(ny,jj+1);
 jm=max(1,jj-1);
 for ii=1:nx;
  if periodic==1
   ip=mod(ii,nx)+1;
   im=mod(ii-2+nx,nx)+1;
  else
   ip=min(nx,ii+1);
   im=max(1,ii-1);
  end

  hp=H(ii,jj);
  hm=H(im,jj);
  ho=hx(ii,jj);
  if ho>max(hm,hp)
   if hm<hp
    H(ii,jj)=ho;
   else
    H(im,jj)=ho;
   end
  end

  hp=H(ii,jj);
  hm=H(ii,jm);
  ho=hy(ii,jj);
  if ho>max(hm,hp)
   if hm<hp
    H(ii,jj)=ho;
   else
    H(ii,jm)=ho;
   end
  end

 end
end
1	function [H]=drop_thinwalls(h,hx,hy,periodic)
2	%
3	% Returns the depths at T points based on T, U and V points.
4	%
5	% e.g.
6	% H=drop_thinwalls(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	[nx,ny]=size(h);
13
14	% Fill-up cell-centers where all sides are higher
15	Hsides=min( hx, hx([2:end 1],:,:) );
16	Hsides=min( Hsides, hy );
17	Hsides=min( Hsides, hy(:,[2:end end],:) );
18	H=max( h, Hsides );
19
20	%H=h;
21	for jj=1:ny;
22	jp=min(ny,jj+1);
23	jm=max(1,jj-1);
24	for ii=1:nx;
25	if periodic==1
26	ip=mod(ii,nx)+1;
27	im=mod(ii-2+nx,nx)+1;
28	else
29	ip=min(nx,ii+1);
30	im=max(1,ii-1);
31	end
32
33	hp=H(ii,jj);
34	hm=H(im,jj);
35	ho=hx(ii,jj);
36	if ho>max(hm,hp)
37	if hm<hp
38	H(ii,jj)=ho;
39	else
40	H(im,jj)=ho;
41	end
42	end
43
44	hp=H(ii,jj);
45	hm=H(ii,jm);
46	ho=hy(ii,jj);
47	if ho>max(hm,hp)
48	if hm<hp
49	H(ii,jj)=ho;
50	else
51	H(ii,jm)=ho;
52	end
53	end
54
55	end
56	end