X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?p=unres.git;a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fenergy_p_new_barrier.F;h=3703b07d990641bc33b36092912d6ca369dad53d;hp=e3a17dd1d69b8645e431f6d8003e0c3c5bda09f5;hb=e12c8db731e266c6e38b2f5883f3ae7c15bbfb98;hpb=e64f51efa7130a5e8e2ade165f9a7615e084bba0

diff --git a/source/unres/src_MD-M/energy_p_new_barrier.F b/source/unres/src_MD-M/energy_p_new_barrier.F
index e3a17dd..3703b07 100644
--- a/source/unres/src_MD-M/energy_p_new_barrier.F
+++ b/source/unres/src_MD-M/energy_p_new_barrier.F
@@ -3692,13 +3692,13 @@ c        if ((zmedi.gt.((0.5d0)*boxzsize)).or.
 c     &       (zmedi.lt.((-0.5d0)*boxzsize))) then
 c        go to 196
 c        endif
-          xmedi=mod(xmedi,boxxsize)
+          xmedi=dmod(xmedi,boxxsize)
           if (xmedi.lt.0) xmedi=xmedi+boxxsize
-          ymedi=mod(ymedi,boxysize)
+          ymedi=dmod(ymedi,boxysize)
           if (ymedi.lt.0) ymedi=ymedi+boxysize
-          zmedi=mod(zmedi,boxzsize)
+          zmedi=dmod(zmedi,boxzsize)
           if (zmedi.lt.0) zmedi=zmedi+boxzsize
-          zmedi2=mod(zmedi,boxzsize)
+          zmedi2=dmod(zmedi,boxzsize)
           if (zmedi2.lt.0) zmedi2=zmedi2+boxzsize
        if ((zmedi2.gt.bordlipbot)
      &.and.(zmedi2.lt.bordliptop)) then
@@ -3757,11 +3757,11 @@ c     &  .or. itype(i-1).eq.ntyp1
         xmedi=c(1,i)+0.5d0*dxi
         ymedi=c(2,i)+0.5d0*dyi
         zmedi=c(3,i)+0.5d0*dzi
-          xmedi=mod(xmedi,boxxsize)
+          xmedi=dmod(xmedi,boxxsize)
           if (xmedi.lt.0) xmedi=xmedi+boxxsize
-          ymedi=mod(ymedi,boxysize)
+          ymedi=dmod(ymedi,boxysize)
           if (ymedi.lt.0) ymedi=ymedi+boxysize
-          zmedi=mod(zmedi,boxzsize)
+          zmedi=dmod(zmedi,boxzsize)
           if (zmedi.lt.0) zmedi=zmedi+boxzsize
        if ((zmedi.gt.bordlipbot)
      &.and.(zmedi.lt.bordliptop)) then
@@ -3965,6 +3965,7 @@ C lipbufthick is thickenes of lipid buffore
        enddo
        enddo
        if (isubchap.eq.1) then
+C          print *,i,j
           xj=xj_temp-xmedi
           yj=yj_temp-ymedi
           zj=zj_temp-zmedi
@@ -4173,10 +4174,14 @@ C           print *,"bafter", gelc_long(1,i), gelc_long(1,j)
 C Lipidic part for lipscale
             gelc_long(3,j)=gelc_long(3,j)+
      &     ssgradlipj*eesij/2.0d0*lipscale**2
-
+C           if ((ssgradlipj*eesij/2.0d0*lipscale**2).ne.0.0 )
+C     &     write(iout,*) "WTF",j
             gelc_long(3,i)=gelc_long(3,i)+
      &     ssgradlipi*eesij/2.0d0*lipscale**2
 
+C            if ((ssgradlipi*eesij/2.0d0*lipscale**2).ne.0.0 )
+C     &     write(iout,*) "WTF",i
+
 *
 * Loop over residues i+1 thru j-1.
 *