#!/bin/bash


gnuplot<<EOF
set output "Cu111_rVV10.ps"
set terminal postscript enhanced color
set multiplot layout 3,2
set tics font "Helvetica, 12"
set key font "Helvetica, 12"
set title font "Helvetica, 12"
set ylabel font "Helvetica, 12"
set xlabel font "Helvetica, 12"

set bmargin 2.0
set tmargin 2.0
set lmargin 3.5
set rmargin 3.5

set ylabel"reaction probability"
set key left
set title font "Helvetica, 12"
set title"pure H_2"
plot "pureH2" w lp lw 2 title"exp.", "MS.B86b-pureH2" w lp lw 2 title"MS-B86bl", "MS.B86b-rVV10-pureH2" w lp lw 2 title"MS-B86bl-rVV10", "MS.RPBE-rVV10-pureH2" w lp lw 2 title"MS-RPBEl-rVV10", "SRP48-pureH2" w lp lw 2 title"SRP48"
unset ylabel

set title font "Helvetica, 12"
set title"pure D_2"
plot "pureD2" w lp lw 2 title"", "MS.B86b-pureD2" w lp lw 2 title"", "MS.B86b-rVV10-pureD2" w lp lw 2 title"", "MS.RPBE-rVV10-pureD2" w lp lw 2 title"", "SRP48-pureD2" w lp lw 2 title""

set title font "Helvetica, 12"
set title"H_2 T_n = 1740K"
set ylabel"reaction probability"
set ylabel font "Helvetica, 12"
plot "Tn1740H2" w lp lw 2 title"", "MS.B86b-Tn1740H2" w lp lw 2 title"", "MS.B86b-rVV10-Tn1740H2" w lp lw 2 title"", "MS.RPBE-rVV10-Tn1740H2" w lp lw 2 title"", "SRP48-Tn1740H2" w lp lw 2 title""
unset ylabel

set title font "Helvetica, 12"
set title"H_2 T_n = 2000K"
plot "Tn2000H2" w lp lw 2 title"", "MS.B86b-Tn2000H2" w lp lw 2 title"", "MS.B86b-rVV10-Tn2000H2" w lp lw 2 title"", "MS.RPBE-rVV10-Tn2000H2" w lp lw 2 title"", "SRP48-Tn2000H2" w lp lw 2 title""

set title font "Helvetica, 12"
set title"D_2 T_n = 2100K" 
set ylabel font "Helvetica, 12"
set ylabel"reaction probability"
set xlabel font "Helvetica, 12"
set xlabel"average translational energy [kJ/mol]"
plot "Tn2100D2" w lp lw 2 title"", "MS.B86b-Tn2100D2" w lp lw 2 title"", "MS.B86b-rVV10-Tn2100D2" w lp lw 2 title"", "MS.RPBE-rVV10-Tn2100D2" w lp lw 2 title"", "SRP48-Tn2100D2" w lp lw 2 title""
unset ylabel

set title font "Helvetica, 12"
set title"H_2 Rendulic"
set xlabel font "Helvetica, 12"
plot "rendulic" w lp lw 2 title"", "MS.B86b-rendulic" w lp lw 2 title"", "MS.B86b-rVV10-rendulic" w lp lw 2 title"", "MS.RPBE-rVV10-rendulic" w lp lw 2 title"", "SRP48-rendulic" w lp lw 2 title""
EOF

ps2pdf Cu111_rVV10.ps
evince Cu111_rVV10.pdf