vasp_tools/e_fc

   1 #!/bin/bash
   2
   3 ready=`grep General\ timing $1/OUTCAR | awk '{ print $1 }'`
   4
   5 if [ "$ready" = "General" ]; then
   6         echo "$i seems to be finisched ..."
   7
   8 e=`./e_form plain_si_2333_nosym_sp plain_c_2333_nosym_sp $1 | grep M2 | \
   9         awk '{ print $4 }'`
  10 echo "E_f = $e"
  11
  12 ec100=3.15438
  13 r=`echo $e $ec100 $ec100 | awk '{ print $1-$2-$3 }'`
  14 echo "C 100: $r"
  15
  16 ec110=3.59738
  17 r=`echo $e $ec100 $ec110 | awk '{ print $1-$2-$3 }'`
  18 echo "C 110: $r"
  19
  20 ecsub=1.38891
  21 r=`echo $e $ec100 $ecsub | awk '{ print $1-$2-$3 }'`
  22 echo "C sub: $r"
  23
  24 ecbc=4.09638
  25 r=`echo $e $ec100 $ecbc | awk '{ print $1-$2-$3 }'`
  26 echo "C bc: $r"
  27
  28 esvac=3.62905
  29 r=`echo $e $ec100 $esvac | awk '{ print $1-$2-$3 }'`
  30 echo "Si vac: $r"
  31
  32 es110=3.38999
  33 r=`echo $e $ec100 $es110 | awk '{ print $1-$2-$3 }'`
  34 echo "Si 110: $r"
  35
  36 estet=3.76899
  37 r=`echo $e $ec100 $estet | awk '{ print $1-$2-$3 }'`
  38 echo "Si tet: $r"
  39
  40 eshex=3.41999
  41 r=`echo $e $ec100 $eshex | awk '{ print $1-$2-$3 }'`
  42 echo "Si hex: $r"
  43
  44 else
  45         echo "$i not yet finished ..."
  46 fi
  47