implicit none
      include 'DIMENSIONS'
      include 'COMMON.CHAIN'
      include 'COMMON.INTERACT'
      include 'COMMON.IOUNITS'
      integer i,ii,iii,ires
      integer tperm,iperm
      iout=6
      read (*,*) nres,(itype(i),i=1,nres)
      call seq2chains(nres,itype,nchain,chain_length,chain_border,
     &  ireschain)
      print *,"nres",nres," nchain",nchain
      do i=1,nchain
        print *,"chain",i,chain_length(i),chain_border(1,i),
     &    chain_border(2,i)
      enddo
      print *,"ireschain"
      do i=1,nres
        print *,i,ireschain(i)
      enddo
      call chain_symmetry(nchain,nres,itype,chain_border,
     &    chain_length,npermchain,tabpermchain)
      print *,"ireschain permutations"
      do i=1,nres
        print '(60i4)',i,(tperm(ireschain(i),ii,tabpermchain),
     &   ii=1,npermchain)
      enddo    
      print *,"residue permutations"
      do i=1,nres
        print '(60i4)',i,(iperm(i,ii),ii=1,npermchain)
      enddo    
      stop
      end