g=delft3d_io_grd('read','island.grd')
d=delft3d_io_dep('read','island.dep',g,'location','cor')
x=d.cor.x;
y=d.cor.y;
zb=d.cor.dep;
figure(1)
plot(x,y,'k',x',y','k');
axis equal
hold on;
pcolor(x,y,zb)
colorbar
save 'island.x' x -ascii
save 'island.y' y -ascii
save 'island.zb' zb -ascii

