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
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
include 'DIMENSIONS'
include 'COMMON.CHAIN'
include 'COMMON.IOUNITS'
- include 'COMMON.DISTFIT'
+ include 'COMMON.FRAG'
include 'COMMON.VAR'
include 'COMMON.GEO'
include 'COMMON.CONTROL'