102
6 Ratchet Dimer Brownian Motor with Hydrodynamic Interactions
integer n
parameter ( n = 2 )
character cnfile*(*)
integer itape
real rx(n), ry(n), rz(n)
c ***************************************************************
open(unit=itape,file=cnfile,form=’formatted’,status=’unknown’)
rewind(itape)
do i = 1,n
end do
return
end
subroutine force ( consii, consij, v, r, t)
common / block1 / rx, ry, rz
common / block2 / d, xic
common / block3 / fx, fy, fz
common / block4 / A, Fload, gamma, w, Fo
c ***************************************************************
c ** routine to compute systematic forces and the diffusion tensor **
c ** **
c ** principal variables: **
c ** **
c ** integer n number of atoms **
c ** integer n3 number of degrees of freedom **
c ** real rx(n),ry(n),rz(n) positions **
c ** real fx(n),fy(n),fz(n) forces **
c ** real d(n3,n3) the diffusion tensor **
c ** real xic(n3) correlated random normal deviates **
c ** real consii constant in the diffusion tensor **
c ** real consij constant in the diffusion tensor **
c ** real v the potential energy **
c ** **
c ** usage: **
c ** **
c ** force is called in a brownian dynamics program to calculate **
c ** the systematic force on each atom and the elements of the **
c ** diffusion tensor. **
c ***************************************************************
integer n, n3
parameter ( n = 2, n3 = n * 3 )
real consii, consij, v, DPI
real rx(n), ry(n), rz(n), fx(n), fy(n), fz(n)
real d(n3,n3), xic(n3)
integer ic, jc, i, j
real rxi, ryi, rzi, fxij, fyij, fzij, fij