calc_delta_e

   1 #!/bin/bash
   2
   3 e0=`grep ^0 $1/energy | awk '{ print $4 }'`
   4 e1=`tail -n 1 $1/energy | awk '{ print $4 }'`
   5 ev=`grep ^19.000 $1/energy | awk '{ print $4 }'`
   6 en=`grep ^20.000 $1/energy | awk '{ print $4 }'`
   7
   8 echo "$e0 $e1 $2 $ev $en" | \
   9         awk '{ print $5*($3+1)-$4*$3-($2*($3+1)-$1*$3) " eV" }'