radial profiles of mean streamwise velocity at X/D=3

load aver_ux_array.dat;
load z_array.dat;
r=z_array(:,1);
r=r.‘
r_j=0.00125;
r_nor=r/d;
ux_x_3d=aver_ux_array[41,:];  % extract 41 arrow, which is X/D=3
ux_x_3d=ux_x_3d.‘;  % transpose
u0=1102.9;
ux_nor=ux_x_3d/u0;
ux_vs_r_3d=[r_nor ux_nor];
save -ascii ux_vs_r_3d_nor.dat ux_vs_r_3d;
plot(r_nor,ux_nor);
xlabel(‘R/R_{j}‘);
ylabel(‘U/U_{j}‘);
legend({‘X/D_{j}=3‘},‘Location‘,‘northeast‘);
print(‘ux_vs_r_3d‘,‘-depsc‘);  %Encapsulated PostScript® file
print(‘ux_vs_r_3d‘,‘-dpng‘);  %png image