%ar_1203.m
cd E:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
bin=3;

figure

subplot(231)
title('GOLTS Décembre 2003 - bin 3 (16m)')
load g1203tr0t1at.mat
lon=xyt(:,1);
lat=xyt(:,2);
u=uv(1:2:end-1,:);
v=uv(2:2:end,:);

load g1203tr1t3at.mat
lon2=xyt(:,1);
lat2=xyt(:,2);
u2=uv(1:2:end-1,:);
v2=uv(2:2:end,:);

axis([4.8 6.1 42.2 43.45]);
asp=[0.3 0.3 1];
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

% AJOUT DE LA BATHYMETRIE
%isobathes à partir de la bathymétrie réaliste
plotclaude
%TRACET DE LA CAMPAGNE GOLTS
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];

%  Stations sur la route externe
stat=[
	5.13	43.25
   5.13	43.07
   5.13	42.9275
   5.13	42.785
   5.13	42.6425
   5.13	42.5
5.813	43.04
5.6764	43.082
5.5398	43.124
5.4032	43.166
5.2666	43.208
5.13	43.25
];

%  Stations internes
statii=[
   5.3088   	42.6503
   5.5230  	 42.8047
    5.6613  	 42.8982
	5.7097	42.9544 
    5.2811   	42.7860
    5.3295   	42.9591
    5.2880  	 43.0667
 ];
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')
plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
plot(statii(:,1),statii(:,2),'xg',statii(:,1),statii(:,2),'og')

plot([stat(2,1),stat(11,1)],[stat(2,2),stat(11,2)],'g')
plot([stat(3,1),stat(10,1)],[stat(3,2),stat(10,2)],'g')
plot([stat(4,1),stat(9,1)],[stat(4,2),stat(9,2)],'g')
plot([stat(5,1),stat(8,1)],[stat(5,2),stat(8,2)],'g')

% TRACET DES COURANTS
quiver(lon(1:2:end),lat(1:2:end),u(bin,1:2:end)*asp(1),v(bin,1:2:end)*asp(1),0,'.k')
hold on
quiver(lon2(1:2:end),lat2(1:2:end),u2(bin,1:2:end)*asp(1),v2(bin,1:2:end)*asp(1),0,'.k')


subplot(232)
clear uv xyt lon lat u v u2 v2 lon2 lat2

cd E:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
load g1203tr0t1bt.mat
lon=xyt(:,1);
lat=xyt(:,2);
u=uv(1:2:end-1,:);
v=uv(2:2:end,:);
clear uv xyt

load g1203tr1t6at.mat
lon2=xyt(:,1);
lat2=xyt(:,2);
u2=uv(1:2:end-1,:);
v2=uv(2:2:end,:);


axis([4.8 6.1 42.2 43.45]);
asp=[0.3 0.3 1];
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

% AJOUT DE LA BATHYMETRIE
%isobathes à partir de la bathymétrie réaliste
plotclaude
%TRACET DE LA CAMPAGNE GOLTS
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];

%  Stations sur la route externe
stat=[
	5.13	43.25
   5.13	43.07
   5.13	42.9275
   5.13	42.785
   5.13	42.6425
   5.13	42.5
5.813	43.04
5.6764	43.082
5.5398	43.124
5.4032	43.166
5.2666	43.208
5.13	43.25
];

%  Stations internes
statii=[
   5.3088   	42.6503
   5.5230  	 42.8047
    5.6613  	 42.8982
	5.7097	42.9544 
    5.2811   	42.7860
    5.3295   	42.9591
    5.2880  	 43.0667
 ];
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')
plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
plot(statii(:,1),statii(:,2),'xg',statii(:,1),statii(:,2),'og')

plot([stat(2,1),stat(11,1)],[stat(2,2),stat(11,2)],'g')
plot([stat(3,1),stat(10,1)],[stat(3,2),stat(10,2)],'g')
plot([stat(4,1),stat(9,1)],[stat(4,2),stat(9,2)],'g')
plot([stat(5,1),stat(8,1)],[stat(5,2),stat(8,2)],'g')

% TRACET DES COURANTS
quiver(lon(1:2:end),lat(1:2:end),u(bin,1:2:end)*asp(1),v(bin,1:2:end)*asp(1),0,'.k')
hold on
quiver(lon2(1:2:end),lat2(1:2:end),u2(bin,1:2:end)*asp(1),v2(bin,1:2:end)*asp(1),0,'.k')


subplot(233)
clear u v uv xyt lon lat

cd E:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
load g1203tr4t1at.mat
lon=xyt(:,1);
lat=xyt(:,2);
u=uv(1:2:end-1,:);
v=uv(2:2:end,:);
clear uv xyt

axis([4.8 6.1 42.2 43.45]);
asp=[0.3 0.3 1];
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

% AJOUT DE LA BATHYMETRIE
%isobathes à partir de la bathymétrie réaliste
plotclaude
%TRACET DE LA CAMPAGNE GOLTS
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];

%  Stations sur la route externe
stat=[
	5.13	43.25
   5.13	43.07
   5.13	42.9275
   5.13	42.785
   5.13	42.6425
   5.13	42.5
5.813	43.04
5.6764	43.082
5.5398	43.124
5.4032	43.166
5.2666	43.208
5.13	43.25
];

