+ else if (molnum(j).eq.2) then
+ itj=iabs(itype(j,2))
+ if (dist(nres+i,nres+j).lt.redfac*sigma_nucl(iti,itj)) then
+ overlap=.true.
+
+! print *,'overlap, SC-SC: i=',i,' j=',j,
+! & ' dist=',dist(nres+i,nres+j),' rcomp=',
+! & rcomp
+ return
+ endif
+
+ endif