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

diff --git a/source/wham/src-HCD-5D/molread_zs.F b/source/wham/src-HCD-5D/molread_zs.F
index 878e4dd..e1a64f3 100644
--- a/source/wham/src-HCD-5D/molread_zs.F
+++ b/source/wham/src-HCD-5D/molread_zs.F
@@ -79,7 +79,7 @@ C Convert sequence to numeric code
         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
+      if (nchain.gt.1) chain_border1(1,nchain)=chain_border(1,nchain)-1
       chain_border1(2,nchain)=nres
       write(iout,*) "nres",nres," nchain",nchain
       do i=1,nchain
@@ -260,6 +260,7 @@ c       call flush(iout)
       endif
 
       call setup_var
+      write (iout,*) "Calling init_int_table"
       call init_int_table
       if (ns.gt.0) then
         write (iout,'(/a,i3,a)') 'The chain contains',ns,
@@ -332,6 +333,12 @@ C Read bridging residues.
       read (inp,*) ns,(iss(i),i=1,ns)
       print *,'ns=',ns
       write (iout,*) 'ns=',ns,' iss:',(iss(i),i=1,ns)
+c 5/24/2020 Adam: Added a table to translate residue numbers to cysteine
+c     numbers
+      icys=0
+      do i=1,ns
+        icys(iss(i))=i
+      enddo
 C Check whether the specified bridging residues are cystines.
       do i=1,ns
 	if (itype(iss(i)).ne.1) then