source/maxlik/src_MD_T_maxlik-NEWCORR-PMF-5/testseqchains.f

   1       implicit none
   2       include 'DIMENSIONS'
   3       integer nres,itype(2000),nchain,chain_length(2000),
   4      & chain_border(2,100),npermchain,tabpermchain(maxchain,maxperm)
   5       integer i
   6       read (*,*) nres,(itype(i),i=1,nres)
   7       call seq2chains(nres,itype,nchain,chain_length,chain_border)
   8       print *,"nres",nres," nchain",nchain
   9       do i=1,nchain
  10         print *,"chain",i,chain_length(i),chain_border(1,i),
  11      &    chain_border(2,i)
  12       enddo
  13       call chain_symmetry(nchain,nres,itype,chain_border,
  14      &    chain_length,npermchain,tabpermchain)
  15       end