From d24930a9135ffd5927fbf65a71494f5cf772f264 Mon Sep 17 00:00:00 2001
From: hackbard <hackbard@sage.physik.uni-augsburg.de>
Date: Tue, 23 Mar 2010 17:39:25 +0100
Subject: [PATCH] mus out of mu Si and mu SiC (Tersoff way)

---
 calc_delta_e | 17 ++++++++++++++++-
 1 file changed, 16 insertions(+), 1 deletion(-)

diff --git a/calc_delta_e b/calc_delta_e
index 2b4130e..eff07c9 100755
--- a/calc_delta_e
+++ b/calc_delta_e
@@ -1,8 +1,21 @@
 #!/bin/bash
 
+# guess potential
+pot=`grep ^potential $1/config | awk '{ print $2 }'`
+
+echo
+if [ "$pot" = "tersoff" ]; then
+echo "potential: $pot"
+musi=-4.629595
+muc=-7.421824
+music=-12.37363
+else
+echo "potential: $pot"
 musi=-4.628414
 muc=-7.373091
 music=-12.679618
+fi
+echo
 
 file=`ls $1/atomic_conf_* | tail -1`
 atom_cnt=`grep '# \[P\]' $file | awk '{ print $3 }'`
@@ -21,5 +34,7 @@ echo "$si_cnt $c_cnt $e1 $e0" | \
 	awk '{ print "   Gao: "($3-$4)*($1+$2) }'
 echo "$si_cnt $c_cnt $musi $muc $e1 $music" | \
 	#awk '{ print "   Posselt: "$5*($1+$2)-0.5*($1+$2)*$6-0.5*($1-$2)*($3-$4) }'
-	awk '{ print "   Posselt: "$5*($1+$2)-$1*$3-$2*$4) }'
+	awk '{ print "   Posselt: "$5*($1+$2)-$1*$3-$2*$4 }'
+echo "$si_cnt $c_cnt $musi $music $e1 $music" | \
+	awk '{ print "   Tersoff: "$5*($1+$2)-$1*$3-$2*($4-$3) }'
 
-- 
2.39.2