function y=omegalmc(x)

% OMEGALMC: estimates wind log profile correction due to surface waves.

% y=OMEGALMC(x) computes the log profile correction function due to wind

% distortion associated with surface waves. Input is x=za/Hw, where za

% is the measurement height and Hw is the dominant surface wave height. 

% Functional form is simplified (analytic) version of empirical omega

% curves shown in Fig. 9b of Large, Morzel, and Crawford (1995), J. Phys.

% Oceanog., 25, 2959-2971, with the wave-induced roughness length xr=0.15. 

% Assumes x is a vector with all elements greater than zr. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 3/8/97: version 1.0

% 8/5/99: version 2.0

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



xr=0.15;

ylimit=-6;

y=ylimit.*ones(size(x));

i=find(x<3.2967);

% polynomial fit

a=-2.6;

p1=-0.0199;

p2=0.0144;

p3=0.7660;

p4=0.0654;

x2=x(i).^2;x3=x2.*x(i);

y(i)=a.*log(x(i)./xr)+p1.*x3+p2.*x2+p3.*x(i)+p4;



  



