-
Notifications
You must be signed in to change notification settings - Fork 12
/
Copy pathwdnotes.m
37 lines (35 loc) · 1.7 KB
/
wdnotes.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
% WDNOTES: notes on estimating wave effects on wind measurements.
% The following set of mfiles can be used to correct the wind speed
% Ua measured at height za for the effects of the wave boundary layer
% following the empirical model presented by Large, Morzel, and
% Crawford (1995), J. Phys. Oceanog., 25, 2959-2971. In particular,
% an analytic expression was found for the omega function (omegalmc.m)
% shown in their Fig. 9b, which allows the 'true' wind speed (Ut10)
% and stress at 10m (assumed above the wave boundary layer height)
% to be computed using wavedist.m and the true wind speed (Uta) at the
% measurement height za using wavedis1.m. The Large et al model assumes
% neutral stability (reasonable for high winds and wave conditions)
% and uses a 10-m neutral drag law (cdnve.m) based on Vera (1983;
% unpublished manuscript). This drag law follows Large and Pond (1982)
% for winds above 10 m/s but increases at lower wind speeds like
% Smith (1987). The wave field is specified by the significant wave
% height Hw.
%
% To compute 'true' wind speed Uta at za given Hw, use
% Uta=wavedis1(Ua,za,Hw).
%
% To compute 'true' wind speed Ut at 10m given Hw, use
% [Ut10,(Ut10-U10)]=wavedist(Ua,za,Hw).
%
% To plot the predicted effects of wave distortion on the wind Ua
% measured at the height za for a range of significant wave heights
% Hw=[0:2:8] in m, use
% y=wavedis2(za).
%
% Subroutines called:
% y=omegalmc(x)
% cd10=cdnve(u10)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 7/28/99: version 1.1
% 8/5/99: version 2.0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%