#set term postscript color enhanced
#set output "sigma.ps"
set multiplot
set size 0.5,0.25
set origin 0.,0.
set ylabel "<x>[mm]"
set xlabel "turn"
set yrange [-20:20]
set ytics -20,10,20
set xrange [0:]
plot "fort.16" u 2:($3*1000) w l notitle
#
set origin 0.5,0.
set ylabel "sigma_x[mm]"
set yrange [0:50]
set ytics 0, 10,50
plot "fort.16" u 2:($6**.5*1000) w l not
#
set origin 0.,0.25
unset xlabel
set ylabel "<y>[mm]"
set yrange [-4:4]
set ytics -4,2,4
plot "fort.16" u 2:($4*1000) w l not
#
set origin 0.5,0.25
set ylabel "sigma_y[mm]"
set yrange [10:25]
set ytics 10,5,25
set autoscale
set xrange [0:]
plot "fort.16" u 2:($9**.5*1000) w l not 
#
set origin 0.,0.50
unset ylabel
set ylabel "<time>[us]"
incr = GPVAL_X_MAX *0.15/4
set autoscale
set yrange [0:]
set ytics 0,incr
set xrange [0:]
plot "fort.16" u 2:($5*1.e6) w l not
#
set origin 0.5,0.50
set ylabel "length[ns]"
set autoscale
#set ytics 4
plot "fort.16" u 2:($12**.5*1e9) w l not 
#

set autoscale
set xrange [1:]
set origin 0.,0.75
unset ylabel
set ylabel "<dp/p>[%]"
set ytics autofreq
plot "fort.16" u 2:($15*100) w l not 
#
set ylabel "sigma_p [%]"
set origin 0.5,0.75
plot "fort.16" u 2:($14**0.5*100) w l not