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

diff --git a/source/unres/REMD.F90 b/source/unres/REMD.F90
index 03147da..699631c 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
@@ -623,11 +627,11 @@
         ii = ind+m
         mnum=molnum(i)
         iti=itype(i,mnum)
-        if (mnum.eq.5) then
-        mscab=0.0
-        else
+!        if (mnum.eq.5) then
+!        mscab=0.0
+!        else
         mscab=msc(iabs(iti),mnum)
-        endif
+!        endif
         massvec(ii)=mscab
 
         if (iti.ne.10 .and. iti.ne.ntyp1_molec(mnum) .and. mnum.ne.5) then