BOOLEAN output "Output 1D data for debugging" STEERABLE=recover
{
} no

SCHEDULE GROUP HydroToyAMReX_OutputGroup AT poststep
{
} "Output HydroToy quantities"

  SCHEDULE HydroToyAMReX_Pressure IN HydroToyAMReX_OutputGroup

if (output) {
  SCHEDULE GROUP HydroToyAMReX_OutputGroup AT poststep

    LANG: C
    READS: conserved(everywhere)
    WRITES: primitive(everywhere)
  } "Calculate pressure"

  } "Output HydroToy quantities"

  SCHEDULE HydroToyAMReX_Fluxes IN HydroToyAMReX_OutputGroup AFTER HydroToyAMReX_Pressure
  {
    LANG: C
    READS: conserved(everywhere)
    READS: primitive(everywhere)
    WRITES: flux_x(interior) flux_y(interior) flux_z(interior)
  } "Calculate the hydro fluxes"

    SCHEDULE HydroToyAMReX_Pressure IN HydroToyAMReX_OutputGroup
    {
      LANG: C
      READS: conserved(everywhere)
      WRITES: primitive(everywhere)
    } "Calculate pressure"
    SCHEDULE HydroToyAMReX_Fluxes IN HydroToyAMReX_OutputGroup AFTER HydroToyAMReX_Pressure
    {
      LANG: C
      READS: conserved(everywhere)
      READS: primitive(everywhere)
      WRITES: flux_x(interior) flux_y(interior) flux_z(interior)
    } "Calculate the hydro fluxes"

  SCHEDULE HydroToyAMReX_Output IN HydroToyAMReX_OutputGroup AFTER HydroToyAMReX_Fluxes
  {
    LANG: C
    READS: conserved(interior)
    READS: flux_x(interior) flux_y(interior) flux_z(interior)
    READS: AMReX::regrid_error(interior)
  } "Output grid data"

    SCHEDULE HydroToyAMReX_Output IN HydroToyAMReX_OutputGroup AFTER HydroToyAMReX_Fluxes
    {
      LANG: C
      READS: conserved(interior)
      READS: flux_x(interior) flux_y(interior) flux_z(interior)
      READS: AMReX::regrid_error(interior)
    } "Output grid data"
} # if output

  // const CCTK_REAL t = cctk_time;
  // const CCTK_REAL dt = CCTK_DELTA_TIME;

#if 0

  const Loop::GF3D<const CCTK_REAL, 1, 1, 1> regrid_error_(cctkGH,
                                                           regrid_error);

  Loop::loop_all<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
    if (p.j == 0 && p.k == 0) {

  int levfac = cctk_levfac[0];
  int lev = 0;
  while (levfac > 1) {
    levfac >>= 1;
    lev += 1;
  }

    const int oldprec = cout.precision(17);
    cout << "iteration " << cctk_iteration << " level " << lev << ":\n";
    Loop::loop_all<1, 1, 1>(cctkGH, [&](const Loop::PointDesc &p) {
      if (p.j == 0 && p.k == 0) {

             << p.y << "," << p.z << ") rho=" << rho_(p.I)
             << " momx=" << momx_(p.I) << " etot=" << etot_(p.I) << "\n";
      }

             << p.y << "," << p.z << ") rho="
             << rho_(p.I)
             // << " momx=" << momx_(p.I) << " etot=" << etot_(p.I)
             // << " err=" << regrid_error_(p.I)
             << "\n";

#pragma omp critical(HydroToyAMReX_Output)
      {

             << p.y << "," << p.z << ") fxrho=" << fxrho_(p.I)
             << " fxmomx=" << fxmomx_(p.I) << " fxetot=" << fxetot_(p.I)

             << p.y << "," << p.z << ") fxrho="
             << fxrho_(p.I)
             // << " fxmomx=" << fxmomx_(p.I) << " fxetot=" << fxetot_(p.I)

    }

    cout.precision(oldprec);
  }

#endif