wavdist2

 WAVDIST2: plots wave distortion effects on wind at za.
 WAVDIST2(za) plots the effects of wave distortion on the 
 wind Ua measured at height za for the following significant 
 wave heights Hw=[0:2:8] in m.

 INPUT:   za - wind measurement height [m]

0001 function y=wavdist2(za)
0002 % WAVDIST2: plots wave distortion effects on wind at za.
0003 % WAVDIST2(za) plots the effects of wave distortion on the
0004 % wind Ua measured at height za for the following significant
0005 % wave heights Hw=[0:2:8] in m.
0006 %
0007 % INPUT:   za - wind measurement height [m]
0008 
0009 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0010 % 5/5/97: version 1.0
0011 % 4/10/98: version 1.1
0012 % 8/5/99: version 2.0
0013 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0014 
0015 Hw=[0:2:8];
0016 Ua=[0.01:.01:20]';
0017 
0018 N=length(Hw);
0019 M=length(Ua);
0020 Ut=zeros(M,N);
0021 
0022 clg
0023 hold on
0024 for n=1:N
0025   Ut=wavdist1(Ua,za,Hw(n));
0026   plot(Ua,Ut)
0027 end
0028 
0029 title(['Predicted effects of wave distortion on wind speed at height ',num2str(za),' m'])
0030 xlabel('Measured wind speed Ua (m/s)')
0031 ylabel('Predicted wind speed Ut (m/s)')
0032 text(10,2,'Wave height increment = 2 m')
0033 grid
0034 
0035