function [vm] = Heave_VM(omega,m,rho,g,S,w,a33)

am = interp1(w,a33,omega);
vm = -(m+am).*omega.^2+rho*g*S;

return;

end


