3 integer nres,itype(2000),nchain,chain_length(2000),
4 & chain_border(2,100),npermchain,tabpermchain(maxchain,maxperm)
6 read (*,*) nres,(itype(i),i=1,nres)
7 call seq2chains(nres,itype,nchain,chain_length,chain_border)
8 print *,"nres",nres," nchain",nchain
10 print *,"chain",i,chain_length(i),chain_border(1,i),
13 call chain_symmetry(nchain,nres,itype,chain_border,
14 & chain_length,npermchain,tabpermchain)