function depcoruv = depthcor_uv(depth, elev, k, omega )
% DEPTHCOR_UV - Linear transfer function for depth-attenuated velocities
% depcoruv = depthcor_uv(depth, elev, k, omega )
%
% Method from Guza and Thornton (1980)
%
% Input
%   depth  average water depth over burst [m]
%   elev   depth of sensor, measured positive up from still water level [m]
%   k      wavenumber [k = 2*pi/L, L in m]
%   omega  angular frequency velocity [omega = 2pi/T, T in s]
%
% Output 
%   depcoruv	depth correction

% Modified by Chris Sherwood, USGS
% Last modified 17 Sept 2001

kh = k*depth;
kz = k*elev;
omega = omega(:);
nf = length(omega);
depcoruv = ones(nf,1);

% change wavenumber at 0 Hz to NaN to avoid divide by zero
depcoruv(2:nf,1) = omega(2:nf,1) .* ( cosh(kz(2:nf,1)) ./ sinh(kh(2:nf,1)));