%
function [y]=spladd(x1,y1,timevec,y);
%
%  function [y]=spladd(x1,y1,timevec,y);
%
%  purpose: to spline fit a series to a common time vector
%  first checking to rid of NaNs in the original series
%  and then adding the splinefitted series to a common matrix called "y"
%
 
    id=find(isnan(y1)~=1);
    sy1=spline(yd(id),y1(id),timevec);
    y=[y sy1];
