% hauteur d'eau amont
%
% Auteurs : Thibault Prévost.
% NOMA(s)   :0525 1500.
% Date    : 06/05/2021

clc;
clear all;
close all;

k=0;
Moyenne=0;

for j=1:20

%récupération des données de post traitements 
CheminDonnees=strcat(fullfile('C:\Users\Admin\Desktop\memoire\MatlabMesure hauteur\WEIR partie 1\WEIR0\weir0i\postProcessing\linesample',num2str(2*j),'data_alpha.water.xy'));
Donnees=load(CheminDonnees);

y=Donnees(:,2);
alpha=Donnees(:,4);

y2=Donnees(:,2);
alphaWater2=Donnees(:,4);

%calcul de l'intégral pour connaitre la hauteur de l'eau
Integral=sum(alphaWater2);
Hauteur2(j)=y2(round(Integral))+ (Integral-round(Integral))*(y2(round(Integral))-y2(round(Integral-1)));

%je prend la hauteur moyenne une fois l'écoulement stabilisé car il y a
%toujours de petites turbulences
if j>14 && j<=20
    %     if Hauteur2(j)-Hauteur2(j-1)<0.01
    Moyenne=Moyenne+Hauteur2(j)
    %     end
    k=k+1;
end
end

Hprint2=Moyenne/k;

figure
plot(Hauteur2)
legend(['Hauteur finale ', num2str(Hprint2),'[m]'],'Location','southeast')