gnuplot

set logscale x 10

set format x "10^{%L}"

set title 'Turbulence Boundary Layer'

set xlabel 'y^+'

set ylabel 'U^+'

set xrange [1:1e+04]

set grid

set arrow 1 from 5,1 to 5,18 nohead linetype 1 linewidth 2

set arrow 2 from 50,1 to 50,18 nohead linetype 1 linewidth 2

set arrow 3 from 600,1 to 600,26 nohead linetype 1 linewidth 2

show arrow

set label 1 "Visc." at 1.8,7

set label 2 "Buffer" at 15,7

set label 3 "Log." at 150,7

show label

set key top left

set key box

plot 'boundary_vlaw.dat' title 'Van Driest profile' with line

replot 'boundary_line.dat' title 'Visc.' with line

replot 'boundary_log.dat' title 'Log.' with line

replot 'Uplus-zplus_000.dat' title 'NO.1' with points pointtype 4

replot 'Uplus-zplus_067.dat' title 'NO.2' with points pointtype 5

replot 'Uplus-zplus_134.dat' title 'NO.3' with points pointtype 6

replot 'Uplus-zplus_201.dat' title 'NO.4' with points pointtype 7

replot 'Uplus-zplus_268.dat' title 'NO.5' with points pointtype 8