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	dimitri	1.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