#!/bin/bash


gnuplot<<EOF

set output "test.ps"
set terminal postscript enhanced color landscape
set multiplot layout 1,4 
#title "E_0({/Symbol n},J) Au(111)" font ",12"
#set bmargin 0.5
#set tmargin 2.5
#set lmargin 1.5
#set rmargin 1.5
set tics font "Helvetica, 10"
#
set title "H_2 {/Symbol n} = 0" font ",10" 
unset key
set ylabel "E_0({/Symbol n},J) [ev]
set xlabel "J"
set yrange[0.6:1.35]
set xrange[-0.5:9.5]
plot "exp.H2.v0" lc -1 pointtype 7 pointsize 0.5 , "B86SRP68-DF2.H2_v0_10pr" lc 1 pointtype 7 pointsize 0.5, "optPBE-DF2.H2_v0_10pr" lc 2 pointtype 7 pointsize 0.5, "SRP48.H2_v0_10pr" lc 3 pointtype 7 pointsize 0.5, "MS.PBE.H2_v0_10pr" lc 4 pointtype 7 pointsize 0.5, "PBEa57-DF2.H2_v0_10pr" lc 5 pointtype 7 pointsize 0.5, "PBE.H2_v0_10pr" lc 8 pointtype 8 pointsize 0.5
unset ylabel
#
set title "H_2 {/Symbol n} = 1"
unset key
set xlabel "J"
set yrange[0.6:1.35]
set xrange[-0.5:7.5]
set ytics format " "
plot "exp.H2.v1" lc -1 pointtype 7 pointsize 0.5 , "B86SRP68-DF2.H2_v1_10pr" lc 1 pointtype 7 pointsize 0.5, "optPBE-DF2.H2_v1_10pr" lc 2 pointtype 7 pointsize 0.5, "SRP48.H2_v1_10pr" lc 3 pointtype 7 pointsize 0.5, "MS.PBE.H2_v1_10pr" lc 4 pointtype 7 pointsize 0.5, "PBEa57-DF2.H2_v1_10pr" lc 5 pointtype 7 pointsize 0.5, "PBE.H2_v1_10pr" lc 8 pointtype 8 pointsize 0.5

#
set title"D_2 {/Symbol n} = 0"
unset key
set xlabel"J"
set yrange[0.6:1.35]
set xrange [-0.5:12.5]
plot "exp.D2.v0" lc -1 pointtype 7 pointsize 0.5 , "B86SRP68-DF2.D2_v0_10pr" lc 1 pointtype 7 pointsize 0.5, "optPBE-DF2.D2_v0_10pr" lc 2 pointtype 7 pointsize 0.5, "SRP48.D2_v0_10pr" lc 3 pointtype 7 pointsize 0.5, "MS.PBE.D2_v0_10pr" lc 4 pointtype 7 pointsize 0.5, "PBEa57-DF2.D2_v0_10pr" lc 5 pointtype 7 pointsize 0.5, "PBE.D2_v0_10pr" lc 8 pointtype 8 pointsize 0.5

#
set title"D_2 {/Symbol n} = 1"
set key 
#set key left 
set key font "Helvetica, 7"
set xlabel "J"
set yrange[0.6:1.35]
#set yrange[0.88:1.05]
set xrange [-0.5:9.5]
plot "exp.D2.v1" title "exp." lc -1 pointtype 7 pointsize 0.5 , "B86SRP68-DF2.D2_v1_10pr" title "B86SRP68-DF2" lc 1 pointtype 7 pointsize 0.5 , "optPBE-DF2.D2_v1_10pr" lc 2 pointtype 7 pointsize 0.5 title "optPBE-DF2", "SRP48.D2_v1_10pr" lc 3 pointtype 7 pointsize 0.5 title "SRP48", "MS.PBE.D2_v1_10pr" lc 4 pointtype 7 pointsize 0.5 title"MS.PBE", "PBEa57-DF2.D2_v1_10pr" lc 5 pointtype 7 pointsize 0.5 title"PBEa57-DF2", "PBE.D2_v1_10pr" lc 8 pointtype 8 pointsize 0.5 title"PBE"

EOF
ps2pdf test.ps
evince test.pdf