% make_wind - Script to make up wind forcing and write it
   
   % start w/ wind, then wind stops
   if(exist('GRDNAME')~=1),GRDNAME = 'ocean_grd.nc';end
   if(exist('WINDFNAME')~=1),WINDFNAME ='ocean_frc.nc';end
   wtime = [0; 1/24; 18/24; 19/24;  1; 10]; % days
   tauy =  [0;  0.2;  .2;       0;       0;  0]; % N/m2
   taux = 0.*tauy;
   % write uniform forcing using
   ier=write_roms_wind_stress_forcing(GRDNAME,WINDFNAME,wtime,taux,tauy)