% direction_princ_ADCP2.m
% sous-pg de flux_ADCP2.m
% calcule la direction moyenne des vitesses le long de la radiale
% fait appel à
%    - anglesection.m

% CALCUL DE L'ANGLE POUR DETERMINER LES COURANTS A GAUCHE DE LA SECTION
% BATEAU
if choixs==1 

    % asect = angle de la section totale /Nord (en deg)
    cd d:\cluster\matlab
    [asect,Dlonsect,Dlatsect]=anglesection(xyt(indt(1),1),xyt(indt(1),2),xyt(indt(end-1),1),xyt(indt(end-1),2));
    cd (ici)
    
    % direction orthogonale à la section bateau clockwise /Nord (en deg)
    directionsect=asect+90;
    
    %angle vitesse/direction orthogonale à la section bateau (en rad)
    alphasect=pi/180*(tt-directionsect);
end

% CALCUL DE LA DIRECTION MOYENNE DES COURANTS
tu=0;
moy=0;
ss=0;
for i=1:length(indt)-1
    ind2=indt(i);
    if choixs==1 % pour les courants à gauche de la section bateau
       indz=find(cos(alphasect(:,ind2))<0 & r(:,ind2)>Vmin);
    elseif choixs==2 % pour les courants ouest
       indz=find(u(:,ind2)<0 & r(:,ind2)>Vmin);
    elseif choixs==3 % pour les courants nord 
       indz=find(v(:,ind2)>0 & r(:,ind2)>Vmin);   
    elseif choixs==4 % aucun critère 
       indz=find(r(:,ind2)>Vmin);
    end
    %indz=find(u(:,ind2)>0 & r(:,ind2)>Vmin);
    inter=find(indz(:)>=indzc(1) & indz(:)<=indzc(end));
    indz=indz(inter);
    for p=1:length(indz)
        ind1=indz(p);
        tu=tt(ind1,ind2);
        moy=moy+tu;
        ss=ss+1;
    end
end
direction=moy/ss;

