projects
/
physik
/
posic.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
Merge branch 'leadoff'
[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: "
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"]
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)
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
echo