nosearchsc in src_MIN as in other versions

[unres.git] / source / unres / src_MIN / readrtns_min.F

diff --git a/source/unres/src_MIN/readrtns_min.F b/source/unres/src_MIN/readrtns_min.F
index dee929a..6352622 100644 (file)
--- a/source/unres/src_MIN/readrtns_min.F
+++ b/source/unres/src_MIN/readrtns_min.F
@@ -118,7 +118,8 @@ C Set up the time limit (caution! The time must be input in minutes!)
       dccart=(index(controlcard,'CART').gt.0)
       overlapsc=(index(controlcard,'OVERLAP').gt.0)
       overlapsc=.not.overlapsc
-      searchsc=(index(controlcard,'SEARCHSC').gt.0)
+      searchsc=(index(controlcard,'NOSEARCHSC').gt.0)
+      searchsc=.not.searchsc
       sideadd=(index(controlcard,'SIDEADD').gt.0)
       energy_dec=(index(controlcard,'ENERGY_DEC').gt.0)
       outpdb=(index(controlcard,'PDBOUT').gt.0)
@@ -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.