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

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

%Data from Newman:  pg 300;
w_data = [0    1    2   3    4    5    6];
a_data = [0.75 0.5 0.4 .35  .36  .37   .38];
b_data = [0.0  .31 .2  .15  .13  .1  .08];
 
a1 = interp1(a_data,CL_data,alpha,'spline'); 
a2 =