X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fcluster%2Fwham%2Fsrc-M-NEWCORR%2Freadrtns.F;h=326d4ec8dc79e61fb5e0efe804bbc643c7b02828;hb=a09fdb5f3686a0c1242c8dd3b6dfadf0d3678aaf;hp=321e11e719d3f9b0ecf4025afaeaf2f0df7250e3;hpb=5467f8060e885151d3415a9897913776c337d88a;p=unres.git

diff --git a/source/cluster/wham/src-M-NEWCORR/readrtns.F b/source/cluster/wham/src-M-NEWCORR/readrtns.F
index 321e11e..326d4ec 100644
--- a/source/cluster/wham/src-M-NEWCORR/readrtns.F
+++ b/source/cluster/wham/src-M-NEWCORR/readrtns.F
@@ -64,6 +64,7 @@ C
       write (iout,*) 'beta_h',(beta_h(i),i=1,nT)
       lprint_cart=index(controlcard,"PRINT_CART") .gt.0
       lprint_int=index(controlcard,"PRINT_INT") .gt.0
+      with_dihed_constr = index(controlcard,"WITH_DIHED_CONSTR").gt.0
       if (min_var) iopt=1
       return
       end
@@ -87,6 +88,7 @@ C
       include 'COMMON.CONTROL'
       include 'COMMON.CONTACTS'
       include 'COMMON.TIME1'
+      include 'COMMON.TORCNSTR'
 #ifdef MPL
       include 'COMMON.INFO'
 #endif
@@ -193,7 +195,11 @@ C 12/1/95 Added weight for the multi-body term WCORR
         bad(i,2)=scalscp*bad(i,2)
       enddo
 
+#ifdef AIX
+      call flush_(iout)
+#else
       call flush(iout)
+#endif
       print *,'indpdb=',indpdb,' pdbref=',pdbref
 
 C Read sequence if not taken from the pdb file.
@@ -233,6 +239,28 @@ C Convert sequence to numeric code
 
       print *,'Call Read_Bridge.'
       call read_bridge
+
+      if (with_dihed_constr) then
+
+      read (inp,*) ndih_constr
+      write (iout,*) "ndih_constr",ndih_constr
+      if (ndih_constr.gt.0) then
+        read (inp,*) ftors
+        write (iout,*) 'FTORS',ftors
+        read (inp,*) (idih_constr(i),phi0(i),drange(i),i=1,ndih_constr)
+        write (iout,*)
+     &   'There are',ndih_constr,' constraints on phi angles.'
+        do i=1,ndih_constr
+          write (iout,'(i5,2f8.3)') idih_constr(i),phi0(i),drange(i)
+        enddo
+        do i=1,ndih_constr
+          phi0(i)=deg2rad*phi0(i)
+          drange(i)=deg2rad*drange(i)
+        enddo
+      endif
+
+      endif
+
       nnt=1
       nct=nres
       print *,'NNT=',NNT,' NCT=',NCT