X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fstochfric.F;h=10531b4db9d25cb2ee43e9bef9270209f48b427c;hb=5ba2a6014fbd859202ed28885b66520ac09bb77d;hp=cda93c394537d95b91e019f15e3d35ab3fbfb5d6;hpb=2acc9917f7452c3a59a242e6033bea16f846cbfe;p=unres.git

diff --git a/source/unres/src_MD-M/stochfric.F b/source/unres/src_MD-M/stochfric.F
index cda93c3..10531b4 100644
--- a/source/unres/src_MD-M/stochfric.F
+++ b/source/unres/src_MD-M/stochfric.F
@@ -39,7 +39,7 @@
         ind=ind+3
       enddo
       do i=nnt,nct
-        if (itype(i).ne.10) then
+        if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
           do j=1,3
             d_t_work(ind+j)=d_t(j,i+nres)
           enddo
@@ -68,7 +68,7 @@
         ind=ind+3
       enddo
       do i=nnt,nct
-        if (itype(i).ne.10) then
+        if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
           do j=1,3
             friction(j,i+nres)=fric_work(ind+j)
           enddo
@@ -234,7 +234,7 @@ c Compute the stochastic forces acting on virtual-bond vectors.
         stochforc(j,0)=ff(j)+force(j,nnt+nres)
       enddo
       do i=nnt,nct
-        if (itype(i).ne.10) then
+        if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
           do j=1,3
             stochforc(j,i+nres)=force(j,i+nres)
           enddo
@@ -252,7 +252,7 @@ c Compute the stochastic forces acting on virtual-bond vectors.
         ind=ind+3
       enddo
       do i=nnt,nct
-        if (itype(i).ne.10) then
+        if ((itype(i).ne.10).and.(itype(i).ne.ntyp1)) then
           do j=1,3
             stochforcvec(ind+j)=stochforc(j,i+nres)
           enddo
@@ -365,6 +365,7 @@ c  Load the friction coefficients corresponding to peptide groups
 c  Load the friction coefficients corresponding to side chains
       m=nct-nnt
       ind=0
+C      gamsc(ntyp1)=1.0d0
       do i=nnt,nct
         ind=ind+1
         ii = ind+m
@@ -556,7 +557,7 @@ c
       include 'COMMON.INTERACT'
       include 'COMMON.IOUNITS'
       include 'COMMON.NAMES'
-      double precision radius(maxres2),gamvec(maxres6)
+      double precision radius(maxres2),gamvec(maxres2)
       parameter (twosix=1.122462048309372981d0)
       logical lprn /.false./
 c