X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fenergy_p_new_barrier.F;h=0af23e5475e6002018cd9dd0b69f1f8db84c8184;hb=3176770f951fd8b689c14df6260e55fa4f479689;hp=102862aad1d19b7119803dfb86d7475c331bca8c;hpb=39bffec2e960cdaa362c075425c0aa94bbd673c9;p=unres.git

diff --git a/source/unres/src_MD/energy_p_new_barrier.F b/source/unres/src_MD/energy_p_new_barrier.F
index 102862a..0af23e5 100644
--- a/source/unres/src_MD/energy_p_new_barrier.F
+++ b/source/unres/src_MD/energy_p_new_barrier.F
@@ -4279,13 +4279,13 @@ C 24/11/03 AL: SS bridges handled separately because of introducing a specific
 C    distance and angle dependent SS bond potential.
 cmc        if (ii.gt.nres .and. itype(iii).eq.1 .and. itype(jjj).eq.1) then
 C 18/07/06 MC: Use the convention that the first nss pairs are SS bonds
-c          if (.not.dyn_ss .and. i.le.nss) then
-C 15/02/13 CC dynamic SSbond
-        if (.not.dyn_ss.and.
-     &   ii.gt.nres .and. itype(iii).eq.1 .and. itype(jjj).eq.1) then 
+        if (.not.dyn_ss .and. i.le.nss) then
+C 15/02/13 CC dynamic SSbond - additional check
+         if (ii.gt.nres 
+     &       .and. itype(iii).eq.1 .and. itype(jjj).eq.1) then 
           call ssbond_ene(iii,jjj,eij)
           ehpb=ehpb+2*eij
-c          write (iout,*) "eij",eij
+         endif
         else if (ii.gt.nres .and. jj.gt.nres) then
 c Restraints from contact prediction
           dd=dist(ii,jj)