dodanie czytania wagi wsccor do clustrowania-M
[unres.git] / source / cluster / unres / src / intcor.f
1 C
2 C------------------------------------------------------------------------------
3 C
4       double precision function alpha(i1,i2,i3)
5 c
6 c  Calculates the planar angle between atoms (i1), (i2), and (i3).
7 c
8       include 'DIMENSIONS'
9       include 'COMMON.GEO'
10       include 'COMMON.CHAIN'
11       if (i1.lt.1 .or. i2.lt.1 .or. i3.lt.1 .or. i4.lt.1) then
12         alpha=0.0
13         return
14       endif
15       x12=c(1,i1)-c(1,i2)
16       x23=c(1,i3)-c(1,i2)
17       y12=c(2,i1)-c(2,i2)
18       y23=c(2,i3)-c(2,i2)
19       z12=c(3,i1)-c(3,i2)
20       z23=c(3,i3)-c(3,i2)
21       vnorm=dsqrt(x12*x12+y12*y12+z12*z12)
22       wnorm=dsqrt(x23*x23+y23*y23+z23*z23)
23       scalar=(x12*x23+y12*y23+z12*z23)/(vnorm*wnorm)
24       alpha=arcos(scalar)
25       return
26       end
27 C
28 C------------------------------------------------------------------------------
29 C
30       double precision function beta(i1,i2,i3,i4)
31 c
32 c  Calculates the dihedral angle between atoms (i1), (i2), (i3) and (i4)
33 c
34       include 'DIMENSIONS'
35       include 'COMMON.GEO'
36       include 'COMMON.CHAIN'
37       if (i1.lt.1 .or. i2.lt.1 .or. i3.lt.1 .or. i4.lt.1) then
38         beta=0.0
39         return
40       endif
41       x12=c(1,i1)-c(1,i2)
42       x23=c(1,i3)-c(1,i2)
43       x34=c(1,i4)-c(1,i3)
44       y12=c(2,i1)-c(2,i2)
45       y23=c(2,i3)-c(2,i2)
46       y34=c(2,i4)-c(2,i3)
47       z12=c(3,i1)-c(3,i2)
48       z23=c(3,i3)-c(3,i2)
49       z34=c(3,i4)-c(3,i3)
50 cd    print '(2i3,3f10.5)',i1,i2,x12,y12,z12
51 cd    print '(2i3,3f10.5)',i2,i3,x23,y23,z23
52 cd    print '(2i3,3f10.5)',i3,i4,x34,y34,z34
53       wx=-y23*z34+y34*z23
54       wy=x23*z34-z23*x34
55       wz=-x23*y34+y23*x34
56       wnorm=dsqrt(wx*wx+wy*wy+wz*wz)
57       vx=y12*z23-z12*y23
58       vy=-x12*z23+z12*x23
59       vz=x12*y23-y12*x23
60       vnorm=dsqrt(vx*vx+vy*vy+vz*vz)
61       if (vnorm.gt.1.0D-13 .and. wnorm.gt.1.0D-13) then
62       scalar=(vx*wx+vy*wy+vz*wz)/(vnorm*wnorm)
63       if (dabs(scalar).gt.1.0D0) 
64      &scalar=0.99999999999999D0*scalar/dabs(scalar)
65       angle=dacos(scalar)
66 cd    print '(2i4,10f7.3)',i2,i3,vx,vy,vz,wx,wy,wz,vnorm,wnorm,
67 cd   &scalar,angle
68       else
69       angle=pi
70       endif 
71 c     if (angle.le.0.0D0) angle=pi+angle
72       tx=vy*wz-vz*wy
73       ty=-vx*wz+vz*wx
74       tz=vx*wy-vy*wx
75       scalar=tx*x23+ty*y23+tz*z23
76       if (scalar.lt.0.0D0) angle=-angle
77       beta=angle
78       return
79       end
80 C
81 C------------------------------------------------------------------------------
82 C
83       function dist(i1,i2)
84 c
85 c  Calculates the distance between atoms (i1) and (i2).
86 c
87       include 'DIMENSIONS'
88       include 'COMMON.GEO'
89       include 'COMMON.CHAIN'
90       if (i1.lt.1 .or. i2.lt.1) then
91         dist=0.0
92         return
93       endif
94       x12=c(1,i1)-c(1,i2)
95       y12=c(2,i1)-c(2,i2)
96       z12=c(3,i1)-c(3,i2)
97       dist=dsqrt(x12*x12+y12*y12+z12*z12)
98       return
99       end
100 C