MITgcm_contrib/tides/intp_line_safe.m

% 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 
1	% function utide = intp_line_safe(utide, lat_bc)
2	%
3	% Interpolate to all the BC grid
4	%
5	% Input: utide, lat_bc, utide has NaN
6	%
7	% Output: utide, does not have Nan
8	%
9	% XC WANG /08/12/2012
10	%
11	function utide = intp_line_safe(utide, lat_bc)
12
13	badp = ~isfinite(utide);
14	goodp = isfinite(utide) ;
15	utide(badp) = interp1(lat_bc(goodp), utide(goodp), lat_bc(badp), 'spline');
16
17	badp = ~isfinite(utide);
18	goodp = isfinite(utide) ;
19	utide(badp) = interp1(lat_bc(goodp), utide(goodp), lat_bc(badp), 'nearest') ;
20
21	return
22	end