close all
clear all

Hs = 2;

Type = 1;  %Type refers monochromatic(Type = 1) or PM spectrum (Type = 0) wave
A = Hs/2;
%wind_spd = sqrt(Hs*g/0.21); %Wind appropriate to wave-height in PM spectrum
wave = IncWave;
wave = set(wave,'A',A);
wave = set(wave,'Type',Type);
wave = set(wave,'beta',-1);

disp('Computing incident wave profile');
t= 0:.1:100;
eta0 = eta(wave,0,t);

figure
plot(t,eta0);
