% calculflux_ADCP.m
% sous programme appele par flux_ADCP.m
% donne le flux (en Sv) et le flux surfacique (en Sv.km^-2) qui traverse 
% une section

% 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
    FF=FF+r(ind1,ind2)*cos(alpha)*dz*dll*abs(cos(angs))/1000000;
    %surface en km^2
    SS=SS+dz*dll*abs(cos(angs))/1000000;
end

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