#!/bin/bash
# generate eavg.svg plot 
#
# assemble gnuplot plot command

gnuplot -persist << EOF

set term postscript enhanced eps color font "Helvetica, 18"
set output "eavg.eps"

#set title "E_{KS} (avg)" font "Helvetica,24"

set key off
set grid
set xrange [0:43.53993]
set ylabel 'E_{KS}+1100 (a.u.)' font "Helvetica,22"
set xlabel 'time (ps)' font "Helvetica,22"
set ytics format "%8.2f"
set tics font "Helvetica,18"
# ps = duration of one run in ps
# one run is 2000 steps of 10 au_t
# 20000 au_t = 0.483777 ps
fac=0.483777
e0=1100
plot "eagnu.dat" i 0 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 1 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 2 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 3 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 4 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 5 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 6 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 7 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 8 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 9 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 10 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 11 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 12 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 13 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 14 u (\$0*fac):(\$1+e0) w l lt 1 lw 4, \
              "" i 15 u (\$0*fac):(\$1+e0) w l lt 1 lw 4 

quit
EOF