X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-HCD%2Fmolread_zs.F;h=878e4ddc18c822583d497994f296929282eb758a;hb=df61eb1a0ad598d7c6c5b4037dee4bec200d5bc2;hp=d7f586da02e8668cf3b4eb080705094985eb66fa;hpb=5aac27f7c6317fa0975184363614b42716034716;p=unres.git

diff --git a/source/wham/src-HCD/molread_zs.F b/source/wham/src-HCD/molread_zs.F
index d7f586d..878e4dd 100644
--- a/source/wham/src-HCD/molread_zs.F
+++ b/source/wham/src-HCD/molread_zs.F
@@ -73,10 +73,18 @@ C Convert sequence to numeric code
       nct=nres
       call seq2chains(nres,itype,nchain,chain_length,chain_border,
      &  ireschain)
+      chain_border1(1,1)=1
+      chain_border1(2,1)=chain_border(2,1)+1
+      do i=2,nchain-1
+        chain_border1(1,i)=chain_border(1,i)-1
+        chain_border1(2,i)=chain_border(2,i)+1
+      enddo
+      chain_border1(1,nchain)=chain_border(1,nchain)-1
+      chain_border1(2,nchain)=nres
       write(iout,*) "nres",nres," nchain",nchain
       do i=1,nchain
         write(iout,*)"chain",i,chain_length(i),chain_border(1,i),
-     &    chain_border(2,i)
+     &    chain_border(2,i),chain_border1(1,i),chain_border1(2,i)
       enddo
       call chain_symmetry(nchain,nres,itype,chain_border,
      &    chain_length,npermchain,tabpermchain)