function [Hx,Hy]=gen_hxhy(He)
%
% Returns the depths at U and V points.
%
% e.g.
%      [Hx,Hy,xa,ya]=gen_hxhy(He,xe,ye);

Hx=max( He([end 1:end-1],:,:) , He );
Hy=max( He(:,[end 1:end-1],:) , He );