X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2FMD_A-MTS.F;h=d4d6be688b71f0664cb0610d2dab964b8b41a35c;hb=ff73ddbe323f5af85e31af14caa0ee6fc162a94a;hp=30d72a7dd64d8c7039ea7b064648af280b5d78ee;hpb=746b71e825ca1fada9877a347562459b8390629d;p=unres.git diff --git a/source/unres/src_MD/MD_A-MTS.F b/source/unres/src_MD/MD_A-MTS.F index 30d72a7..d4d6be6 100644 --- a/source/unres/src_MD/MD_A-MTS.F +++ b/source/unres/src_MD/MD_A-MTS.F @@ -1681,8 +1681,7 @@ c if the friction coefficients do not depend on surface area stdforcp(i)=stdfp*dsqrt(gamp) enddo do i=nnt,nct - stdforcsc(i)=stdfsc(iabs(itype(i))) - & *dsqrt(gamsc(iabs(itype(i)))) + stdforcsc(i)=stdfsc(itype(i))*dsqrt(gamsc(itype(i))) enddo endif c Open the pdb file for snapshotshots @@ -1858,7 +1857,7 @@ c Removing the velocity of the center of mass call chainbuild_cart call kinetic(EK) if (tbf) then - call verlet_bath(EK) + call verlet_bath endif kinetic_T=2.0d0/(dimen3*Rb)*EK if(me.eq.king.or..not.out1file)then