remd with homol_nset>1 working with FG parallelization

[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 b731cf2..69c7fe8 100644 (file)
--- a/source/unres/src_MD/readrtns.F
+++ b/source/unres/src_MD/readrtns.F
@@ -2508,33 +2508,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