%% Solution du test2122 - Partie 1
%% Lecture des données et affectation à des listes
data= load("vel_11000.prof");

% Extraction de la colonne 2 : valeurs tabulées de y+
yp= data(:,2);

% Extraction de la colonne 13 : valeurs tabulées de S+= dU+/dy+
Sp= data(:,13);

% Extraction de la colonne 7 : valeurs tabulées de -tauxy+
% Ne pas oublier le signe - !
tauxyp= -data(:,7);

% Contrainte totale
taup= Sp + tauxyp;

%% Q1.2.2Figure1
figure(1)
hold on
plot( yp, Sp, "b","LineWidth",1);
plot( yp, tauxyp, "r","LineWidth",1);
plot( yp, taup, "k","LineWidth",1);
axis([0 2500 0 1])
xlabel("y+"); ylabel("tau+");
legend("dU+/dy+","tauxy+","tau+");
pbaspect([2.5 1 1])
nfig = "fig1test2122.eps";
print('-depsc',nfig); system("epstopdf " + nfig);

%% Q1.2.4Figure2v0 pour voir
delete(gca)
figure(2)
hold on
plot( yp, Sp, "b","LineWidth",2);
plot( yp, tauxyp, "r","LineWidth",2);
plot( yp, taup, "k","LineWidth",2);
yline(0.9, "k-","LineWidth",2)
axis([0 1000 0 1])
xlabel("y+"); ylabel("tau+");
legend("dU+/dy+","tauxy+","tau+");

%% Q1.2.4Calcul de ypg= 1ère valeur de yp pour laquelle taup < 0,9
% I= liste des valeurs d'indice-position dans la liste des taup pour
% lesquelles taup < 0,9
I= find(taup<0.9);

% indg= 1ère valeur dans I = 1er indice-position pour lequel le critère
% n'est plus vérifié
indg= I(1)

% Valeur de y+ correspondante, exacte puis arrondie à l’entier le + proche
ypg= yp(indg)
ypg= round(ypg)

% Verification : on passe bien en dessous de 0,9 au point d'indice indg
taupavant= taup(indg-1)
taupapres= taup(indg)

%% Q1.2.5Figure2 avec le ypg calculé
delete(gca)
figure(2)
hold on
plot( yp, Sp, "b","LineWidth",2);
plot( yp, tauxyp, "r","LineWidth",2);
plot( yp, taup, "k","LineWidth",2);
yline(0.9, "k--","LineWidth",2)
axis([0 ypg 0 1])
ax = gca; ax.FontSize = 13;
xlabel("y+", "FontSize",15); ylabel("tau+", "FontSize",15);
legend("dU+/dy+","tauxy+","tau+", "Location","southeast", "FontSize",15);
nfig = "fig2test2122.eps";
print('-depsc',nfig); system("epstopdf " + nfig);

%% End !
system('rmeps');