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 |