% plotbathy adds the bathymetry to the map plot window...

col=jet(100);

if(~exist('xb'))
    fname=MAP_FILE_UNIX;
    tt='error if you get here?';
  try
    xb = getnc(fname,'lon');
  catch
    pwd
    fname=MAP_FILE_PC;
    xb = getnc(fname,'lon');
  end
    yb = getnc(fname,'lat');
end

v=[xlimit ylimit];

if(handles.keys.map180)
    xb180=xb;
    gg=find(xb>180);
    xb180(gg)=-(360-xb(gg));
%    gg=find(xb<=180);
%    xb180(end+1:)=xb(gg);
    
    xb=xb180;    
end
ix = find(xb > v(1) & xb < v(2));
iy = find(yb > v(3) & yb < v(4));
hb = -1*getnc(fname,'height',[min(iy) min(ix)],[max(iy) max(ix)],[4 4]);
xb2 = xb(ix);
yb2 = yb(iy);

contourf(xb2(1:4:end),yb2(1:4:end),hb,[0:100:2000],'k');
%contourf(xb2,yb2,hb,[0:100:2000],'k')
caxis([0,2000]);
%colorbar

hold on