if (acavtub(iti).eq.0.0d0) go to 667
denominator=(1.0+dcavtub(iti)*rdiff6*rdiff6)
enecavtube(i)=
- & acavtub(iti)*(rdiff+bcavtub(iti)*sqrt(rdiff)+cavtub(iti))
+ & acavtub(iti)*(rdiff+bcavtub(iti)*sqrt(rdiff)+ccavtub(iti))
& /denominator
faccav=(acavtub(iti)*(1.0+bcavtub(iti)/2.0/sqrt(rdiff))
& *denominator-acavtub(iti)*(rdiff+bcavtub(iti)*sqrt(rdiff)
- & +cavtub(iti))*rdiff6**2.0d0/rdiff*dcavtub(iti))
+ & +ccavtub(iti))*rdiff6**2.0d0/rdiff*dcavtub(iti))
& /denominator**2.0d0
fac=fac+faccav
667 continue