spshftlp

 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]

0001 function [sp2,ustar]=spshfttc(sp1,z1,z2)
0002 % SPSHFTTC: adjusts wind speed from z1 to z2 following Large&Pond (1981).
0003 % sp2 = SPSHFTLP(sp1,z1,z2) shifts the wind speed sp1 measured at z1 to
0004 % z2 using the neutral drag coefficient given the wind speed and air
0005 % temperature at height z following Large and Pond (1981), J. Phys. Oceanog.,
0006 % 11, 324-336.
0007 %
0008 %   INPUT:  sp1 - measured wind speed  [m/s]
0009 %           z1 - measurement height [m]
0010 %           z2 - desired height [m]
0011 %
0012 %   OUTPUT: sp2 - predicted wind speed  [m/s]
0013 %           ustar - friction velocity  [m/s]
0014 
0015 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0016 % 3/8/97: version 1.0
0017 % 8/27/98: version 1.1 (revised to use CDNLP efficiently by RP)
0018 % 8/5/99: version 2.0
0019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0020 
0021 % get constants
0022 as_consts;
0023 
0024 % find cd and ustar
0025 [cd10,sp10]=cdnlp(sp1,z1);
0026 
0027 ustar=sqrt(cd10).*sp10;
0028 
0029 sp2=sp10+ustar.*log(z2./10)/kappa;
0030