set term postscript color enhanced
set output "sigma.ps"
set multiplot
set size 0.5,0.34
set origin 0.,0.
set ylabel "<x>[mm]"
set xlabel "turn"
set yrange [-20:20]
set ytics -20,10,20
set xrange [0:]
plot "Beam_moments_by_element.dat" u 2:($3*1000) w l notitle
#
set origin 0.5,0.
set ylabel "{/Symbol s}_x[mm]"
set yrange [0:30]
set ytics 0, 5,30
plot "Beam_moments_by_element.dat" u 2:($6**.5*1000) w l not
#
set origin 0.,0.34
set size 0.5,0.32
unset xlabel
set ylabel "<y>[mm]"
set yrange [-10:10]
set ytics -10,5,10
plot "Beam_moments_by_element.dat" u 2:($4*1000) w l not
#
set origin 0.5,0.34
set ylabel "{/Symbol s}_y[mm]"
set yrange [0:30]
set ytics 0,5,30
set xrange [0:]
plot "Beam_moments_by_element.dat" u 2:($9**.5*1000) w l not 
#exit
#
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 "Beam_moments_by_element.dat" u 2:($5*1.e6) w l not
#
set origin 0.5,0.50
set ylabel "length[ns]"
set autoscale
#set ytics 4
plot "Beam_moments_by_element.dat" 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 "Beam_moments_by_element.dat" u 2:($15*100) w l not 
#
set ylabel "sigma_p [%]"
set origin 0.5,0.75
plot "Beam_moments_by_element.dat" u 2:($14**0.5*100) w l not