spshftve

 SPSHFTVE: adjusts wind speed from z1 to z2 following Vera (1983).
 sp2 = SPSHFTVE(sp1,z1,z2) shifts the wind speed sp1 measured at z1 to
 z2 using the neutral drag law of Vera (1983) [see Large, Morzel, 
 and Crawford (1995), J. Phys. Oceanog., 25, 2959-2971 (eqn. 8)]. 
 Assumes z1 and z2 scalars.

   INPUT:  sp1 - measured wind speed  [m/s] 
           z1 - measurement height  [m]
           z2 - desired height of sp2  [m]

   OUTPUT: sp2 - predicted wind speed  [m/s]
           ustar - friction velocity  [m/s]

0001 function [sp2,ustar]=spshftve(sp1,z1,z2)
0002 % SPSHFTVE: adjusts wind speed from z1 to z2 following Vera (1983).
0003 % sp2 = SPSHFTVE(sp1,z1,z2) shifts the wind speed sp1 measured at z1 to
0004 % z2 using the neutral drag law of Vera (1983) [see Large, Morzel,
0005 % and Crawford (1995), J. Phys. Oceanog., 25, 2959-2971 (eqn. 8)].
0006 % Assumes z1 and z2 scalars.
0007 %
0008 %   INPUT:  sp1 - measured wind speed  [m/s]
0009 %           z1 - measurement height  [m]
0010 %           z2 - desired height of sp2  [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 CDNVE efficiently by RP)
0018 % 8/5/99: version 2.0
0019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0020 
0021 % set constants
0022 as_consts;
0023 
0024 % find cd and ustar
0025 [cd10,sp10]=cdnve(sp1,z1);
0026 
0027 ustar=sqrt(cd10).*sp10;
0028 
0029 sp2=sp10+ustar.*log(z2./10)/kappa;
0030