implicit none include 'DIMENSIONS' integer nres,itype(2000),nchain,chain_length(2000), & chain_border(2,100),npermchain,tabpermchain(maxchain,maxperm) integer i read (*,*) nres,(itype(i),i=1,nres) call seq2chains(nres,itype,nchain,chain_length,chain_border) print *,"nres",nres," nchain",nchain do i=1,nchain print *,"chain",i,chain_length(i),chain_border(1,i), & chain_border(2,i) enddo call chain_symmetry(nchain,nres,itype,chain_border, & chain_length,npermchain,tabpermchain) end