X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fsc_move.F;h=2082e98ef02e27bca34940773cd2327619f6e162;hb=49b4d93d38312bdd6e5e9c6bb845431396c9b0f5;hp=a7a4f6484f12125a730f9c331cb48cf2fb0aef31;hpb=478a9d9a1c99eb3f4bc4ca676ff3162bdd01d633;p=unres.git

diff --git a/source/unres/src_MD-M/sc_move.F b/source/unres/src_MD-M/sc_move.F
index a7a4f64..2082e98 100644
--- a/source/unres/src_MD-M/sc_move.F
+++ b/source/unres/src_MD-M/sc_move.F
@@ -211,7 +211,7 @@ c     Define what is meant by "neighbouring side-chain"
 
 c     Don't do glycine or ends
       i=itype(res_pick)
-      if (i.eq.10 .or. i.eq.21) return
+      if (i.eq.10 .or. i.eq.ntyp1) return
 
 c     Freeze everything (later will relax only selected side-chains)
       mask_r=.true.
@@ -253,7 +253,7 @@ cd      print *,'new       ',(energy(k),k=0,n_ene)
       n_try=0
       do while (n_try.lt.n_maxtry .and. orig_e-cur_e.lt.e_drop)
 c     Move the selected residue (don't worry if it fails)
-        call gen_side(itype(res_pick),theta(res_pick+1),
+        call gen_side(iabs(itype(res_pick)),theta(res_pick+1),
      +       alph(res_pick),omeg(res_pick),fail)
 
 c     Minimize the side-chains starting from the new arrangement
@@ -717,8 +717,8 @@ c     if (icall.eq.0) lprn=.true.
       do i=iatsc_s,iatsc_e
 
 
-        itypi=itype(i)
-        itypi1=itype(i+1)
+        itypi=iabs(itype(i))
+        itypi1=iabs(itype(i+1))
         xi=c(1,nres+i)
         yi=c(2,nres+i)
         zi=c(3,nres+i)
@@ -733,7 +733,7 @@ C
           do j=istart(i,iint),iend(i,iint)
           IF (mask_side(j).eq.1.or.mask_side(i).eq.1) THEN
             ind=ind+1
-            itypj=itype(j)
+            itypj=iabs(itype(j))
             dscj_inv=dsc_inv(itypj)
             sig0ij=sigma(itypi,itypj)
             chi1=chi(itypi,itypj)