%
% USAGE: [Tx, Ty]=ra_windstr(u,v)
%
% DESCRIPTION: Function to compute wind stress from wind field data Based on Gill, 1982
% Formula and a non-linear Cd based on Large and Pond (1981), modified for low wind
% speeds (Trenberth et al., 1990)
%
% INPUTS:
% u = Zonal wind component [m/s], must be 2D
% v = Meridional wind component [m/s], must be 2D
%
% OUTPUT:
% Tx = Zonal wind stress [N/m^2]
% Ty = Meridional wind stress [N/m^2]
%
% DISCLAIMER:
% Albeit this function is designed only for academic purpose, it can be implemented in
% research. Nonetheless, author does not guarantee the accuracy.
%
% REFERENCE:
% A.E. Gill, 1982, “Atmosphere-Ocean Dynamics”, Academy Press, Vol. 30.
% W. G. Large & S. Pond., 1981,“Open Ocean Measurements in Moderate to Strong Winds”,
% J. Physical Oceanography, Vol. 11, pp. 324 - 336.
% K.E. Trenberth, W.G. Large & J.G. Olson, 1990, “The Mean Annual Cycle in Global Ocean
% Wind Stress”, J.Physical Oceanography, Vol. 20, pp. 1742 – 1760.
% ***********************************************************************************************%
Ramkrushn Patel (2021). Wind Stresses computation (https://www.mathworks.com/matlabcentral/fileexchange/53391-wind-stresses-computation), MATLAB Central File Exchange. Retrieved .
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Create scripts with code, output, and formatted text in a single executable document.
Hello again and sorry for the double post,
There is no guideline on the height of the velocity in input. I am assuming it is 10m above sea level ?
Hi,
I agree with the previous comment (U > 1 && U <= 3) should be written instead of (U > 1 || U <= 3).
Let us know...
Regards,
Sam
there are two lines of this code that have an error:
on this line: else if (U > 1 || U <= 3)
it should say: else if (U > 1 && U <= 3)
on this line: else if (U > 3 || U < 10)
it should say: else if (U > 3 && U < 10)