X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fposic.git;a=blobdiff_plain;f=vasp_tools%2Fcreate_intf_crude;fp=vasp_tools%2Fcreate_intf_crude;h=f5e0c3d66638bb32e145687e8c5b522d5c65a1c0;hp=0000000000000000000000000000000000000000;hb=6e6d7126ea9a845f11637d8e1b8eb2b570ac4dc9;hpb=97dc63eb6a519b8e1f4fbfaa9760dd94539436b0

diff --git a/vasp_tools/create_intf_crude b/vasp_tools/create_intf_crude
new file mode 100755
index 0000000..f5e0c3d
--- /dev/null
+++ b/vasp_tools/create_intf_crude
@@ -0,0 +1,69 @@
+#!/bin/bash
+
+si_file=$1
+sic_file=$2
+
+cnt_si_si=`sed -n 6p $si_file`
+cnt_si_sic=`sed -n 6p $sic_file | awk '{ print $1 }'`
+cnt_c_sic=`sed -n 6p $sic_file | awk '{ print $2 }'`
+
+lc_si=`sed -n 2p $si_file`
+
+x_si=`sed -n 3p $si_file | awk '{ print $1 }'`
+x_sic=`sed -n 3p $sic_file | awk '{ print $1 }'`
+
+h_si=`sed -n 5p $si_file | awk '{ print $3 }'`
+h_sic=`sed -n 5p $sic_file | awk '{ print $3 }'`
+h_sic=`echo $h_sic $x_si $x_sic | awk '{ print $1*$2/$3 }'`
+
+((tot_si=cnt_si_si+cnt_si_sic))
+
+h_tot=`echo $h_si + $h_sic | bc`
+
+((offset=8))
+
+# header:
+sed -n 1,4p $si_file
+echo " 0.00000 0.00000 $h_tot"
+echo " $tot_si $cnt_c_sic"
+sed -n 7,8p $si_file
+
+# silicon atoms of si lattice (bottom)
+((cnt=1))
+while [ $cnt -le $cnt_si_si ]; do
+	((gl=cnt+offset))
+	line=`sed -n ${gl}p $si_file`
+	x=`echo $line | awk '{ print $1 }'`
+	y=`echo $line | awk '{ print $2 }'`
+	z=`echo $line | awk '{ print $3 }'`
+	nz=`echo $z $h_tot $h_si | awk '{ print $1*$3/$2 }'`
+	echo " $x $y $nz T T T"
+	((cnt+=1))
+done
+
+# silicon atoms of sic lattice (top / 2nd fcc)
+((cnt=1))
+while [ $cnt -le $cnt_si_sic ]; do
+	((gl=cnt+offset+cnt_c_sic))
+	line=`sed -n ${gl}p $sic_file`
+	x=`echo $line | awk '{ print $1 }'`
+	y=`echo $line | awk '{ print $2 }'`
+	z=`echo $line | awk '{ print $3 }'`
+	nz=`echo $z $h_tot $h_si $h_sic | awk '{ print $1*$4/$2 + $3/$2 }'`
+	echo " $x $y $nz T T T"
+	((cnt+=1))
+done
+
+# carbon atoms of sic lattice (top / 1st fcc <- odd counts)
+((cnt=1))
+while [ $cnt -le $cnt_c_sic ]; do
+	((gl=cnt+offset))
+	line=`sed -n ${gl}p $sic_file`
+	x=`echo $line | awk '{ print $1 }'`
+	y=`echo $line | awk '{ print $2 }'`
+	z=`echo $line | awk '{ print $3 }'`
+	nz=`echo $z $h_tot $h_si $h_sic | awk '{ print $1*$4/$2 + $3/$2 }'`
+	echo " $x $y $nz T T T"
+	((cnt+=1))
+done
+