% plotADCPwavets.m
rot = 2;
zADCP = 2.25; % elevation of ADCP transducers (mab)
fn='wavets.txt';
if(exist('w')~=1),w=load(fn);,end

datenum_wave = datenum(w(:,1)+2000,w(:,2),w(:,3),w(:,4),w(:,5),w(:,6));
Hs = w(:,7); % Sig. height (m)
Tp = w(:,8); % Average period (s)
Dp = w(:,9)+rot; Dp = Dp-( Dp>=360. )*360.; % true direction (degrees)
Depth = w(:,10)+zADCP; % water depth (m)
cspd = w(:,11); % current speed (m/s)
cdir = w(:,12)+rot; cdir = cdir - (cdir>=360.)*360. ; % current direction

sqrt2 = sqrt(2);
omega = 2*pi./Tp;
kh = qkhf( omega, 9 );
do = (Hs./sqrt2)./sinh(kh);
ub = omega.*(do./2);
clear w;