% test aspect ratio

cd e:\Gatti2\golts\Survey\Survey081203\ADCPTethys\tronca
load g1203tr1t3at.mat
lon1=xyt(:,1);
lat1=xyt(:,2);
u1=uv(1:2:end-1,:);
v1=uv(2:2:end,:);
load g1203tr3t10at.mat
lon2=xyt(:,1);
lat2=xyt(:,2);
u2=uv(1:2:end-1,:);
v2=uv(2:2:end,:);

figure
axis([4.8 5.9 42.2 43.45]);
%asp=[0.3 0.3 1];
asp=[1 cos(42.5*pi/180) 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 DES COURANTS
% rajout de l'échelle de courant
[i,j]=size(u1);
u1(5,j+1)=0.5;
v1(5,j+1)=0;
lon1(j+1)=5.6;
lat1(j+1)=43.24;
u1(5,j+2)=0;
v1(5,j+2)=0.5;
lon1(j+2)=5.6;
lat1(j+2)=43.24;
% rajout de l'échelle de courant
[i2,j2]=size(u2);
u2(5,j2+1)=0.5;
v2(5,j2+1)=0;
lon2(j2+1)=5.6;
lat2(j2+1)=43.22;
u2(5,j2+2)=0;
v2(5,j2+2)=0.5;
lon2(j2+2)=5.58;
lat2(j2+2)=43.24;


hold on

quiver(lon1,lat1,u1(5,:)*asp(1),v1(5,:)*asp(2),0,'.k')
quiver(lon2,lat2,u2(5,:)*asp(1),v2(5,:)*asp(2),0,'.k')

hold off

figure
subplot(211)
axis([4.8 5.9 42.2 43.3]);
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

%asp=[1 cos(42.5*pi/180) 1]
asp=[0.3 0.3 1]
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
ue=0.5;
ve=-0.5;
lone=5.4;
late=42.8;
quiver(lone,late,ue*asp(1),ve*asp(2),0)

subplot(212)
axis([4.8 5.9 42.2 43.3]);
% TRACET DU TRAIT DE COTE
cd e:\Gatti2\Maryse;
continenta
hold on

asp=[1 cos(42.5*pi/180) 1]
%asp=[0.3 0.3 1]
set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
ue=0.5;
ve=-0.5;
lone=5.4;
late=42.8;
quiver(lone,late,ue*asp(1),ve*asp(2),0)