added acos to angle calc

[physik/posic.git] / vasp_tools / angle_calc

diff --git a/vasp_tools/angle_calc b/vasp_tools/angle_calc
index da01dfb..b23d51d 100755 (executable)
--- a/vasp_tools/angle_calc
+++ b/vasp_tools/angle_calc
@@ -41,7 +41,7 @@ yc=`echo $temp | awk '{ print $2 }'`
 zc=`echo $temp | awk '{ print $3 }'`
 
 echo -en "angle: "
-echo "$xa $ya $za $xb $yb $zb $xc $yc $zc $scale" | \
+foo=`echo "$xa $ya $za $xb $yb $zb $xc $yc $zc $scale" | \
        awk ' \
        BEGIN {
                X1=ENVIRON["X1"]; X2=ENVIRON["X2"]; X3=ENVIRON["X3"]
@@ -92,8 +92,10 @@ echo "$xa $ya $za $xb $yb $zb $xc $yc $zc $scale" | \
                sp=dxt*Dxt+dyt*Dyt+dzt*Dzt
                d=sqrt(dxt^2+dyt^2+dzt^2)
                D=sqrt(Dxt^2+Dyt^2+Dzt^2)
-               print 180.0/pi*acos2(sp/(d*D))
-       }'
+               print sp/(d*D)
+       }'`
+
+./acos $foo
 
 echo