enddo
enddo
write(iout,*) "tube param"
- read(itube,*) epspeptube,sigmapeptube,tubetranenepep
+ read(itube,*) epspeptube,sigmapeptube,acavtubpep,bcavtubpep,
+ & ccavtubpep,dcavtubpep,tubetranenepep
sigmapeptube=sigmapeptube**6
sigeps=dsign(1.0D0,epspeptube)
epspeptube=dabs(epspeptube)