C Torsional constants of the rotation about virtual-bond dihedral angles
      double precision v1,v2,vlor1,vlor2,vlor3,v0
      integer itortyp,ntortyp,nterm,nlor,nterm_old
      common/torsion/v0(maxtor,maxtor),v1(maxterm,maxtor,maxtor),
     &    v2(maxterm,maxtor,maxtor),vlor1(maxlor,maxtor,maxtor),
     &    vlor2(maxlor,maxtor,maxtor),vlor3(maxlor,maxtor,maxtor),
     &    itortyp(ntyp),ntortyp,nterm(maxtor,maxtor),
     &    nlor(maxtor,maxtor),nterm_old
C 6/23/01 - constants for double torsionals
      double precision v1c,v1s,v2c,v2s
      integer ntermd_1,ntermd_2
      common /torsiond/ v1c(2,maxtermd_1,maxtor,maxtor,maxtor),
     &    v1s(2,maxtermd_1,maxtor,maxtor,maxtor),
     &    v2c(maxtermd_2,maxtermd_2,maxtor,maxtor,maxtor),
     &    v2s(maxtermd_2,maxtermd_2,maxtor,maxtor,maxtor),
     &    ntermd_1(maxtor,maxtor,maxtor),ntermd_2(maxtor,maxtor,maxtor)
C 9/18/99 - added Fourier coeffficients of the expansion of local energy 
C           surface
      double precision b1,b2,cc,dd,ee,ctilde,dtilde,b1tilde
      integer nloctyp
      common/fourier/ b1(2,maxtor),b2(2,maxtor),cc(2,2,maxtor),
     &    dd(2,2,maxtor),ee(2,2,maxtor),ctilde(2,2,maxtor),
     &    dtilde(2,2,maxtor),b1tilde(2,maxtor),nloctyp
      double precision b
      common /fourier1/ b(13,maxtor)