1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
#if HAVE_CONFIG_H
# include "config.fh"
#endif
subroutine force
#include "common.fh"
c
double precision tbeg,wraptime
integer i,j
integer inode
c
tbeg = wraptime()
c
c clear forces
c
do 100 j = 1, 3
do 200 i = 1, antot
ra(i,j,4) = 0.0d00
200 continue
100 continue
c
do 300 i = 1, antot
asev(i) = 0.0d00
300 continue
c
esvr = 0.0d00
esvd = 0.0d00
esvoh = 0.0d00
esva = 0.0d00
presf = 0.0d00
presfx = 0.0d00
presfy = 0.0d00
presfz = 0.0d00
c
call gather
c
c clear buffers
c
do i = 1, btot
xfrc(i) = 0.0d00
yfrc(i) = 0.0d00
zfrc(i) = 0.0d00
end do
c
c assign atoms to individual molecules
c
call pairs
call scatter
c
nrg(6) =esvr + esvd + esvoh + esva
nrg(13) = esvr
nrg(14) = esvd
nrg(17) = esvoh
nrg(21) = esva
nrg(10) = presf
nrg(18) = presfx
nrg(19) = presfy
nrg(20) = presfz
presf = presf / (xbox * ybox * zbox)
nrg(15) = presf/ 3.0d00
c
tmstat(2) = tmstat(2) + wraptime() - tbeg
c
return
end
|
