pc script
[physik/posic.git] / vasp_tools / pc_script
1 #!/bin/bash
2
3 if [ ! -d $1 ]; then
4         echo "not a valid directory ..."
5         exit
6 fi
7
8 if [ ! -f $1/OUTCAR ]; then
9         echo "no vasp data found in $1 ..."
10         exit
11 fi
12
13 if [ ! -d $1/video ]; then
14         echo "no video directory found ..."
15         echo "forgot to run outcar2moldyn?"
16         exit
17 fi
18
19 echo
20 echo "WARNING !!!!"
21 echo "-> reading in sx sy and sz is not correct ..."
22 echo "WARNING !!!!"
23 echo
24
25 line=`sed -n -e '/length of/,/FORCES/p' $1/OUTCAR | grep ^\ \ \  | tail -1`
26
27 sx=`echo $line | awk '{ print $1 }'`
28 sy=`echo $line | awk '{ print $2 }'`
29 sz=`echo $line | awk '{ print $3 }'`
30
31 file=`ls $1/video/atomic_conf*.xyz | tail -1`
32
33 echo "running on $file using $sx $sy $sz ..."
34
35 ./pc_calc $file $sx $sy $sz | grep ^pc | \
36         awk '{ print $2 " " $3 " " $4 " "$5 }' > pc.txt
37
38 echo "done"
39