parameter(n1=32,n2=32,n3=32)
       implicit real*8(a-h,o-z), integer*4(i-n)


       parameter(pi=3.14159)

       common /grid/  x1a(-1:n1+2), x2a(-1:n2+2), x3a(-1:n3+2),
     +                x1b(-1:n1+1), x2b(-1:n2+1), x3b(-1:n3+1),
     +               dx1a(-1:n1+1),dx2a(-1:n2+1),dx3a(-1:n3+1),
     +               dx1b(-1:n1+0),dx2b(-1:n2+0),dx3b(-1:n3+0),
     +               dv(n1,n2,n3)

       common /variables/  d(-1:n1+2,-1:n2+2,-1:n3+2),
     +                    v1(-1:n1+2,-1:n2+2,-1:n3+2),
     +                    v2(-1:n1+2,-1:n2+2,-1:n3+2),
     +                    v3(-1:n1+2,-1:n2+2,-1:n3+2),
     +                     p(-1:n1+2,-1:n2+2,-1:n3+2),
     +                     e(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   phi(-1:n1+2,-1:n2+2,-1:n3+2)


       common /ds/       ds1(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   ds2(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   ds3(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fm1(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fm2(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fm3(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fe1(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fe2(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   Fe3(-1:n1+2,-1:n2+2,-1:n3+2),  
     +                   vs(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   es1(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   es2(-1:n1+2,-1:n2+2,-1:n3+2),
     +                   es3(-1:n1+2,-1:n2+2,-1:n3+2)


       common /bdary/ var(-1:n1+2,-1:n2+2,-1:n3+2)

       common /potential/ Glm(-n1:n1-1,-n2:n2-1,-n3:n3-1),eps
       common /pot test/  Plm(-1:n1+2,-1:n2+2,-1:n3+2)


       common /information/ pk,gamma,time,dt,cf,cv,trun,tsmin,
     +   time_restart,outint,rm,x1min,x1max,x2min,x2max,x3min,x3max
       common /information2/ istart,nstep,nsteps,istep,ndim,istrt,
     +                       istart_restart,iacc