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_aa_tube=4.0d0*epspeptube*sigmapeptube**2
pep_bb_tube=-sigeps*4.0d0*epspeptube*sigmapeptube
- write(iout,*) pep_aa_tube,pep_bb_tube
+ 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)
sc_aa_tube_par(i)=4.0d0*epssctube*sigmasctube**2
sc_bb_tube_par(i)=-sigeps*4.0d0*epssctube*sigmasctube
- write(iout,*) sc_aa_tube_par(i), sc_bb_tube_par(i)
+ write(iout,*) sc_aa_tube_par(i), sc_bb_tube_par(i),tubetranene(i)
enddo
#ifdef OLDSCP