function [cb,hb]=topography(area)
% TOPOGRAPHY	Map with background bathymetry
%
% hb = topography(area)
%
% area	= [lonmin lonmax latmin latmax]
% 
% hb	= handle to bathymetic patches. 
%	  Use this to change appearance afterwards (clim, etc.) 
%
% See also M_MAP

%Time-stamp:<Last updated on 01/10/25 at 10:10:37 by even@gfi.uib.no>
%File:<d:/home/matlab/topography.m>

if nargin<1|isempty(area), area=[-20 12 60 72]; end

m_proj('Albers Equal-Area Conic',...
       'lon',area(1:2),'lat',area(3:4),'rec','off');  
[cb,hb]=m_tbase('contourf',fliplr(-[0:100:1000 1500:500:5000]),...
	'EdgeColor','none');
%[cb,hb]=m_tbase('contourf',fliplr(-[200:100:3000 3500:500:5000]),...
%	'EdgeColor',[.7 .7 .7]);
% [bz,blo,bla]=m_tbase(area); 
% bz=bz(1:2:end,1:2:end);blo=blo(1:2:end,1:2:end);bla=bla(1:2:end,1:2:end);
% [blo,bla]=m_ll2xy(blo,bla);
% hold on; hb=surf(blo,bla,bz); hold off;
colormap gray; cm=colormap; colormap(flipud(cm));
%ecolorbar(fliplr(-[0:100:1000 1500:500:3500]),...
%	  [.78 .32 .02 .6],'Bottom depth [m]',[],[],'EdgeColor','none');
caxis([-3500 0]); brighten(.8);% shading interp; 
%m_gshhs_l('patch','k'); hl=findobj(gcf,'tag','m_gshhs_l');
m_grid; hmg=findobj(gcf,'tag','m_grid_color');
% set(hmg,'visible','off');
%owsm;

%caxis(-[5000 0]); 

%m_gshhs_l patch k