%
% load already gridded bathymetry  
%
  load bathgrid.dat
%
  blon=bathgrid(1:100,1)';
  blat=bathgrid(1:100:10000,2);
  bz=reshape(bathgrid(:,3),100,100)';
%
  clear bathgrid
%
% define bathymetry contour levels
%
%  cint=[60;200];
  cint=input(' Enter depth level: ');
  clab=contour(flipud(blon),blat,bz,cint);
%  clabel(clab);
