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
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