X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc-HCD-5D%2Freadrtns.F;h=057f1ac9836f9a02c13944e3179f63707dd5b024;hb=1441ec2b7d68b4fd94014649e4bd9574db0bfc36;hp=a3229a6a437c1e9f67e0003095471e3325bcc917;hpb=79fd686e9d9aaa23e4fda802ee16379470e961d3;p=unres.git

diff --git a/source/cluster/wham/src-HCD-5D/readrtns.F b/source/cluster/wham/src-HCD-5D/readrtns.F
index a3229a6..057f1ac 100644
--- a/source/cluster/wham/src-HCD-5D/readrtns.F
+++ b/source/cluster/wham/src-HCD-5D/readrtns.F
@@ -242,8 +242,8 @@ C Read weights of the subsequent energy terms.
       do i=1,maxres
         dyn_ss_mask(i)=.false.
       enddo
-      do i=1,maxres-1
-        do j=i+1,maxres
+      do i=1,max_cyst-1
+        do j=i+1,max_cyst
           dyn_ssbond_ij(i,j)=1.0d300
         enddo
       enddo
@@ -716,6 +716,12 @@ C Read information about disulfide bridges.
 C Read bridging residues.
       read (inp,*) ns,(iss(i),i=1,ns)
 c      print *,'ns=',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