% read OK mesonet wind data files (for 5-minute intervals)

% v. 1.1 - updated 3-6-12

% v. 1.2 - updated 8-24-16

% v. 1.3 - updated 8-25-16

% download the data files (into the same directory as your program) using a web browser:
http://www.inscc.utah.edu/~krueger/5270/windclassfiles_meso1.tar
http://www.inscc.utah.edu/~krueger/5270/windclassfiles_meso2.tar 

% each file contains data for one year for one variable at 3 stations

% see sgp05okmX1.a1.pdf for metadata

% array sizes:
% 3x105120 (or 3x105408 for a leap year)

% station IDs:
% group 1: 92 44 52
% group 2: 29 108 50

% *** COPY CODE HERE ***

load groupone_wspd97.mat

% whos

% Name                 Size                  Bytes  Class     Attributes
% 
%   groupone_wspd97      3x105120            2522880  double 

% *** COPY CODE HERE ***

% find and count missing data (NaN)

% *** COPY CODE HERE ***

% annual average wind speed and standard deviation

% *** COPY CODE HERE ***

% mean wind power density: use density = 1 kg/m^3 for now  (Eq. 10-8)

% *** YOUR NEW CODE GOES HERE ***

% wind speed frequency distribution (for all directions)

% *** COPY CODE HERE ***

% mean daily wind speed versus day of year (plot)

% 1. calculate & store day of year for every array element: first array 
% element is 5-minute average ending at 00:05  (hh:mm), second at 00:10, etc

% 2. Use the above array to get mean wind speed for day.

% *** YOUR NEW CODE GOES HERE ***


% mean hourly wind speed versus hour of day (plot)

% 1. calculate & store hour of day for every array element: first array 
% element is 5-minute average ending at 00:05  (hh:mm), second at 00:10, etc

% 2. Use the above array to get mean wind speed for each hour of day.

% *** YOUR NEW CODE GOES HERE ***


% mean turbulence intensity (3-sec samples over 5-min interval) versus wind
% speed (Eq. 10-7)

% *** YOUR NEW CODE GOES HERE ***


% *** COPY CODE HERE ***