X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fsc_move.F;h=02a329db5ef03bc8a6c2a25c6ddfa6b737bc98ac;hb=5316f1d82cb17d390be734cefbfca809cc65dcc4;hp=74e9bf2d00a941d45451cd8ba21f24b9f4dbb78b;hpb=746b71e825ca1fada9877a347562459b8390629d;p=unres.git diff --git a/source/unres/src_MD/sc_move.F b/source/unres/src_MD/sc_move.F index 74e9bf2..02a329d 100644 --- a/source/unres/src_MD/sc_move.F +++ b/source/unres/src_MD/sc_move.F @@ -213,7 +213,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. @@ -255,7 +255,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(iabs(itype(res_pick)),theta(res_pick+1), + call gen_side(itype(res_pick),theta(res_pick+1), + alph(res_pick),omeg(res_pick),fail) c Minimize the side-chains starting from the new arrangement @@ -719,8 +719,8 @@ c if (icall.eq.0) lprn=.true. do i=iatsc_s,iatsc_e - itypi=iabs(itype(i)) - itypi1=iabs(itype(i+1)) + itypi=itype(i) + itypi1=itype(i+1) xi=c(1,nres+i) yi=c(2,nres+i) zi=c(3,nres+i) @@ -735,7 +735,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=iabs(itype(j)) + itypj=itype(j) dscj_inv=dsc_inv(itypj) sig0ij=sigma(itypi,itypj) chi1=chi(itypi,itypj)