-! write (iout,*) "j",j," k",k
+ write (iout,*) "j",j," k",k,nres
+ j1=j
+ k1=k
+ if (j.gt.nres) j1=j-nres
+ if (k.gt.nres) k1=k-nres
+ mnumkk=molnum(k1)
+ mnumjj=molnum(j1)
+
+ if ((itype(k1,mnumkk).gt.ntyp_molec(mnumkk)).or.&
+ (itype(j1,mnumjj).gt.ntyp_molec(mnumjj))) cycle
+ write (iout,*) "j",j," k",k,itype(k1,mnumkk),itype(j1,mnumjj)