ncont=0
kkk=3
do i=nnt+kkk,nct
- iti=itype(i)
+ iti=iabs(itype(i))
do j=nnt,i-kkk
- itj=itype(j)
+ itj=iabs(itype(j))
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
- write (iout,*) "Hairpins:"
+ write (iout,*) "Hairpins:",nharp
do i=1,nharp
i1=iharp(1,i)
j1=iharp(2,i)