cd print *,'nnt=',nnt,' nct=',nct
do j=nnt,i-1
itj=iabs(itype(j))
+ if (itj.eq.ntyp1) cycle
if (j.lt.i-1 .or. ipot.ne.4) then
rcomp=sigmaii(iti,itj)
else
fail=.true.
do while (fail.and.nsi.le.maxsi)
call gen_side(iti,theta(i+1),alph(i),omeg(i),fail)
+ call sc_coord_rebuild(i)
nsi=nsi+1
enddo
if(fail) goto 999
endif
enddo
- call chainbuild_extconf
+c write (iout,*) "before chaincuild overlap_sc_list: dc0",dc(:,0)
+c call chainbuild_extconf
+c write (iout,*) "after chaincuild overlap_sc_list: dc0",dc(:,0)
call overlap_sc_list(ioverlap,ioverlap_last)
write (iout,*) 'Overlaping residues ',ioverlap_last,
& (ioverlap(j),j=1,ioverlap_last)