X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2FMREMD.F90;h=37c772052894c92c3cace6e7e22607a2384b3dc3;hb=bc23440fbe68672d430f71f22f46b11265f003db;hp=8c4a4b5057badfc9a01f5592961f70728be55d94;hpb=f458de69c692cc132fdb9adfa1a130fc33b35782;p=unres4.git

diff --git a/source/unres/MREMD.F90 b/source/unres/MREMD.F90
index 8c4a4b5..37c7720 100644
--- a/source/unres/MREMD.F90
+++ b/source/unres/MREMD.F90
@@ -108,7 +108,7 @@
       integer :: k,il,il1,i,j,nharp,ii,ierr,itime_master,irr,iex,&
             i_set_temp,itmp,i_temp,i_mult,i_iset,i_mset,i_dir,i_temp1,&
             i_mult1,i_iset1,i_mset1,ierror,itime,mnum
-      integer,dimension(4,nres/3) :: iharp	!(4,nres/3)(4,maxres/3)
+      integer,dimension(4,nres) :: iharp	!(4,nres/3)(4,maxres/3)
 !deb      imin_itime_old=0
       integer :: nres2	!el
       WRITE(iout,*) "JUST AFTER CALL"
@@ -411,11 +411,13 @@
            do i=nnt,nct-1
              mnum=(molnum(i))
              stdforcp(i)=stdfp(mnum)*dsqrt(gamp(mnum))
+!             write(iout,*) "stdforcp=",stdforcp(i),itype(i,mnum),i
            enddo
            do i=nnt,nct
             mnum=molnum(i)
              if (itype(i,mnum).ne.ntyp1) stdforcsc(i)=stdfsc(iabs(itype(i,mnum)),mnum)&
                         *dsqrt(gamsc(iabs(itype(i,mnum)),mnum))
+!              write(iout,*) "stdforcsc=",stdforcsc(i),itype(i,mnum),i
            enddo
          endif
        call rescale_weights(t_bath)