X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc-HCD-5D%2Felecont.f;h=f43b037eed0e82878622ed86d21bc19e352fa36d;hb=c711143ad3fffb04d27b55aa823f399b8343c4c5;hp=73325f2fa65e9a380e1e9406c1a826983c405528;hpb=76ef494efde78d2d85d0e72d936c13166961256c;p=unres.git diff --git a/source/unres/src-HCD-5D/elecont.f b/source/unres/src-HCD-5D/elecont.f index 73325f2..f43b037 100644 --- a/source/unres/src-HCD-5D/elecont.f +++ b/source/unres/src-HCD-5D/elecont.f @@ -82,6 +82,8 @@ c data epp / 0.3045d0, 0.3649d0, 0.3649d0, 0.5743d0/ zj=mod(zj,boxzsize) if (zj.lt.0) zj=zj+boxzsize dist_init=(xj-xi)**2+(yj-yi)**2+(zj-zi)**2 + write (iout,*) "i",i,xi,yi,zi," j",j,xj,yj,xj,"dist", + & dsqrt(dist_init) xj_safe=xj yj_safe=yj zj_safe=zj @@ -141,6 +143,7 @@ c data epp / 0.3045d0, 0.3649d0, 0.3649d0, 0.5743d0/ endif ees=ees+eesij evdw=evdw+evdwij*sss + write (iout,*) "i"," j",j," rij",dsqrt(rij)," eesij",eesij 4 continue 1 continue if (lprint) then @@ -251,7 +254,7 @@ c-------------------------------------------- include 'DIMENSIONS' include 'COMMON.CHAIN' include 'COMMON.IOUNITS' - include 'COMMON.DISTFIT' + include 'COMMON.FRAG' include 'COMMON.VAR' include 'COMMON.GEO' include 'COMMON.CONTROL'