% Calculates wind speed extremes by Lieblein's BLUE for Mean Return 
% Intervals M from NIST database
% V_TH, V_NT = annual wind speed max for thunderstorms and non-thunderstorm
% speed = vector of extreme wind speeds with MRI M 
% M(i)does not have to be an integer 
% mu, sigma = location and scale parameters of Gumbel distribution
addpath('C:\windPRESSURE');
filename = 'Bakersfield_723840'; % Name of Excel file
shtname  = 'Sheet1';             % Name of sheet including data
[V_NT,V_TH] = AnnualWindSpeed(filename,shtname)
M = [10 20 25 50 100 300 700 1700 3000]
[speed, mu, sigma] = blue8(V_NT,M)
[speed, mu, sigma] = blue8(V_TH,M)

