6.5 Program 6.1, Fortran Code
95
tape9 *lastposi*tape11 *tape12 *tape13 *tape14 *tape15 *tape16 *
nstep
:0001000000
isave
:0000000001
iprint
:0000010000
itape
:0000000001
dt
:000.001250
===================================================
c program motor
c this program performs brownian dynamics with hydrodynamic
c interactions (Ermak and Mccammon, J. Chem. Phys. 69, 1352, 1982)
c for a harmonic dimer)
implicit none
common / block1 / rx, ry, rz
common / block2 / d, xic
common / block3 / fx, fy, fz
common / block4 / A, Fload, gamma, w, Fo
integer itape, step, i, nstep, isave, iprint, step1
integer n, n3
parameter ( n = 2, n3 = 3 * n )
real rx(n), ry(n), rz(n), fx(n), fy(n), fz(n)
real rx0(n), ry0(n), rz0(n)
real d(n3,n3), xic(n3)
real drxsq, drysq, drzsq,sigma,StDev
real pi, dt,ts
real consii, consij
real acv, acvsq
real avv, flv
real v, radius, vn, r
real acr,acrsq,avr, avrsq, flr
real rcx0, rcy0, rcz0
real rx210, ry210, rz210, r21sq0, r210
real rcx, rcy, rcz, rc2, rc, rsq
real rx21, ry21, rz21, r21sq, r21
real acdn,acdiff, acrotn, acrot, deff, peclet, cos1
real acd, vcx, vcy, vcz, sumvcx, sumvcy, sumvcz, vc, t
real sumr, sumr2,dlamda,a,w, Ein, t1, Fload, DPI
real rmean, rmeansq, deltar, ranf, dummy,work,efficiency
real Fo , sumFrat, kF, dtW, vcx2, sumvcx2,tau,mass
real gamma, STOKESEFF, ENERGCONVEFF, KINETICEFF, VCX2EQ
real CUBEROEFF, Frat, SUZUKIEFF
real mvcx, mvcx2, sumvc,mvcy, mvcz
real vcy2, vcz2, vcx1, t2, acdnx,DEFFX,acdiffx
parameter (pi = 3.1415927, DPI=pi*2)
character*8 fn(16), tape
tape = ’tape****’