1 c--------------------------------------------------------------------------
3 c Calculating the adjusted velocities for accelerations
4 implicit real*8 (a-h,o-z)
6 include 'COMMON.CONTROL'
9 include 'COMMON.LANGEVIN'
10 include 'COMMON.CHAIN'
11 include 'COMMON.DERIV'
13 include 'COMMON.LOCAL'
14 include 'COMMON.INTERACT'
15 include 'COMMON.IOUNITS'
16 include 'COMMON.NAMES'
17 double precision stochforcvec(MAXRES6),stochforcvecV(MAXRES6)
18 common /stochcalc/ stochforcvec
20 c Compute the stochastic forces which contribute to velocity change
22 call stochastic_force(stochforcvecV)
28 ddt1=ddt1+vfric_mat(i,j)*d_a_work(j)
29 ddt2=ddt2+vrand_mat1(i,j)*stochforcvec(j)+
30 & vrand_mat2(i,j)*stochforcvecV(j)
32 d_t_work(i)=d_t_work_new(i)+0.5d0*ddt1+ddt2
40 d_t(j,i)=d_t_work(ind+j)
45 if (itype(i).ne.10) then
48 d_t(j,inres)=d_t_work(ind+j)