integer nstep,nacc,nblock,ndim,np,nslices,nanal,nspill,iwrap,iwra
      integer ikin,iepot,ietot,irg,ieff,igr,ln,nmovers,nrap,level,ntry
      integer nalg,ntrial,iev,ntrydsp,naccdsp,ilcut,iord,idsp,nav,ilind
      integer nin,nout,nstart,ncumtr,idbg,ifewald,jlindem,ipot,io2,io2nl
      integer isofk,nshlls,kmult(0:mnsh),nvects,nkpair,isk,nlamb,ngvects
      integer iph(iomax,iomax),isg(iomax,iomax),iosg
      real*8 x,temp,beta,dxmax,mratio,epsph,pi,spi,sq2pi,sq2,ax,ay,pbc
      real*8 ell(mdim),tpiell(mdim),el2(mdim),cutr,vol,rho,rs,cutr2
      real*8 alpha,hop,gammadsp,a12,a11,alp,rnn,proton
      real*8 amass,charge,omega,alamb,tau,cke,sigma,sigma2,dpi,dr
      real*8 rdsp(ldsptable),dlogr,deltae,emin
      real*8 spring(mparts,mslices),sold(mmovers),xt
      real*8 soldm(mmovers),action,cpot,epact,ep12,epact1
      real*8 cpotew(mslices),epactew(mslices,mslices)
      real*8 phprop(mslices,iom2nl),dphprop(mslices,iom2nl)
      real*8 phpvar(mslices),dphpvar(mslices),amu,xtr
      real*8 phpvar2(mslices),udsp(mslices)
      real*8 av(mav),anorm(mav),avp(mav),anormp(mav)
      real*8 d1,d2,d3,cost,omt,d4,dspcut
      real*8 xnode(iomax),wg(iomax),up2(iomax),um2(iomax)
      real*8 chi(mdim,mparts,mslices),phi(mdim,mparts,mslices)
      real*8 aja(mdim,mparts,mslices),bja(mdim,mparts,mslices)
      real*8 di(lptable),aself(lptable,2),epself(lptable,2)
      real*8 evself(lptable,2),rptable(6),vsrtab(lptable,2,iom2)
      real*8 ssrtab(lptable,2,iom2),ssreph(lptable,2,iomsg)
      real*8 vsreph(lptable,2,iomsg),rptself(6)
      real*8 cf,wf,agx,ahx,asiguv,costact
      real*8 xb(mdim,mparts,mslices),cxb(mdim,mparts,mslices)
      real*8 xd(mdim,mparts,mslices)
      real*8 xbold(mdim,mparts,mslices),cxbold(mdim,mparts,mslices)
      real*8 xdold(mdim,mparts,mslices)
      real*8 avar(iomax),bvar(iomax)
      real*8 sdvar(iomax,iomax),ssvar(iomax,iomax)
      real*8 cutk,rkcomp(mdim,mnkv),rknorm(mnsh),rhok(mnkv2,mslices)
      real*8 pwmat(mnkv2,mparts,mslices),vk(mnsh),vkconst
      real*8 rhokb(mnkv2,mslices),pwmatb(mnkv2,mparts,mslices)
      real*8 rsites(mdim,mparts),bond,zp
      real*4 tt,ttanal
      character filen*14,title*80
      logical ifex
      common/st/np,nslices,ndim,ln,nblock,nstep,nacc,nmovers,nrap,level
     &         ,ntrial,ntry,nalg,ntrydsp,naccdsp,nin,nout,nstart
      common/pos/x(mdim,mparts,mslices),cke,sigma,action,cpot,epact,ep12
     &          ,hop,gammadsp,a12,a11,sigma2,epact1,ell,tpiell,el2,cutr
     &          ,vol,rho,rs,cpotew,epactew,rsites,cutr2,proton(mdim,2)
     &          ,bond,zp
      common/ave/av,anorm,avp,anormp,tt,ttanal,dpi,dr,rdsp,dlogr
     &          ,deltae,rnn,ikin,iepot,ietot,irg,ieff,igr,iev,idsp,nav
     &          ,ilind,nanal,nspill,jlindem,filen
      common/pph/rptself,aself,epself,evself,phprop,dphprop,alpha
     &         ,charge,omega,beta,temp,tau,alamb,epsph,pi,spi,sq2pi,sq2
     &         ,udsp,di,alp,ilcut,iord,io2,io2nl,iph,isg,iosg
      common/phcost/d1,d2,d3,cost,omt,d4
      common/gauint/xnode,wg,up2,um2,chi,phi,aja,bja
      common/fvar/xb,cxb,xd,xbold,cxbold,xdold,avar,bvar,sdvar,ssvar
     &           ,phpvar2,ncumtr,idbg,cf,wf,phpvar,dphpvar,amu,xtr,agx
     &           ,ahx,asiguv,costact
      common/kspace/rkcomp,rknorm,rhok,pwmat,vk,rhokb,pwmatb,vkconst
     &             ,cutk,nlamb,ngvects,isofk,ifewald,nshlls,nvects
     &             ,nkpair,kmult,isk
      common/tables/vsrtab,ssrtab,ssreph,vsreph,rptable,ipot
      iwrap(iwra)=1+mod(iwra+nslices-1,nslices)
      pbc(ax,ay,iwra)=ax-ay
c     pbc(ax,ay,iwra)=ax-ay-anint((ax-ay)/ell(iwra))*ell(iwra)