calc_delta_e

   1 #!/bin/bash
   2
   3 file=`ls $1/atomic_conf_* | tail -1`
   4 atom_cnt=`grep '# \[P\]' $file | awk '{ print $3 }'`
   5 e0=`awk 'NR==2' $1/energy | awk '{ print $4 }'`
   6 e1=`tail -n 1 $1/energy | awk '{ print $4 }'`
   7
   8 echo "$e0 $e1 $atom_cnt" | \
   9         awk '{ print ($2-$1)*$3" eV" }'