% plotgolts

colormap=[0.9 0.9 0.9];

% bathy_digitalisee 1: lat; 2: long; 3 prof
% pas de 999999 entre les differents niveaux

load bathy0;
load bathy100;
load bathy200;
load bathy500;
load bathy1000;
load bathy1500;
load bathy2000;
load bathy2500;


golts=[5.13	43.07 0.];

%  Triangle externe
rr=[5.13	43.25,
   5.13		43.07,
   5.13		42.5,
   5.813   	43.04,
   5.13		43.25];

hold on
plot(bathy0(:,2),bathy0(:,1),'k','linewidth',2);
plot(bathy100(:,2),bathy100(:,1));
plot(bathy200(:,2),bathy200(:,1),'linewidth',2);
plot(bathy500(:,2),bathy500(:,1));
plot(bathy1000(:,2),bathy1000(:,1),'linewidth',2);
plot(bathy1500(:,2),bathy1500(:,1));
plot(bathy2000(:,2),bathy2000(:,1));
plot(bathy2500(:,2),bathy2500(:,1));
 
% trace de la terre d apres continenta
poly=zeros(length(bathy0)+3,2);
poly(1:length(bathy0),1)=bathy0(:,2);
poly(1:length(bathy0),2)=bathy0(:,1);
poly(length(bathy0)+1:length(bathy0)+3,:)=[bathy0(end,2),43.6;bathy0(1,2),43.6;bathy0(1,2),bathy0(1,1)];
brighten(0.8)
   patch(poly(:,1),poly(:,2),[0.5,0.5,0.5])
   axis([3 6 42 43.6]) 
 

%
asp=[1 cos(43*pi/180) 1];
set(gca,'ylim',[42.2  43.4],'xlim',[5 6],...
	'ytick',[42.5 42.666 42.833 43.0 43.166 43.333],...
	'yticklabel',['42.30';'42.40';'42.50';'43.00';'43.10';'43.20'],...
	'xtick',[5. 5.333 5.666 6.0],...
	'xticklabel',['5.00';'5.20';'5.40';'6.00'],...
        'fontsize',14,...
        'xgrid','on','ygrid','on','gridlinestyle','-',...
        'box','on',...
        'DataAspectRatioMode','manual','DataAspectRatio',asp);
 
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')


% plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
 

% text isobathes
 
text(5.6,42.7,'Isobathes')
text(5.6,42.65,'100 m')
text(5.6,42.62,'200 m')
text(5.6,42.59,'500 m')
text(5.6,42.56,'1000 m')
text(5.6,42.53,'1500 m')
text(5.6,42.50,'2000 m')

