c..common block declarations 
      integer         nd1,ndp8
      parameter       (nd1  = 1000)
      parameter       (ndp8 = 1008)


      integer          nr, n, np1, np2, np3, np4, np5, np6,
     1                 np7, np8, nstep, nmax, nriem, alpha
      common / intgr / nr, n, np1, np2, np3, np4, np5, np6,
     1                 np7, np8, nstep, nmax, nriem, alpha


c..grid 
      double precision rl(nd1), r(nd1), rr(nd1), 
     2                 xl(ndp8), x(ndp8), xr(ndp8), dx(ndp8),
     3                 xlnu(ndp8), xnu(ndp8), xrnu(ndp8), dxnu(ndp8),
     4                 dm(ndp8), dmnu(ndp8), dtdx(ndp8), dtdm(ndp8)

      common / grd   / rl, r, rr, 
     2                 xl, x, xr, dx,
     3                 xlnu, xnu, xrnu, dxnu,
     4                 dm, dmnu, dtdx, dtdm


c..run time parameters
      double precision t, tmax, dt, dtmin, dtmax, cfl, dy,
     1                 gamma, gamm1, gamp1, gmfc, gamfac, g, gr,
     2                 small, smlrho, smallp, smalle

      common / float / t, tmax, dt, dtmin, dtmax, cfl, dy,
     1                 gamma, gamm1, gamp1, gmfc, gamfac, g, gr,
     2                 small, smlrho, smallp, smalle


c..multi-dimensional arrays
      double precision densty(nd1), vr(nd1),
     1                 energy(nd1), press(nd1)

      common / twod  / densty, vr,
     1                 energy, press


c..1d ppm arrays
      double precision tau(ndp8),   u(ndp8),   e(ndp8),
     1                 p(ndp8), taunu(ndp8), unu(ndp8),
     2                 enu(ndp8), pnu(ndp8), rho(ndp8), rhonu(ndp8), 
     3                 c(ndp8),  ce(ndp8)

      common / oned  / tau,   u,   e,
     1                 p, taunu, unu,
     2                 enu, pnu, rho, rhonu, 
     3                 c,  ce


c..for the riemann solve
      double precision taupls(ndp8), upls(ndp8), ppls(ndp8), cpls(ndp8),
     1                 taumns(ndp8), umns(ndp8), pmns(ndp8), cmns(ndp8),
     2                 wpls(ndp8), wmns(ndp8),
     3                 pstar(ndp8), ustar(ndp8), upstar(ndp8)

      common / riemn / taupls, upls, ppls, cpls,
     1                 taumns, umns, pmns, cmns,
     2                 wpls, wmns,
     3                 pstar, ustar, upstar


      double precision vold(nd1)
      common / old  /  vold



c..ppm interpolation
      double precision c1(ndp8), c2(ndp8), c3(ndp8), c4(ndp8), c5(ndp8),
     1                 dela(ndp8), fshk(ndp8), fshk1(ndp8)

      common / intrp / c1, c2, c3, c4, c5,
     1                 dela, fshk, fshk1


      double precision taul(ndp8), taur(ndp8), dtau(ndp8), tau6(ndp8),
     1                 ul(ndp8), ur(ndp8), du(ndp8), u6(ndp8),
     3                 pl(ndp8), pr(ndp8), dp(ndp8), p6(ndp8),
     4                 el(ndp8), er(ndp8), de(ndp8), e6(ndp8),
     5                 rhol(ndp8), rhor(ndp8), drho(ndp8), rho6(ndp8)

      common / intfc / taul, taur, dtau, tau6,
     1                 ul, ur, du, u6,
     3                 pl, pr, dp, p6,
     4                 el, er, de, e6,
     5                 rhol, rhor, drho, rho6