- xj=boxshift(xj-xi,boxxsize)
- yj=boxshift(yj-yi,boxysize)
- zj=boxshift(zj-zi,boxzsize)
+c write (iout,*) "j",xj,yj,zj
+ xj=boxshift(xj-xmedi,boxxsize)
+ yj=boxshift(yj-ymedi,boxysize)
+ zj=boxshift(zj-zmedi,boxzsize)
+c write (iout,*) "j",xj,yj,zj