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]

- as_consts AS_CONSTS: returns values of many constants used in AIR_SEA TOOLBOX
- cdnlp CDNLP: computes neutral drag coefficient following Large&Pond (1981).

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

Generated on Thu 26-Apr-2007 12:23:13 by