it1=itype(i1,1)
it2=itype(i2,1)
if (mnum.eq.0) mnum=1
+ if ((it1.eq.0).or.(it2.eq.0)) then
+ it1=1
+ it2=1
+ endif
print *,"CONTACT",i1,i2,mnum,it1,it2
write (iout,'(i3,2x,a,i4,2x,a,i4,5f8.3,3f10.5)') &
i,restyp(it1,mnum),i1,restyp(it2,mnum),i2,cscore(i),&