function [sp2,ustar]=spshfttc(sp1,z1,z2)
% SPSHFTTC: adjusts wind speed from z1 to z2 following Large&Pond (1981).
% sp2 = SPSHFTLP(sp1,z1,z2) shifts the wind speed sp1 measured at z1 to
% z2 using the neutral drag coefficient given the wind speed and air
% temperature at height z following Large and Pond (1981), J. Phys. Oceanog.,
% 11, 324-336.
%
% INPUT: sp1 - measured wind speed [m/s]
% z1 - measurement height [m]
% z2 - desired height [m]
%
% OUTPUT: sp2 - predicted wind speed [m/s]
% ustar - friction velocity [m/s]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 3/8/97: version 1.0
% 8/27/98: version 1.1 (revised to use CDNLP efficiently by RP)
% 8/5/99: version 2.0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% get constants
as_consts;
% find cd and ustar
[cd10,sp10]=cdnlp(sp1,z1);
ustar=sqrt(cd10).*sp10;
sp2=sp10+ustar.*log(z2./10)/kappa;