4 include 'COMMON.INTERACT'
5 include 'COMMON.IOUNITS'
9 read (*,*) nres,(itype(i),i=1,nres)
10 call seq2chains(nres,itype,nchain,chain_length,chain_border,
12 print *,"nres",nres," nchain",nchain
14 print *,"chain",i,chain_length(i),chain_border(1,i),
19 print *,i,ireschain(i)
21 call chain_symmetry(nchain,nres,itype,chain_border,
22 & chain_length,npermchain,tabpermchain)
23 print *,"ireschain permutations"
25 print '(60i4)',i,(tperm(ireschain(i),ii,tabpermchain),
28 print *,"residue permutations"
30 print '(60i4)',i,(iperm(i,ii),ii=1,npermchain)