load aver_ux_array.dat;

load z_array.dat;

r=z_array(:,);

r=r.'

r_j=0.00125;

r_nor=r/d;

ux_x_3d=aver_ux_array[,:];  % extract  arrow, which is X/D=

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