stressve

 STRESSVE: computes stress using Vera (1983) neutral drag law.
 tau = STRESSVE(sp,z,rhoa) computes the neutral wind stress given the wind
 speed at height z following Vera (1983) [see Large, Morzel, and Crawford
 (1995), J. Phys. Oceanog., 25, 2959-2971 (eqn. 8)]. Assumes z a fixed
 scalar. Air density is an optional input.

   INPUT:  sp    - wind speed  [m/s]
           z     - measurement height  [m]
           rhoa  - air density (optional)  [kg/m^3]

   OUTPUT: tau   - wind stress  [N/m^2]

0001 function tau=stressve(sp,z,rhoa)
0002 % STRESSVE: computes stress using Vera (1983) neutral drag law.
0003 % tau = STRESSVE(sp,z,rhoa) computes the neutral wind stress given the wind
0004 % speed at height z following Vera (1983) [see Large, Morzel, and Crawford
0005 % (1995), J. Phys. Oceanog., 25, 2959-2971 (eqn. 8)]. Assumes z a fixed
0006 % scalar. Air density is an optional input.
0007 %
0008 %   INPUT:  sp    - wind speed  [m/s]
0009 %           z     - measurement height  [m]
0010 %           rhoa  - air density (optional)  [kg/m^3]
0011 %
0012 %   OUTPUT: tau   - wind stress  [N/m^2]
0013 
0014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0015 % 3/8/97: version 1.0
0016 % 8/26/98: version 1.1 (revised by RP)
0017 % 4/2/99: version 1.2 (air density option added by AA)
0018 % 8/5/99: version 2.0
0019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0020 
0021 % load constants
0022 as_consts;
0023 
0024 if nargin == 2
0025   rhoa = rho_air; 
0026 end
0027 
0028 [cd,u10]=cdnve(sp,z);
0029 tau=rhoa*(cd.*u10.^2);
0030 
0031 
0032