#!/bin/bash
# generate msd.png with average curve 
#
title="MSD All + avg"
# assemble gnuplot plot command
cmd="plot \"all_msd.dat\" u 0:1 w l lt 1"
for n in {1..16}
do 
  cmd=$cmd",\"\" u 0:$n w l lt 1"
done
cmd=$cmd", \"msdavg.dat\" u 0:1 lc rgb 'black'"
echo $cmd

gnuplot -persist << EOF
#set term svg size 640,480 fsize 18 linewidth 2
set term svg dynamic enhanced size 640,480 fsize 16 
set output "msd_with_avg.svg"
#set title "E_{KS} (avg)"
#set data style lines
set key off
set grid
set ylabel 'MSD (Bohr^2)'
set xtics ("0" 0, "10" 40000, "20" 80000, "30" 120000, "40" 160000)
set xlabel 'time (ps)'
$cmd

set term png large font arial 18
set output "msd_with_avg.png"
$cmd

set term postscript enhanced color font "Helvetica, 18"
set output "msd_with_avg.eps"
$cmd
quit
EOF