PRIVATE:

CCTK_REAL xmin "Domain size"
{
  *:* :: ""
} -1.0
CCTK_REAL xmax "Domain size"
{
  *:* :: ""
} 1.0
CCTK_REAL ymin "Domain size"
{
  *:* :: ""
} -1.0
CCTK_REAL ymax "Domain size"
{
  *:* :: ""
} 1.0
CCTK_REAL zmin "Domain size"
{
  *:* :: ""
} -1.0
CCTK_REAL zmax "Domain size"
{
  *:* :: ""
} 1.0
CCTK_INT ncells_x "Number of grid cells"
{
  1:*
} 100
CCTK_INT ncells_y "Number of grid cells"
{
  1:*
} 100
CCTK_INT ncells_z "Number of grid cells"
{
  1:*
} 100
CCTK_INT ghost_size "Number of ghost zones"
{
  0:* :: ""
} 1
CCTK_REAL dtfac "The standard timestep condition dt = dtfac*min(delta_space)"
{
  *:* :: ""
} 0.5 

  // Parameters that should not be here
  const int ncells = 100;    // grid size
  const int nghostzones = 1; // number of ghost zones

  DECLARE_CCTK_PARAMETERS;

  IntVect dom_hi(
      AMREX_D_DECL(ghext->ncells - 1, ghext->ncells - 1, ghext->ncells - 1));

  IntVect dom_hi(AMREX_D_DECL(ncells_x - 1, ncells_y - 1, ncells_z - 1));

  RealBox real_box({AMREX_D_DECL(-1.0, -1.0, -1.0)},
                   {AMREX_D_DECL(1.0, 1.0, 1.0)});

  RealBox real_box({AMREX_D_DECL(xmin, ymin, zmin)},
                   {AMREX_D_DECL(xmax, ymax, zmax)});

    groupdata.mfab =
        MultiFab(ghext->ba, dm, groupdata.numvars * groupdata.numtimelevels,
                 ghext->nghostzones);

    groupdata.mfab = MultiFab(
        ghext->ba, dm, groupdata.numvars * groupdata.numtimelevels, ghost_size);

  DECLARE_CCTK_PARAMETERS;

  int ncells[dim] = {ncells_x, ncells_y, ncells_z};
  CCTK_REAL x0[dim] = {xmin, ymin, zmin};
  CCTK_REAL x1[dim] = {xmax, ymax, zmax};
  CCTK_REAL dx[dim];

    cctkGH->cctk_origin_space[d] = 0.0;

    dx[d] = (x1[d] - x0[d]) / ncells[d];
  CCTK_REAL mindx = 1.0 / 0.0;
  for (int d = 0; d < dim; ++d)
    mindx = fmin(mindx, dx[d]);
  for (int d = 0; d < dim; ++d)
    cctkGH->cctk_origin_space[d] = x0[d];

    cctkGH->cctk_delta_space[d] = 1.0;

    cctkGH->cctk_delta_space[d] = dx[d];
  CCTK_VINFO("x0=[%g,%g,%g]", cctkGH->cctk_origin_space[0],
             cctkGH->cctk_origin_space[1], cctkGH->cctk_origin_space[2]);
  CCTK_VINFO("dx=[%g,%g,%g]", cctkGH->cctk_delta_space[0],
             cctkGH->cctk_delta_space[1], cctkGH->cctk_delta_space[2]);
  // Initialize time stepping

  cctkGH->cctk_delta_time = 0.0;

  cctkGH->cctk_delta_time = dtfac * mindx;
  CCTK_VINFO("t=%g", cctkGH->cctk_time);
  CCTK_VINFO("dt=%g", cctkGH->cctk_delta_time);

  DECLARE_CCTK_PARAMETERS;

  // TODO: Get this from mfab

    cctkGH->cctk_nghostzones[d] = ghext->nghostzones;

    cctkGH->cctk_nghostzones[d] = ghost_size;

  DECLARE_CCTK_PARAMETERS;

  // TODO: Get this from mfab
  int ncells[dim] = {ncells_x, ncells_y, ncells_z};

    cctkGH->cctk_gsh[d] = ghext->ncells;

    cctkGH->cctk_gsh[d] = ncells[d];

  DECLARE_CCTK_PARAMETERS;

  // TODO: Get ghost_size from mfab

                     groupdata.numvars, ghext->nghostzones);

                     groupdata.numvars, ghost_size);

SCHEDULE WaveToyAMReX_Setup AT basegrid
{
  LANG: C
  OPTIONS: level
} "Set up time"

extern "C" void WaveToyAMReX_Setup(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS;
  DECLARE_CCTK_PARAMETERS;
  const CCTK_REAL *restrict const x0 = ghext->geom.ProbLo();
  const CCTK_REAL *restrict const x1 = ghext->geom.ProbHi();
  for (int d = 0; d < dim; ++d) {
    CCTK_REAL dx = (x1[d] - x0[d]) / ghext->ncells;
    cctkGH->cctk_origin_space[d] = x0[d] + 0.5 * dx;
    cctkGH->cctk_delta_space[d] = dx;
  }
  CCTK_VINFO("x0=[%g,%g,%g]", cctkGH->cctk_origin_space[0],
             cctkGH->cctk_origin_space[1], cctkGH->cctk_origin_space[2]);
  CCTK_VINFO("dx=[%g,%g,%g]", cctkGH->cctk_delta_space[0],
             cctkGH->cctk_delta_space[1], cctkGH->cctk_delta_space[2]);

  CCTK_REAL mindx = 1.0 / 0.0;
  for (int d = 0; d < dim; ++d)
    mindx = fmin(mindx, cctkGH->cctk_delta_space[d]);
  cctkGH->cctk_time = 0.0;
  cctkGH->cctk_delta_time = 0.5 * mindx;
  CCTK_VINFO("t=%g", cctkGH->cctk_time);
  CCTK_VINFO("dt=%g", cctkGH->cctk_delta_time);
}