%  Stations internes
statii=[
   5.3088   	42.6503
   5.5230  	 42.8047
    5.6613  	 42.8982
	5.7097	42.9544 
    5.2811   	42.7860
    5.3295   	42.9591
    5.2880  	 43.0667
 ];
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')
plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
plot(statii(:,1),statii(:,2),'xg',statii(:,1),statii(:,2),'og')

plot([stat(2,1),stat(11,1)],[stat(2,2),stat(11,2)],'g')
plot([stat(3,1),stat(10,1)],[stat(3,2),stat(10,2)],'g')
plot([stat(4,1),stat(9,1)],[stat(4,2),stat(9,2)],'g')
plot([stat(5,1),stat(8,1)],[stat(5,2),stat(8,2)],'g')

% TRACET DES COURANTS
quiver(lon(1:2:end),lat(1:2:end),u(bin,1:2:end)*asp(1),v(bin,1:2:end)*asp(1),0,'.k')


subplot(234)
clear u v uv xyt lon lat

cd E:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
load g1203tr1t5at.mat
lon=xyt(:,1);
lat=xyt(:,2);
u=uv(1:2:end-1,:);
v=uv(2:2:end,:);

axis([4.8 6.1 42.2 43.45]);
asp=[0.3 0.3 1];
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

% AJOUT DE LA BATHYMETRIE
%isobathes à partir de la bathymétrie réaliste
plotclaude
%TRACET DE LA CAMPAGNE GOLTS
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];

%  Stations sur la route externe
stat=[
	5.13	43.25
   5.13	43.07
   5.13	42.9275
   5.13	42.785
   5.13	42.6425
   5.13	42.5
5.813	43.04
5.6764	43.082
5.5398	43.124
5.4032	43.166
5.2666	43.208
5.13	43.25
];

%  Stations internes
statii=[
   5.3088   	42.6503
   5.5230  	 42.8047
    5.6613  	 42.8982
	5.7097	42.9544 
    5.2811   	42.7860
    5.3295   	42.9591
    5.2880  	 43.0667
 ];
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')
plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
plot(statii(:,1),statii(:,2),'xg',statii(:,1),statii(:,2),'og')

plot([stat(2,1),stat(11,1)],[stat(2,2),stat(11,2)],'g')
plot([stat(3,1),stat(10,1)],[stat(3,2),stat(10,2)],'g')
plot([stat(4,1),stat(9,1)],[stat(4,2),stat(9,2)],'g')
plot([stat(5,1),stat(8,1)],[stat(5,2),stat(8,2)],'g')

% TRACET DES COURANTS
quiver(lon(1:2:end),lat(1:2:end),u(bin,1:2:end)*asp(1),v(bin,1:2:end)*asp(1),0,'.k')

subplot(235)
clear u v uv xyt lon lat

cd E:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
load g1203tr5t1at.mat
lon=xyt(:,1);
lat=xyt(:,2);
u=uv(1:2:end-1,:);
v=uv(2:2:end,:);

load g1203tr1t0at.mat
lon2=xyt(:,1);
lat2=xyt(:,2);
u2=uv(1:2:end-1,:);
v2=uv(2:2:end,:);

axis([4.8 6.1 42.2 43.45]);
asp=[0.3 0.3 1];
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

% AJOUT DE LA BATHYMETRIE
%isobathes à partir de la bathymétrie réaliste
plotclaude
%TRACET DE LA CAMPAGNE GOLTS
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];

%  Stations sur la route externe
stat=[
	5.13	43.25
   5.13	43.07
   5.13	42.9275
   5.13	42.785
   5.13	42.6425
   5.13	42.5
5.813	43.04
5.6764	43.082
5.5398	43.124
5.4032	43.166
5.2666	43.208
5.13	43.25
];

%  Stations internes
statii=[
   5.3088   	42.6503
   5.5230  	 42.8047
    5.6613  	 42.8982
	5.7097	42.9544 
    5.2811   	42.7860
    5.3295   	42.9591
    5.2880  	 43.0667
 ];
 
plot(golts(:,1),golts(:,2),'gx')
plot(rr(:,1),rr(:,2),'--')
plot(stat(:,1),stat(:,2),'xb',stat(:,1),stat(:,2),'ob')
plot(statii(:,1),statii(:,2),'xg',statii(:,1),statii(:,2),'og')

plot([stat(2,1),stat(11,1)],[stat(2,2),stat(11,2)],'g')
plot([stat(3,1),stat(10,1)],[stat(3,2),stat(10,2)],'g')
plot([stat(4,1),stat(9,1)],[stat(4,2),stat(9,2)],'g')
plot([stat(5,1),stat(8,1)],[stat(5,2),stat(8,2)],'g')

% TRACET DES COURANTS
quiver(lon(1:2:end),lat(1:2:end),u(bin,1:2:end)*asp(1),v(bin,1:2:end)*asp(1),0,'.k')
hold on
quiver(lon2(1:2:end),lat2(1:2:end),u2(bin,1:2:end)*asp(1),v2(bin,1:2:end)*asp(1),0,'.k')