function [taux,tauy]=wstress2(WINDX,WINDY)
%
% WSTRESS2  converts u and v wind velocity components to 
%  stress using the Garrett formulation.
%
  if(nargin==1),WINDY=0.,end;
  WINDMAG=sqrt(WINDX.*WINDX+WINDY.*WINDY);
            WDRAGCO = min(0.001*(0.75+0.067*WINDMAG),0.003);
            taux=0.001293*WDRAGCO.*WINDX.*WINDMAG*1e4;
            tauy=0.001293*WDRAGCO.*WINDY.*WINDMAG*1e4;
