% calculflux_ADCP3.m
% sous programme appele par flux_ADCP3.m
% 04/06/07 : modif ajout de valeur absolue a cos(alpha)
% donne le flux (en Sv) et le flux surfacique (en Sv.km^-2) qui traverse 
% une section
% diffère de calculflux_ADCP.m dans la définition de l'angle alpha et il
% calcule aussi l'incertitude du calcul de flux

% calcul de angs (en rad), l'angle du transect local/surface de projection
% ie surface perpendiculaire au courant principal
angs=pi/180*(ang-angtr);

if isnan(angs)==0
    %calcul de la distance entre deux points consécutifs (ind2 et ind2+1)du
    %transect
    dll=sqrt(Dloni^2+Dlati^2);
    %flux en Sv
    % MODIF on rajoute valeur absolue au cos(alpha)
    FF=FF+r(ind1,ind2)*abs(cos(alpha(ind1)))*dz*dll*abs(cos(angs))/1000000;
    % MODIF FIN
    %surface en km^2
    SS=SS+dz*dll*abs(cos(angs))/1000000;
    %incertitude du calcul
    vinc=0.02;
    II=II+vinc*abs(cos(alpha(ind1)))*dz*dll*abs(cos(angs))/1000000;
    
    %trouver bin qui rentrent dans calcul de flux % modif 16/11/2007
    matc(ind1,ind2)=1;
end

%flux surfacique en Sv.km^-2
FU=FF/SS;
