!rm sigma_xy_t.ps
set term postscript color enhanced
set output "sigma_xy_t.ps"
set multiplot
set size 1.,0.50
#set origin 0.,0.
#set ylabel "<x>[mm]"
#set xlabel "time [us]"
#set yrange [-30:30]
#set ytics -30,10,30
set xrange [0:tmax]
#plot "Beam_moments_by_element.dat" u ($5*1.e6):($3*1000) w l notitle
#
set origin 0.,0.
set ylabel "{/Symbol s}_x[mm]"
set yrange [0:30]
set ytics 0, 5,30
plot "Beam_moments_by_element.dat" u ($5*1.e6):($6**.5*1000) w l not, 'rfquad.dat'  u ($1*1.e6):($2*100+12) w l lc 3 t 'rf_h'
#plot 'rfquad.dat'  u ($1*1.e6):($2*100+12) w l lc 3 t 'rf_h'
#
#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 ($5*1.e6):($4*1000) w l not
#
set origin 0.0,0.5
set ylabel "{/Symbol s}_y[mm]"
set yrange [0:30]
set ytics 0,5,30
set xrange [0:tmax]
plot "Beam_moments_by_element.dat" u ($5*1.e6):($9**.5*1000) w l not, 'rfquad.dat' u ($1*1.e6):($3*500+14) w l lc 3 t 'rf_v'
#plot  'rfquad.dat' u ($1*1.e6):($3*100+14) w l lc 3 t 'rf_v'
#exit