function x=plotland(col,lt);
% PLOTLAND (COL,LT)   plot world map
%           COL       colour, defaults to .8*[1 1 1]
%                     or specify a linetype, such as '-' for a contour plot
%               LT    line thickness defaults to 2

if nargin < 1, col=.8*[1 1 1]; end

load topo,
ix=find(topo>=0); topo(ix)=ones(size(ix));
ix=find(topo<0); topo(ix)=zeros(size(ix));
lon=0.5:359.5; lat=-89.5:89.5;

hold on
if isstr(col)
  if nargin < 2, lt=2; end
  contour(lon,lat,topo,[1 1],col);
else
  shadeland(lon,lat,topo,col);
end
clear ix topo topomap1 topomap2
xlabel('Longitude East'),ylabel('Latitude North'),

%$$$ load topo,
%$$$ colormap([1 1 1; .5 .5 .5])
%$$$ contourf(0:359,-89:90,topo,[0 0]);
%$$$ clear topo topomap1 topomap2
%$$$ xlabel('Longitude East'),ylabel('Latitude North'),

%$$$ % alternative code
%$$$ clf reset
%$$$ pcolor(0:359,-89:90,topo),colormap(gray(2)),shading('interp')
%$$$ 
%$$$ % Kirill's method
%$$$ clf reset,colormap gray
%$$$ axis([0 360 -90 90]),hold on
%$$$ im=image('xdata',[0 360],'ydata',[-90 90]);
%$$$ tmp=interp2(0:359,-89:90,topo,0:.5:359,-89:.5:90);
%$$$ cdat = (tmp<0)+64*(tmp>=0);
%$$$ set(im,'cdata',cdat)
