% make_B6_storm_force - build the current stress for ROMS at B6
%
% BF 5/8/8

if ~exist('GRDNAME','var')
    GRDNAME = 'ocean_grd.nc';
end

if ~exist('WINDFNAME','var')
    WINDFNAME ='B6_storm.nc';
end

load D:\bf\STRATAFORM\patsmodel\sstm\PlotFigures\profiles_PCADP_6cells jtime us6 zo6 r26
duration=15;                %run time in days
dt=mean(diff(jtime))*4;       %time step in days
wtime=(0:dt:duration)';     %define time vector
jt=jtime-jtime(1);

ustlp=us6/100;
z0lp=zo6/100;
r2=r26;
ok = find( (ustlp>0)&(r2>.99)&(z0lp>0) );

ust = interp1(jt(ok),ustlp(ok),wtime,'nearest','extrap');

taux=ust.^2*1027;
tauy=0.*taux;

figure
plot(wtime,taux)
ylabel('taux')

%% write uniform forcing using
ier=write_roms_wind_stress_forcing(GRDNAME,WINDFNAME,wtime,taux,tauy)