const CCTK_REAL x = x0[0] + (lbnd[0] + i + CCTK_REAL(CI - 1) / 2) * dx[0];
      const CCTK_REAL y = x0[1] + (lbnd[1] + j + CCTK_REAL(CJ - 1) / 2) * dx[1];
      const CCTK_REAL z = x0[2] + (lbnd[2] + k + CCTK_REAL(CK - 1) / 2) * dx[2];

      const CCTK_REAL x = x0[0] + (lbnd[0] + i + CCTK_REAL(CI) / 2) * dx[0];
      const CCTK_REAL y = x0[1] + (lbnd[1] + j + CCTK_REAL(CJ) / 2) * dx[1];
      const CCTK_REAL z = x0[2] + (lbnd[2] + k + CCTK_REAL(CK) / 2) * dx[2];

    const CCTK_REAL origin_space = global_x0[d];

    const CCTK_REAL origin_space =
        global_x0[d] + (1 - nghostzones[d] + CCTK_REAL(1) / 2) * global_dx[d];

  const CCTK_REAL *restrict x0 = geom.ProbLo();
  const CCTK_REAL *restrict dx = geom.CellSize();

  const CCTK_REAL *restrict const x0 = geom.ProbLo();
  const CCTK_REAL *restrict const dx = geom.CellSize();

    cctkGH->cctk_origin_space[d] = x0[d];
    cctkGH->cctk_delta_space[d] = dx[d];

    cctkGH->cctk_origin_space[d] = NAN;
    cctkGH->cctk_delta_space[d] = NAN;

  // The refinement factor over the top level (coarsest) grid
  for (int d = 0; d < dim; ++d)
    cctkGH->cctk_levfac[d] = 1 << leveldata.level;
  // Offset between this level's and the coarsest level's origin as multiple
  // of the grid spacing

  const amrex::Geometry &geom = ghext->amrcore->Geom(0);
  const CCTK_REAL *restrict const global_x0 = geom.ProbLo();
  const CCTK_REAL *restrict const global_dx = geom.CellSize();

    cctkGH->cctk_levoff[d] = 1 - 2 * cctkGH->cctk_nghostzones[d];
    cctkGH->cctk_levoffdenom[d] = 2;

    // The refinement factor over the top level (coarsest) grid
    const int levfac = 1 << leveldata.level;
    cctkGH->cctk_levfac[d] = levfac;
    // Offset between this level's and the coarsest level's origin as multiple
    // of the grid spacing
    const int levoff = 1 - 2 * cctkGH->cctk_nghostzones[d];
    const int levoffdenom = 2;
    cctkGH->cctk_levoff[d] = levoff;
    cctkGH->cctk_levoffdenom[d] = levoffdenom;
    // Coordinates
    const CCTK_REAL origin_space =
        global_x0[d] +
        (1 - cctkGH->cctk_nghostzones[d] + CCTK_REAL(1) / 2) * global_dx[d];
    const CCTK_REAL delta_space = global_dx[d];
    cctkGH->cctk_delta_space[d] = delta_space / levfac;
    cctkGH->cctk_origin_space[d] =
        origin_space + cctkGH->cctk_delta_space[d] * levoff / levoffdenom;

  for (int d = 0; d < dim; ++d)

  for (int d = 0; d < dim; ++d) {

  for (int d = 0; d < dim; ++d)

  }
  for (int d = 0; d < dim; ++d) {
    cctkGH->cctk_origin_space[d] = NAN;
    cctkGH->cctk_delta_space[d] = NAN;
  }

// TODO: Have separate cctkGH for each level, each local box, and each tile

        x0[d] = cctkGH->cctk_origin_space[d];

        x0[d] = cctkGH->cctk_origin_space[d] - dx[d] / 2;