1 subroutine seq2chains(nres,itype,nchain,chain_length,chain_border,
4 c Split the total UNRES sequence, which has dummy residues separating
5 c the chains, into separate chains. The length of chain ichain is
6 c contained in chain_length(ichain), the first and last non-dummy
7 c residues are in chain_border(1,ichain) and chain_border(2,ichain),
8 c respectively. The lengths pertain to non-dummy residues only.
12 integer nres,itype(nres),nchain,chain_length(nres),
13 & chain_border(2,nres),ireschain(nres)
18 chain_length(ichain)=0
21 if (itype(ii).eq.ntyp1) then
22 if (.not.new_chain) then
24 chain_border(2,ichain)=ii-1
26 chain_border(1,ichain)=ii+1
27 chain_length(ichain)=0
31 chain_border(1,ichain)=ii
34 chain_length(ichain)=chain_length(ichain)+1
38 if (itype(nres).eq.ntyp1) then
41 chain_length(ichain)=chain_length(ichain)+1
43 if (chain_length(ichain).gt.0) then
44 chain_border(2,ichain)=ii
51 do j=chain_border(1,i),chain_border(2,i)