function [windTime, tauU, tauV] = mblt_windstress(dataDirectory);
%
%   [windTime, tauU, tauV] = mblt_windstress(dataDirectory);
%
%   Function which takes the concatenated wind data (lt_wind.nc) in the
%   dataDirectory of the MBLT report, insures that it is on a constant
%   hourly averaged time base, and calculates the low-pass filtered wind
%   stress.

%Soupy Alexander, 10/29/03

%Load in the wind data
windFile = netcdf(fullfile(dataDirectory, 'lt_wind.nc'), 'nowrite');
theTime = singlejd(windFile{'time'}(:), windFile{'time2'}(:));
windU = windFile{'WU_422'}(:);
windV = windFile{'WV_423'}(:);

%Set up the new time base & interpolate the data
newTime = min(theTime):1/24:max(theTime);
[windU_i] = smart_interp(theTime, windU, newTime, 5);
[windV_i] = smart_interp(theTime, windV, newTime, 5);
[windU_lp, windTime] = plfilt(windU_i, newTime, 6);
[windV_lp] = plfilt(windV_i, newTime, 6);

%Estimate the 10 m wind speed (buoy 44013 is at 5 m)
[tauX, tauY, windU_10, windV_10] = wstress(windU_lp, windV_lp, 5);
[tauU, tauV] = wstress(windU_10, windV_10);