double precision aa,bb,augm,aad,bad,app,bpp,ale6,ael3,ael6,
     &aa_lip,bb_lip,aa_aq,bb_aq,sc_aa_tube_par,sc_bb_tube_par,
     & pep_aa_tube,pep_bb_tube
      double precision wdti,wdti2,wdti4,wdti8,
     &                 wdtii,wdtii2,wdtii4,wdtii8 
      common /nosehoover_dt/
     &   wdti(maxyosh),wdti2(maxyosh),wdti4(maxyosh),wdti8(maxyosh),
     &   wdtii(maxyosh),wdtii2(maxyosh),wdtii4(maxyosh),wdtii8(maxyosh)
      integer expon,expon2
      integer nnt,nct,nint_gr,istart,iend,itype,itel,itypro,
     & ielstart,ielend,ielstart_vdw,ielend_vdw,nscp_gr,iscpstart,
     & iscpend,iatsc_s,iatsc_e,
     & iatel_s,iatel_e,iatscp_s,iatscp_e,iatel_s_vdw,iatel_e_vdw,
     & ispp,iscp
      common /interact/aa_aq(ntyp,ntyp),bb_aq(ntyp,ntyp),
     & aa_lip(ntyp,ntyp),bb_lip(ntyp,ntyp),
     & sc_aa_tube_par(ntyp),sc_bb_tube_par(ntyp),
     & pep_aa_tube,pep_bb_tube,
     & augm(ntyp,ntyp),
     & aad(ntyp,2),bad(ntyp,2),app(2,2),bpp(2,2),ael6(2,2),ael3(2,2),
     & expon,expon2,nnt,nct,nint_gr(maxres),istart(maxres,maxint_gr),
     & iend(maxres,maxint_gr),itype(maxres),itel(maxres),itypro,
     & ielstart(maxres),ielend(maxres),ielstart_vdw(maxres),
     & ielend_vdw(maxres),nscp_gr(maxres),
     & iscpstart(maxres,maxint_gr),iscpend(maxres,maxint_gr),
     & iatsc_s,iatsc_e,iatel_s,iatel_e,iatel_s_vdw,iatel_e_vdw,
     & iatscp_s,iatscp_e,ispp,iscp
C 12/1/95 Array EPS included in the COMMON block.
      double precision eps,epslip,sigma,sigmaii,rs0,chi,chip,alp,
     & sigma0,sigii,
     & rr0,r0,r0e,r0d,rpp,epp,elpp6,elpp3,eps_scp,rscp
      common /body/eps(ntyp,ntyp),epslip(ntyp,ntyp),
     & sigma(0:ntyp1,0:ntyp1),
     & sigmaii(ntyp,ntyp),
     & rs0(ntyp,ntyp),chi(ntyp,ntyp),chip(ntyp),alp(ntyp),sigma0(ntyp),
     & sigii(ntyp),rr0(ntyp),r0(ntyp,ntyp),r0e(ntyp,ntyp),r0d(ntyp,2),
     & rpp(2,2),epp(2,2),elpp6(2,2),elpp3(2,2),eps_scp(ntyp,2),
     & rscp(ntyp,2)
c 12/5/03 modified 09/18/03 Bond stretching parameters.
      double precision vbldp0,vbldpDUM,
     &  vbldsc0,akp,aksc,abond0,distchainmax
      integer nbondterm
      common /stretch/ vbldp0,vbldpDUM,
     & vbldsc0(maxbondterm,ntyp),akp,
     & aksc(maxbondterm,ntyp),abond0(maxbondterm,ntyp),
     & distchainmax,nbondterm(ntyp)
C 01/29/15 Lipidic parameters
      double precision   pepliptran,liptranene,
     &tubetranene, tubetranenepep
      common /lipid/ pepliptran,liptranene(ntyp)
      common /tubepar/ tubetranene(ntyp), tubetranenepep