#!/bin/bash
gnuplot<<EOF
set output "vdw_sites.pdf"
set terminal pdf enhanced size 3in,3in
set xtics nomirror
set ytics nomirror
set ylabel "interaction energy [meV]" font "Helvetica, 12"
set xlabel "Z [\305]" font "Helvetica, 12"
set yrange [-50:50]
plot "Cu111_MS.PBE-rVV10_C3v_TOP_T90.vdw-well" w l lw 3 title"top site", "Cu111_MS.PBE-rVV10_C3v_BtH_T90.vdw-well" w l lw 3 title"bridge site", "Cu111_MS.PBE-rVV10_C3v_FCC_T90.vdw-well" w l lw 3 title"fcc site"
EOF