set term postscript color enhanced
set output "sigma_xyz.ps"
set multiplot
set size 0.5,0.30
set origin 0.,0.
set ylabel "<x>[mm]"
set xlabel "turn"
set yrange [-20:20]
set ytics -20,10,20
set xrange [0:turns]
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.3
set size 0.5,0.3
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.3
set ylabel "{/Symbol s}_y[mm]"
set yrange [0:30]
set ytics 0,5,30
set xrange [0:turns]
plot "Beam_moments_by_element.dat" u 2:($9**.5*1000) w l not 
#exit
#
set origin 0.5,0.60
set ylabel "length[ns]"
set autoscale
set xrange [0:turns]
set ytics autofreq
plot "Beam_moments_by_element.dat" u 2:($12**.5*1e9) w l not 
#

set autoscale
set xrange [0.1:turns]
set xtics 0,turns/5
set origin 0.,0.6
unset ylabel
#set format y "%0.3t 10^{%L}"
set ylabel "<dp/p>[%]"
set ytics autofreq
plot "Beam_moments_by_element.dat" u 2:($15*100) w l not