Fixes to compile WHAM and CLUSTER on AIX (nostromo.icm.edu.pl)

[unres.git] / source / unres / src_MD / readrtns.F

diff --git a/source/unres/src_MD/readrtns.F b/source/unres/src_MD/readrtns.F
index 636c6a4..918c761 100644 (file)
--- a/source/unres/src_MD/readrtns.F
+++ b/source/unres/src_MD/readrtns.F
@@ -1087,6 +1087,7 @@ C     Juyong:READ read_info
 C     READ fragment information!!
 C     both routines should be in dfa.F file!!
 
+#ifdef DFA
       if (.not. (wdfa_dist.eq.0.0 .and. wdfa_tor.eq.0.0 .and.
      &            wdfa_nei.eq.0.0 .and. wdfa_beta.eq.0.0)) then
        call init_dfa_vars
@@ -1094,6 +1095,7 @@ C     both routines should be in dfa.F file!!
        call read_dfa_info
        print*, 'read_dfa_info finished!'
       endif
+#endif
 C
       if (pdbref) then
         if(me.eq.king.or..not.out1file)
@@ -2508,33 +2510,36 @@ c-------------------------------------------------------------------------------
       include 'COMMON.IOUNITS'
       include 'COMMON.MD'
       include 'COMMON.CONTROL'
+      integer iset1
       read(inp,*) nset,nfrag,npair,nfrag_back
       if(me.eq.king.or..not.out1file)
      & write(iout,*) "nset",nset," nfrag",nfrag," npair",npair,
      &  " nfrag_back",nfrag_back
-      do iset=1,nset
-         read(inp,*) mset(iset)
+      do iset1=1,nset
+         read(inp,*) mset(iset1)
        do i=1,nfrag
-         read(inp,*) wfrag(i,iset),ifrag(1,i,iset),ifrag(2,i,iset), 
-     &     qinfrag(i,iset)
+         read(inp,*) wfrag(i,iset1),ifrag(1,i,iset1),ifrag(2,i,iset1), 
+     &     qinfrag(i,iset1)
          if(me.eq.king.or..not.out1file)
-     &    write(iout,*) "R ",i,wfrag(i,iset),ifrag(1,i,iset),
-     &     ifrag(2,i,iset), qinfrag(i,iset)
+     &    write(iout,*) "R ",i,wfrag(i,iset1),ifrag(1,i,iset1),
+     &     ifrag(2,i,iset1), qinfrag(i,iset1)
        enddo
        do i=1,npair
-        read(inp,*) wpair(i,iset),ipair(1,i,iset),ipair(2,i,iset), 
-     &    qinpair(i,iset)
+        read(inp,*) wpair(i,iset1),ipair(1,i,iset1),ipair(2,i,iset1), 
+     &    qinpair(i,iset1)
         if(me.eq.king.or..not.out1file)
-     &   write(iout,*) "R ",i,wpair(i,iset),ipair(1,i,iset),
-     &    ipair(2,i,iset), qinpair(i,iset)
+     &   write(iout,*) "R ",i,wpair(i,iset1),ipair(1,i,iset1),
+     &    ipair(2,i,iset1), qinpair(i,iset1)
        enddo 
        do i=1,nfrag_back
-        read(inp,*) wfrag_back(1,i,iset),wfrag_back(2,i,iset),
-     &     wfrag_back(3,i,iset),
-     &     ifrag_back(1,i,iset),ifrag_back(2,i,iset)
+        read(inp,*) wfrag_back(1,i,iset1),wfrag_back(2,i,iset1),
+     &     wfrag_back(3,i,iset1),
+     &     ifrag_back(1,i,iset1),ifrag_back(2,i,iset1)
         if(me.eq.king.or..not.out1file)
-     &   write(iout,*) "A",i,wfrag_back(1,i,iset),wfrag_back(2,i,iset),
-     &   wfrag_back(3,i,iset),ifrag_back(1,i,iset),ifrag_back(2,i,iset)
+     &   write(iout,*) "A",i,wfrag_back(1,i,iset1),
+     &   wfrag_back(2,i,iset1),
+     &   wfrag_back(3,i,iset1),ifrag_back(1,i,iset1),
+     &   ifrag_back(2,i,iset1)
        enddo 
       enddo
       return
@@ -2841,7 +2846,7 @@ c    &       sigma_odl_temp(i+nnt-1,j+nnt-1,k)
 c         enddo
 c 1401   continue
 c         close (ientin)
-        if (waga_dist.gt.0.0d0) then
+        if (waga_dist.ne.0.0d0) then
           ii=0
           do i = nnt,nct-2 ! right? without parallel.
             do j=i+2,nct ! right?
@@ -2881,8 +2886,13 @@ c
               sigma_odl(k,ii)=rescore(k,i)+rescore(k,j) ! other exprs possible
 c             sigma_odl(k,ii)=hmscore(k)*rescore(k,i)*rescore(k,j)
             else
+#ifdef OLDSIGMA
               sigma_odl(k,ii)=(rescore(k,i)+rescore(k,j))* ! sigma ~ rescore ~ error 
      &                      dexp(0.5d0*(odl(k,ii)/dist_cut)**2)
+#else
+              sigma_odl(k,ii)=(rescore(k,i)+rescore(k,j))* ! sigma ~ rescore ~ error 
+     &                      dexp(0.5d0*(odl(k,ii)/dist_cut)**2-0.5d0)
+#endif
 
 c   Following expr replaced by a positive exp argument
 c             sigma_odl(k,ii)=(rescore(k,i)+rescore(k,j))*
@@ -3006,7 +3016,7 @@ c              read (ientin,*) sigma_d(k,i) ! 1st variant
         endif
         close(ientin)
       enddo
-      if (waga_dist.gt.0.0d0) lim_odl=ii
+      if (waga_dist.ne.0.0d0) lim_odl=ii
       if (constr_homology.gt.0) call homology_partition
       if (constr_homology.gt.0) call init_int_table
 cd      write (iout,*) "homology_partition: lim_theta= ",lim_theta,