%carte_courant_old.m
% trace les courants mesurés par l'ADCP de coque à une profondeur choisie le long
% d'un nombre de transects choisis, sur une carte du GdL avec bathy et trajet GOLTS.
% fait appel à :
%    - continenta.m  pour tracer le trait de cote
%    - plotclaude.m  pour tracer les isobathes

clear
ici=pwd;
path_mat='d:\golts\Survey\Survey110602\ADCPTethys\tronca';
%path_mat='e:\Gatti2\golts\Valor\Valor070405\ADCPTethys\transect';

%choix des transects à tracer
choix=input('combien de transects à représenter? choix= ');
for ii=1:choix
    eval(['[f_mat',int2str(ii),',p_mat',int2str(ii),']=uigetfile([path_mat,''\*.mat'']);'])
    eval(['path_mat_file',int2str(ii),'=[p_mat',int2str(ii),',f_mat',int2str(ii),'];'])
    eval(['load (path_mat_file',int2str(ii),')'])
    eval(['lon',int2str(ii),'=xyt(:,1);'])
    eval(['lat',int2str(ii),'=xyt(:,2);'])
    eval(['u',int2str(ii),'=uv(1:2:end-1,:);'])
    eval(['v',int2str(ii),'=uv(2:2:end,:);'])
end

% Liste de toutes les profondeurs :
for jj=8:4:244
    if jj<=9
        p(jj/4-1,1:3)=sprintf('00%d',jj);
    elseif jj>=10 & jj<=99
        p(jj/4-1,1:3)=sprintf('0%d',jj);
    elseif jj>=100 & jj<=999
        p(jj/4-1,1:3)=sprintf('%d',jj);
    end
end
p=cellstr(p);

% choix de la profondeur à choisir
boucle=0;
kk=1;
while(boucle==0)
    bin=menu('Quelle profondeur ?',p);
    eval(['bin',int2str(kk),'=bin;'])
    eval(['prof',int2str(kk),'=(bin',int2str(kk),'+1)*4;'])
    
    figure
    axis([4.8 6.1 42.2 43.45]);
    asp=[0.3 0.3 1];
    set(gca,'DataAspectRatioMode','manual','DataAspectRatio',asp);
    
    % TRACER DU TRAIT DE COTE
    cd d:\Maryse;
    continenta
    hold on
    
    % AJOUT DE LA BATHYMETRIE
    %isobathes à partir de la bathymétrie réaliste
    plotclaude
    %TRACER 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.3295   	42.9591
    ];
    
    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')
    
    % TRACER DES COURANTS
    % rajout de l'échelle de courant
    [i,j]=size(u1);
    eval(['u1(bin',int2str(kk),',j+1)=0.5;'])
    eval(['v1(bin',int2str(kk),',j+1)=0;'])
    lon1(j+1)=5.6;
    lat1(j+1)=43.24;
    eval(['u1(bin',int2str(kk),',j+2)=0;'])
    eval(['v1(bin',int2str(kk),',j+2)=0.5;'])
    lon1(j+2)=5.6;
    lat1(j+2)=43.24;
    
    hold on
    for ii=1:choix
        eval(['quiver(lon',int2str(ii),',lat',int2str(ii),',u',int2str(ii),'(bin',int2str(kk),',:)*asp(1),v',int2str(ii),'(bin',int2str(kk),',:)*asp(2),0,''.k'')'])
         %eval(['quiver(lon',int2str(ii),'(300),lat',int2str(ii),'(300),u',int2str(ii),'(bin',int2str(kk),',300)*asp(1),v',int2str(ii),'(bin',int2str(kk),',300)*asp(2),0,''.r'')'])
        eval(['text(6.2,43.2-0.1*',int2str(ii),',''',eval(['f_mat',int2str(ii)']),''')'])
    end
    %hold off
    
    eval(['title([''GOLTS Juin 2002 - bin ',eval(['int2str(bin',int2str(kk),')']),' ( ',eval(['int2str(prof',int2str(kk),')']),' m)''])'])
    text(5.65,43.35,'0.5 ms^{-1}')
    
    choix5=menu('voulez vous visualiser une autre profondeur?','non','oui');
    if choix5==1
        boucle=1;
    elseif choix5==2
        kk=kk+1;
    end   
end % boucle while pour tracer à différentes profondeurs

cd (ici)