integer :: rescale_mode
real(kind=8) :: wsc,wscp,welec,wbond,wstrain,wtor,wtor_d,wang,&
wscloc,wcorr,wcorr4,wcorr5,wcorr6,wsccor,wel_loc,wturn3,wturn4,&
- wturn6,wvdwpp,wliptran,wshield,lipscale
+ wturn6,wvdwpp,wliptran,wshield,lipscale,wtube
#ifdef CLUSTER
real(kind=8) :: scalscp
#endif
! common.interact
! common /interact/
real(kind=8),dimension(:,:),allocatable :: aa_aq,bb_aq,augm,aa_lip,bb_lip !(ntyp,ntyp)
+ real(kind=8),dimension(:),allocatable :: sc_aa_tube_par,sc_bb_tube_par,&
+ acavtub,bcavtub,ccavtub,dcavtub,tubetranene
+ real(kind=8) :: acavtubpep,bcavtubpep,ccavtubpep,dcavtubpep, &
+ tubetranenepep,pep_aa_tube,pep_bb_tube,tubeR0
+ real(kind=8),dimension(3) :: tubecenter
real(kind=8),dimension(:,:),allocatable :: aad,bad !(ntyp,2)
real(kind=8),dimension(2,2) :: app,bpp,ael6,ael3
integer :: expon,expon2, nnt,nct,itypro
ibondp_start,ibondp_end,ivec_start,ivec_end,iset_start,iset_end,&
iturn3_start,iturn3_end,iturn4_start,iturn4_end,iint_start,&
iint_end,iphi1_start,iphi1_end,itau_start,itau_end,&
- ilip_start,ilip_end
+ ilip_start,ilip_end,itube_start,itube_end
integer,dimension(:),allocatable :: ibond_displ,ibond_count,&
ithet_displ,ithet_count,iphi_displ,iphi_count,iphi1_displ,&
iphi1_count,ivec_displ,ivec_count,iset_displ,iset_count,&