]> www.hackdaworld.org Git - physik/posic.git/blob - vasp_tools/e_fc
Merge branch 'leadoff'
[physik/posic.git] / 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