+! 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,cPoise,Rb,surfarea,
+ & reset_fricmat
+! Quantities used in Langevin dynamics calculations
double precision friction(3,0:MAXRES2),stochforc(3,0:MAXRES2),
& fricmat(MAXRES2,MAXRES2),fric_work(MAXRES6),
& stoch_work(MAXRES6),