104
6 Ratchet Dimer Brownian Motor with Hydrodynamic Interactions
* ( oij - 3.0 * rpij ) * ryij * rzij * rrijsq
d( ic+1, jc ) = d( ic , jc+1 )
d( ic+2, jc ) = d( ic , jc+2 )
d( ic+2, jc+1 ) = d( ic+1, jc+2 )
c ** calculate systematic forces **
delta = 0.1
deltasc = delta**2
vij = 0.5*(1./deltasc)*(1.-rij)**2
wij = (1./deltasc)*(1.-rij)*rij
fij = wij * rrijsq
fxij = fij * rxij
fyij = fij * ryij
fzij = fij * rzij
v = v + vij
fxi = fxi + fxij
fyi = fyi + fyij
fzi = fzi + fzij
fx(j) = fx(j) - fxij
fy(j) = fy(j) - fyij
fz(j) = fz(j) - fzij
write(9,*) 1. - RIJ, VIJ
99 continue
fx(i) = fxi
fy(i) = fyi
fz(i) = fzi
100 continue
DPI=6.2832
FAC = A*abs(sin(w*t))
DO I = 1,2
FX(I)= FX(I)+ Fo*(cos(DPI*RX(I)) - 0.5*cos(2*DPI*RX(I))) +
* Fload + FAC
END DO
do 50 i = 1, n
c ** calculate on-diagonal blocks of diffusion tensor **
ic = 3 * ( i - 1 ) + 1
d( ic , ic ) = consii
d( ic+1, ic+1 ) = consii
d( ic+2, ic+2 ) = consii
d( ic , ic+1 ) = 0.0
d( ic , ic+2 ) = 0.0
d( ic+1, ic+2 ) = 0.0
50 continue
c ** fill the lower triangle of the diffusion tensor **
do 70 ic = 1, n3 - 1
do 60 jc = ic + 1, n3