function [coef,nt]=stfn(dol0,hohb)
% **************************************************************
% [coef,nt]=stfn(dol0,hohb)
% **************************************************************
% Spline interpolate Dean's basic stream function coefficients at
%        dol0=d/L0=depth divided by deep water wave length
%        hohb=H/Hb=wave height divided by the breaking wave height
% Output:
% Col vector of stream function coefficients
% nt = number of terms in stream function series
% ***************************************************************
[S1,S2,S3,S4,S5,S6,S7,S8,S9,S10]=sfn1(dol0,hohb);
[S11,S12,S13,S14,S15,S16,S17,S18,S19]=sfn2(dol0,hohb);
nt=numcoef(dol0,hohb);
for j=1:19
   line=['coef(',int2str(j),',1)=interp2d(dol0,hohb,S',int2str(j),');'];
   eval(line)
end