function [taux10,tauy10] = whiskWindNY(wjdtime,wu,wv,data_scale,plot_scale);
%   [taux10,tauy10] = whiskWindNY(time,WU_422,WV_423);
%
%Function which creates a whisker plot of wind stress given a time vector and
%the eastward and northward components of wind velocity.
%
%Soupy Alexander, 1/29/02

%Get rid of bad data points
wu = ridfill(wu);
wv = ridfill(wv);

%Low pass filter the component data
[wulp,wjdtimelp] = plfilt(wu,wjdtime,6);
[wvlp,wjdtimelp] = plfilt(wv,wjdtime,6);

%Compute the wind stress
[taux,tauy,u10,v10] = wstress(wulp,wvlp,5);

%Estimated 10 m wind stress
[taux10,tauy10,u10a,v10a] = wstress(u10,v10);

%Plot the wind stress
NYCurr_Whisk(taux10,tauy10,wjdtimelp,data_scale,plot_scale,[-4 4]);
gregaxd_fooled(wjdtimelp,plot_scale,5)