X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2FREMD.f90;h=3d11d5c69e72e8989bf7be998fe61d4991466302;hb=3cfa75a7cde45537ea8ebf322ae4d396b8e29ddc;hp=e2db47896619eb6a519ebd30fb2bccab65fb94e3;hpb=9dee9f9c465acfd5330503eeba8cf784b19c89f7;p=unres4.git

diff --git a/source/unres/REMD.f90 b/source/unres/REMD.f90
index e2db478..3d11d5c 100644
--- a/source/unres/REMD.f90
+++ b/source/unres/REMD.f90
@@ -452,8 +452,10 @@
         ind=1
         do i=nnt,nct
         mnum=molnum(i)
-       if (iabs(itype(i,1)).ne.10 .and.iabs(itype(i,mnum)).ne.ntyp1 &
-          .and.mnum.eq.5) then
+!       if (iabs(itype(i,1)).ne.10 .and.iabs(itype(i,mnum)).ne.ntyp1 &
+!          .and.mnum.eq.5) then
+       if (iabs(itype(i,1)).ne.10 .and. &
+          iabs(itype(i,mnum)).ne.ntyp1_molec(mnum) .and. mnum.ne.5) then
           DU1(ind)=-isc(iabs(itype(i,1)),1)
           DU1(ind+1)=0.0d0
          ind=ind+2
@@ -587,7 +589,7 @@
 !  Diagonal elements of the dC part of A and the respective friction coefficients
       ind=1
       ind1=0
-      print *,"TUTUTUT?!",nnt,nct-1
+!      print *,"TUTUTUT?!",nnt,nct-1
       do i=nnt,nct-1
         mnum=molnum(i)
         ind=ind+1
@@ -597,7 +599,7 @@
         print *,i,coeff,mp(mnum)
         massvec(ind1)=mp(mnum)
         Gmat(ind,ind)=coeff
-        print *,"i",mp(mnum)
+!        print *,"i",mp(mnum)
         A(ind1,ind)=0.5d0
       enddo
       
@@ -808,6 +810,7 @@
       endif
       deallocate(Gcopy)
 #endif
+!write(iout,*) "end setup_MD_matr"
       return
       end subroutine setup_MD_matrices
 !-----------------------------------------------------------------------------