X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Fmoments.F;h=4f57331bd7623c9939f876b117a40bb5f92df579;hb=5836ecdab5a8b95f079bbf6e07374dee3fce8a26;hp=bbc0f27aed63acdeb3c355e88a111a97e095c4a5;hpb=a05e8996b6fb955dc21aef7a85db21ee73d0a9d3;p=unres.git

diff --git a/source/unres/src-HCD-5D/moments.F b/source/unres/src-HCD-5D/moments.F
index bbc0f27..4f57331 100644
--- a/source/unres/src-HCD-5D/moments.F
+++ b/source/unres/src-HCD-5D/moments.F
@@ -183,7 +183,8 @@ c   Resetting the velocities
 #ifdef FIVEDIAG
       do i=nnt,nct-1
         write (iout,*) itype(i+1),itype(i)
-        if (itype(i+1).ne.ntyp1 .and. itype(i).eq.ntyp1) cycle
+        if (itype(i+1).ne.ntyp1 .and. itype(i).eq.ntyp1 .or.
+     &      itype(i).ne.ntyp1 .and. itype(i+1).eq.ntyp1) cycle
         call vecpr(vrot(1),dc(1,i),vp)  
         do j=1,3
           d_t(j,i)=d_t(j,i)-vp(j)