%ov_morison2  its learning


%% Wave properties
zeta_a=16/2; % half wave height
T=10;

h=[60]; % Water depth 350 30 5

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

omega=2*pi/T;  % Vinkel acceleration
k=2*pi/lambda; % Reduced frequency, wave number

%% Particle motion
x=0;      % [m] x positive in direction of wave. Origin at origin of construction
t=0:T/(100-1):T; % [s] time
z=-30;    % [m] z positive up, origin at mean wave height
[u,w,udot,wdot]=particlemotion(t,x,omega,k,h,z,zeta_a,0);

figure(2)
plot(t,sin(omega.*t).^2,t,cos(omega.*t),t,sin(omega.*t).^2+cos(omega.*t))

%% Geometry of deplaced construction
D=1;
A=pi*(D/2)^2;
zb=h;  % Bottom position of displaced construction
zt=5;  % Top position of displaced construction 
return
%% Forces and moments on construction
dFD=CD*0.5*rho_s*D.*abs(u).*u;
dFM=CM*rho_s*Az.*udot;