function sogbasemap(A)
% sogbasemap(A) draws a simple base map
% with gebco -1000m and dick's coastline,
% for the axes A. if A not given, default
% is used.
% 2/5/02

clf
orient tall

% set box size
Ad=[-75 -59 -70 -61.5]; % default
if nargin==0
   A=Ad;
end

% get data and merge
load wapcoast
load m1000

% plot 
h1=plot(xc,yc,'k');set(h1,'linewidth',2);hold on;
h2=plot(x1000,y1000,'k:');set(h2,'linewidth',2);
axis(A);
ym=(A(3)+A(4))./2;
aspect5(ym);
xlabel('Longitude','fontsize',12,'fontw','bold')
ylabel('Latitude','fontsize',12,'fontw','bold')

% add thick border line
nbox=1;
if nbox==1
   X1=A(1);X2=A(2);Y1=A(3);Y2=A(4);
   box=[X1 Y1;X2 Y1; X2 Y2;X1 Y2;X1 Y1];
   h3=plot(box(:,1),box(:,2),'k');set(h3,'linewidth',2) 
   clear box h3
end
