!rm sigma_xy.ps
set term postscript color enhanced
set output "sigma_xy.ps"
set multiplot
set size 0.5,0.50
set origin 0.,0.
set ylabel "<x>[mm]"
set xlabel "turn"
set yrange [-30:30]
set ytics -30,10,30
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.5
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.5
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