const BoxArray &ba = boxArray(lev);
  const Box &bx = ba.minimalBox();

  const Box &dom = Geom(lev).Domain();

    int sz = bx.bigEnd(d) - bx.smallEnd(d) + 1;
    nbx.setSmall(d, bx.smallEnd(d) + sz / 4 + 1);
    nbx.setBig(d, bx.bigEnd(d) - sz / 4 - 1);

    int md = (dom.bigEnd(d) + dom.smallEnd(d) + 1) / 2;
    int rd = (dom.bigEnd(d) - dom.smallEnd(d) + 1) / 2;
    nbx.setSmall(d, md - rd / (1 << (lev + 1)));
    nbx.setBig(d, md + rd / (1 << (lev + 1)) - 1);

  cout << "lev: " << lev << "\n";

  const BoxArray &ba = boxArray(lev);

  const int maxnumlevels = 2;

  const int maxnumlevels = 3;

      assert(new_finest == level);

      // cout << "level=" << level << " new_finest=" << new_finest << "\n";
      // assert(new_finest == level);