etor & etor_d multichain match prerelease

[unres4.git] / source / unres / energy.f90

diff --git a/source/unres/energy.f90 b/source/unres/energy.f90
index 4bb6804..7960906 100644 (file)
--- a/source/unres/energy.f90
+++ b/source/unres/energy.f90
@@ -5537,6 +5537,7 @@
       etors=0.0D0
       do i=iphi_start,iphi_end
         if (itype(i-2).eq.ntyp1 .or. itype(i-1).eq.ntyp1 &
+             .or. itype(i-3).eq.ntyp1 &
              .or. itype(i).eq.ntyp1) cycle
         etors_ii=0.0D0
          if (iabs(itype(i)).eq.20) then
@@ -5647,6 +5648,7 @@
       do i=iphid_start,iphid_end
         etors_d_ii=0.0D0
         if (itype(i-2).eq.ntyp1 .or. itype(i-1).eq.ntyp1 &
+            .or. itype(i-3).eq.ntyp1 &
             .or. itype(i).eq.ntyp1 .or. itype(i+1).eq.ntyp1) cycle
         itori=itortyp(itype(i-2))
         itori1=itortyp(itype(i-1))