X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MIN%2Freadrtns_min.F;h=5c505655e1ed470e2b56fa1c3b9a9ea7457118b0;hb=refs%2Fheads%2Fhomology;hp=dee929a87dbe8586dcee9d7504dc16702e5e118c;hpb=8893a1b93e875425e4b7ceb78242d488665af7b8;p=unres.git

diff --git a/source/unres/src_MIN/readrtns_min.F b/source/unres/src_MIN/readrtns_min.F
index dee929a..5c50565 100644
--- a/source/unres/src_MIN/readrtns_min.F
+++ b/source/unres/src_MIN/readrtns_min.F
@@ -77,6 +77,7 @@ C
       character*8 diagmeth(0:3) /'Library','EVVRSP','Givens','Jacobi'/
       character*80 ucase
       character*320 controlcard
+      external random_init
 
       nglob_csa=0
       eglob_csa=1d99
@@ -422,6 +423,9 @@ c        print *,'Finished reading pdb data'
 c        call contact(.false.,ncont_ref,icont_ref,co)
 
         if (sideadd) then 
+C Following 2 lines for diagnostics; comment out if not needed
+         write (iout,*) "Before sideadd"
+         call intout
          if(me.eq.king.or..not.out1file)
      &    write(iout,*)'Adding sidechains'
          maxsi=1000
@@ -438,7 +442,12 @@ c        call contact(.false.,ncont_ref,icont_ref,co)
      &              i,' after ',nsi,' trials'
           endif
          enddo
+C 9/29/12 Adam: Recalculate coordinates with new side chain positions
+         call chainbuild
         endif  
+C Following 2 lines for diagnostics; comment out if not needed
+        write (iout,*) "After sideadd"
+        call intout
       endif
       if (indpdb.eq.0) then
 C Read sequence if not taken from the pdb file.