Parent Directory
|
Revision Log
|
Revision Graph
initial checkin of topography and hydrography interpolation scripts for the llc-grid, based on old matlab scripts by Alistair Adcroft Let's hope, they are useful.
1 | function [h,hx,hy]=fill_potholes(h,hx,hy) |
2 | % |
3 | % Fills "pot holes": holes deeper than all neighbouring sides |
4 | % |
5 | % e.g. |
6 | % [H2,Hx2,Hy2]=fill_potholes(H2,Hx2,Hy2); |
7 | |
8 | if size(h)~=size(hx) | size(h)~=size(hy) |
9 | error('Arguments must all be the same size!\n'); |
10 | end |
11 | |
12 | % Fill-up cell-centers where all sides are higher |
13 | Hsides=min( hx, hx([2:end 1],:,:) ); |
14 | Hsides=min( Hsides, hy ); |
15 | Hsides=min( Hsides, hy(:,[2:end end],:) ); |
16 | h=max( h, Hsides ); |
17 |
ViewVC Help | |
Powered by ViewVC 1.1.22 |