%% Solution du test2122 - Partie 2
%% 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);

% Viscosité turbulente adimensionnelle
nup= tauxyp./Sp

%% Q2.2.1 Figure3
figure(3)
hold on
plot( yp, nup, "k","LineWidth",1);
axis([0 2500 0 190])
xlabel("y+"); ylabel("nu+");
pbaspect([2.5 1 1])
nfig = "fig3test2122.eps";
print('-depsc',nfig); system("epstopdf " + nfig);

%% Q2.2.3 Figure4
delete(gca)
chiest= nup ./ yp;
figure(4)
hold on
plot( yp, chiest, "k","LineWidth",2);
axis([0 641 0 0.5])
ax = gca; ax.FontSize = 13;
xlabel("y+", "FontSize",15); ylabel("nu+/y+", "FontSize",15);
nfig = "fig4test2122.eps";
print('-depsc',nfig); system("epstopdf " + nfig);

%% End !
system('rmeps');