% direction_princ_ADCP.m
% sous-pg de flux_ADCP.m
% calcule la direction moyenne des vitesses le long de la radiale
tu=0;
moy=0;
ss=0;
% CALCUL DE LA DIRECTION MOYENNE DES COURANTS
for i=1:length(indt)-1
    ind2=indt(i);
    if choixs==1 % pour les courants ouest
       indz=find(u(:,ind2)<0 & r(:,ind2)>Vmin);
    else % sans 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;

