stresslp

 STRESSLP: computes neutral wind stress following Large and Pond (1981).
 tau = STRESSLP(sp,z,rhoa) computes the neutral wind stress given the wind
 speed at height z following Large and Pond (1981), J. Phys. Oceanog.,
 11, 324-336. Air density is an optional input, otherwise assumed to
 be constant (1.22 kg/m^3). 

   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=stresslp(sp,z,rhoa)
0002 % STRESSLP: computes neutral wind stress following Large and Pond (1981).
0003 % tau = STRESSLP(sp,z,rhoa) computes the neutral wind stress given the wind
0004 % speed at height z following Large and Pond (1981), J. Phys. Oceanog.,
0005 % 11, 324-336. Air density is an optional input, otherwise assumed to
0006 % be constant (1.22 kg/m^3).
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 (optional air density added by AA)
0018 % 8/5/99: version 2.0
0019 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0020 
0021 % get constants
0022 as_consts;
0023 
0024 if nargin == 2
0025   rhoa = rho_air; 
0026 end
0027 
0028 [cd,u10]=cdnlp(sp,z);
0029 tau=rhoa.*(cd.*u10.^2);
0030