% function utide = intp_line_safe(utide, lat_bc)  
%
% Interpolate to all the BC grid    
%
% Input:  utide, lat_bc, utide has NaN  
% 
% Output: utide, does not have Nan 
%
%  XC WANG /08/12/2012 
% 
function utide = intp_line_safe(utide, lat_bc)  

badp =  ~isfinite(utide);  
goodp =  isfinite(utide) ; 
utide(badp) = interp1(lat_bc(goodp), utide(goodp), lat_bc(badp), 'spline');   

badp =  ~isfinite(utide);  
goodp =  isfinite(utide) ; 
utide(badp) = interp1(lat_bc(goodp), utide(goodp), lat_bc(badp), 'nearest') ;     

return 
end