% Oving bolge
clear all, close all,  clc
format compact, format shortg
marineconst; % laster inn konstanter fra annen fil marineconst.m

zeta_a=6/2; % half wave height
T=9;

hv=[5 30 350]; % Water depth 350 30 5

for i=1:length(hv)
h=hv(i);

[lambda,depth]=lambda_cal(T,h)

omega=2*pi/T;  % Vinkel acceleration
k=2*pi/lambda;

z=[0 -0.5*h -0.95*h]; % z positiv opp
x=0; % lar bolgen gå forbi, (t=0 fryser bolgen)

t=T/4; %omega=2*pi/T % Horisontal acceleration storst ved maks vertikal hastighet
[u(i,:),wnan(i,:),udotnan(i,:),wdot(i,:)]=particlemotion(t,x,omega,k,h,z,zeta_a,0);
t=0; %omega=2*pi/T % Vertikal acceleration storst ved maks horisontal hastighet
[unan(i,:),w(i,:),udot(i,:),wdotnan(i,:)]=particlemotion(t,x,omega,k,h,z,zeta_a,0);
end
u
w