% make_B6_wave - Build waves for ROMS at B6
%
% BF - 5/8/8

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

if ~exist('WAVEFNAME','var')
    WAVEFNAME ='B6_waves.nc';
end

load D:\bf\Data\PV_2004\data_proc2007\puv_proc_B6
duration=15;                    %run time in days
dt=mean(diff(jt));           %time step in days
wtime=(0:dt:duration)';         %define time vector

window=ones(6,1)/6;
Hrmsu=filtfilt(window,1,[PUV.Hrmsu]');
ubru=filtfilt(window,1,[PUV.ubr]');
Tr=filtfilt(window,1,[PUV.Tr]');
azr=filtfilt(window,1,[PUV.azr]');

Hwave=interp1(jt-jt(1),Hrmsu,wtime)/100;
Ub=interp1(jt-jt(1),ubru,wtime);
Pwave=interp1(jt-jt(1),Tr,wtime);
Dwave=interp1(jt-jt(1),azr,wtime);
Lwave = 10*ones(length(wtime),1);

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

%% write uniform forcing using
% ier=write_roms_wave_stress_forcing(GRDNAME,WAVEFNAME,wtime,...
%     Hwave,Dwave,Pwave,Lwave)
ier=write_roms_wave_ub_forcing(GRDNAME,WAVEFNAME,wtime,...
    Ub,Dwave,Pwave,Lwave)