econste_cmp.plt 343 Bytes
Newer Older
Francois Gygi committed
1 2 3
#!/bin/bash
# econste_cmp.plt: compare <econst> and <etotal> in two MD simulations
# use: econste_cmp.plt mdrun1.r mdrun2.r
4 5 6 7 8
if [ $1 == "-range" ]
then
  range=$2
  shift 2
fi
Francois Gygi committed
9
gnuplot -persist <<EOF
10
p $range "<grep -h econst $1" u 2 w l, "<grep -h '<etotal>' $1" u 2 w l, \
Francois Gygi committed
11 12 13
  "<grep -h econst $2" u 2 w l, "<grep -h '<etotal>' $2" u 2 w l

EOF