% direction_princ_ADCP3.m
% sous-pg de flux_ADCP3.m
% modif de direction_princ_ADCP2 pour corriger calcul de la direction moy
% avec critères
% calcule la direction moyenne des vitesses le long de la radiale
% fait appel à
%    - anglesection.m

% CALCUL DE L'ANGLE POUR DETERMINER LES COURANTS 
if choixs==1 % A GAUCHE DE LA SECTION BATEAU
    % asect = angle de la section totale /Nord (en deg)
    cd g:\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)
    
    % diection 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);
    
    %% MODIF des angles pour avoir bon calcul du courant moyen
    indg=find(tt<=asect);
    tt(indg)=tt(indg)+360;
elseif choixs==3 % AU NORD 
    indn=find(tt<=90);
    tt(indn)=tt(indn)+360;
    %% MODIF FIN
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);
        if isnan(tt(ind1,ind2))==0 & isnan(r(ind1,ind2))==0
        tu=tt(ind1,ind2);
        moy=moy+tu;
        ss=ss+1;
        end
    end
end
direction=moy/ss;

