pep_bb_tube=-sigeps*4.0d0*epspeptube*sigmapeptube
write(iout,*) pep_aa_tube,pep_bb_tube,tubetranenepep
do i=1,ntyp
- read(itube,*) epssctube,sigmasctube,tubetranene(i)
+ read(itube,*) epssctube,sigmasctube,acavtub(i),bcavtub(i),
+ & ccavtub(i),dcavtub(i),tubetranene(i)
sigmasctube=sigmasctube**6
sigeps=dsign(1.0D0,epssctube)
epssctube=dabs(epssctube)