matlab code:
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