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
       double precision IP,ISC(ntyp+1),mp,
     & msc(ntyp+1)
      common /inertia/ IP,ISC,MP,MSC