%
% load already gridded bathymetry  
%
  load bathgrid.dat
%
  blon=bathgrid(1:100,1)';
  blat=bathgrid(1:100:10000,2);
  id9=find(bathgrid(:,3)>100);
  for kk=1:length(id9);
    bathgrid(id9(kk),3)=NaN;
  end;
  bz=reshape(bathgrid(:,3),100,100)';
%
  clear bathgrid
%
% define bathymetry contour levels
%
   path(path,'/home3/ocn/jmanning/m/extc')
   cint=[40:10:100];
%  cint=input(' Enter depth level: ');
   econtour(flipud(blon),blat,bz,cint,'linewidth',0.1,'fill');
   set(gca,'aspectratio',[1 .78],'box','on');
   colormap(gray)