enddo
enddo
write(iout,*) "tube param"
- read(itube,*) epspeptube,sigmapeptube
+ read(itube,*) epspeptube,sigmapeptube,tubetranenepep
sigmapeptube=sigmapeptube**6
sigeps=dsign(1.0D0,epspeptube)
epspeptube=dabs(epspeptube)
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
+ read(itube,*) epssctube,sigmasctube,tubetranene(i)
sigmasctube=sigmasctube**6
sigeps=dsign(1.0D0,epssctube)
epssctube=dabs(epssctube)