X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2FREMD.F90;h=fc0c9f10fb883b08e6484de9cc9c8fecc3e2e606;hb=d4033f3ee5e94f3f22692d69fea70e36f231ebd1;hp=03147da04528eb208540fb0729f530b832fc909f;hpb=df2469d9ac903d93889867f4e50e9bf6c428c1c6;p=unres4.git

diff --git a/source/unres/REMD.F90 b/source/unres/REMD.F90
index 03147da..fc0c9f1 100644
--- a/source/unres/REMD.F90
+++ b/source/unres/REMD.F90
@@ -357,9 +357,11 @@
          nind=2
        endif	  
        do i=nnt+1,nct-1
+        mnum=molnum(i)
 !         if (iabs(itype(i,1)).eq.ntyp1) cycle
          DM(ind)=2*ip4+mp/2
-         if (iabs(itype(i,1)).eq.10 .or. iabs(itype(i,1)).eq.ntyp1) then
+        if (iabs(itype(i,1)).eq.10 .or. &
+          iabs(itype(i,mnum)).eq.ntyp1_molec(mnum) .or. mnum.eq.5) then
            if (iabs(itype(i,1)).eq.10) DM(ind)=DM(ind)+msc(10)
            ind=ind+1
          else
@@ -426,8 +428,10 @@
        endif     
        do i=nnt+1,nct-1
 !         if (iabs(itype(i,1)).eq.ntyp1) cycle
+         mnum=molnum(i)
          DM(ind)=2*ip4+mp(1)/2
-         if (iabs(itype(i,1)).eq.10 .or. iabs(itype(i,1)).eq.ntyp1) then
+       if (iabs(itype(i,1)).eq.10 .or. &
+          iabs(itype(i,mnum)).eq.ntyp1_molec(mnum) .or. mnum.eq.5) then
            if (iabs(itype(i,1)).eq.10) DM(ind)=DM(ind)+msc(10,1)
            ind=ind+1
          else