function [a, b] = sphere_heave_radiation_coeffs(d,w,g,rho)

% Returns added mass and damping coefficients for half-submerged
% sphere of diameter d. 

V = pi*d^3/12;

%Data from Newman:  pg 300;
w_data = [0  0.2  1    2   3    4    5    6];
a_data = [0.7 0.75 0.55 0.37 .36  .365  .37   .38];
b_data = [0.0 0.1 .31 .22  .14  .12  .11  .1];

a = rho*V*interp1(w_data,a_data,w.^2*d/g,'spline'); 
b = rho*V*w.*interp1(w_data,b_data,w.^2*d/g,'spline'); 
