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 |