if (ipot.ne.4) then
c rcomp=sigmaii(iti,itj)+1.0D0
rcomp=facont*sigmaii(iti,itj)
if (ipot.ne.4) then
c rcomp=sigmaii(iti,itj)+1.0D0
rcomp=facont*sigmaii(iti,itj)
c write (iout,*)'i',i,' iharp',(iharp(k,i),k=1,4)
c enddo
if (lprint) then
c write (iout,*)'i',i,' iharp',(iharp(k,i),k=1,4)
c enddo
if (lprint) then