- double precision friction(3,0:MAXRES2),stochforc(3,0:MAXRES2),
+! Basic Langevin dynamics parameters
+ logical surfarea
+ integer reset_fricmat
+ double precision scal_fric,rwat,etawat,gamp,
+ & gamsc(ntyp1),stdfp,stdfsc(ntyp),stdforcp(MAXRES),
+ & stdforcsc(MAXRES),pstok,restok(ntyp+1),cPoise,Rb
+ common /langevin/ pstok,restok,gamp,gamsc,
+ & stdfp,stdfsc,stdforcp,stdforcsc,rwat,etawat,scal_fric,
+ & cPoise,Rb,surfarea,reset_fricmat
+! Variables used in Langevin dynamics calculations
+ double precision friction(3,0:MAXRES2),stochforc(3,0:MAXRES2),
& fricmat(MAXRES2,MAXRES2),fric_work(MAXRES6),
- & stoch_work(MAXRES6),
- & fricgam(MAXRES6),fricvec(MAXRES2,MAXRES2)
- logical flag_stoch(0:maxflag_stoch)
- common /langforc/ friction,stochforc,
- & fricmat,fric_work,fricgam,stoch_work,fricvec,vrand_mat1,
- & vrand_mat2,prand_mat,vfric_mat,afric_mat,pfric_mat,
- & pfric0_mat,afric0_mat,vfric0_mat,prand0_mat,vrand0_mat1,
- & vrand0_mat2,flag_stoch
- common /langmat/ mt1,mt2,mt3
+ & stoch_work(MAXRES6),fricgam(MAXRES6),fricvec(MAXRES2,MAXRES2)
+ logical flag_stoch(0:maxflag_stoch)
+ common /langforc/ friction,stochforc,fricmat,fric_work,fricgam,
+ & stoch_work,fricvec,flag_stoch