INCLUDES HEADER: AMReX.hxx IN AMReX.hxx

#include <AMReX.hxx>

#include <driver.hxx>

    groupdata.numtimelevels = group.numtimelevels;

    // Allocate grid hierarchy
    groupdata.mfab =
        MultiFab(level.grids, level.dmap,
                 groupdata.numvars * groupdata.numtimelevels, ghost_size);

    // Allocate grid hierarchies
    groupdata.mfab.resize(group.numtimelevels);
    for (int tl = 0; tl < int(groupdata.mfab.size()); ++tl) {
      groupdata.mfab.at(tl) = make_unique<MultiFab>(
          level.grids, level.dmap, groupdata.numvars, ghost_size);
    }

#ifndef AMREX_HXX
#define AMREX_HXX

#ifndef DRIVER_HXX
#define DRIVER_HXX

      int numtimelevels;
      MultiFab mfab;

      vector<unique_ptr<MultiFab> > mfab; // [time level]

#endif // #ifndef AMREX_HXX

#endif // #ifndef DRIVER_HXX

#include <AMReX.hxx>

#include <driver.hxx>

  for (int gi = 0; gi < numgroups; ++gi)

  for (int gi = 0; gi < numgroups; ++gi) {

      if (groupdata.mfab.size() > 0) {
        const int tl = 0;

      string groupname = unique_ptr<char>(CCTK_GroupName(gi)).get();
      groupname = regex_replace(groupname, regex("::"), "-");
      for (auto &c : groupname)
        c = tolower(c);
      ostringstream buf;
      buf << "wavetoy/" << groupname;
      buf << ".rl" << setw(2) << setfill('0') << level.level;
      buf << ".it" << setw(6) << setfill('0') << cctk_iteration;
      string filename = buf.str();

        string groupname = unique_ptr<char>(CCTK_GroupName(gi)).get();
        groupname = regex_replace(groupname, regex("::"), "-");
        for (auto &c : groupname)
          c = tolower(c);
        ostringstream buf;
        buf << "wavetoy/" << groupname;
        buf << ".rl" << setw(2) << setfill('0') << level.level;
        buf << ".it" << setw(6) << setfill('0') << cctk_iteration;
        string filename = buf.str();

      Vector<string> varnames(groupdata.numvars * groupdata.numtimelevels);
      for (int tl = 0; tl < groupdata.numtimelevels; ++tl) {

        Vector<string> varnames(groupdata.numvars);

          varnames.at(tl * groupdata.numvars + n) = buf.str();

          varnames.at(n) = buf.str();

        // TODO: Output all levels into a single file
        // TODO: Output all groups into a single file
        WriteSingleLevelPlotfile(filename, *groupdata.mfab.at(tl), varnames,
                                 level.geom, cctk_time, cctk_iteration);

      // TODO: Output only current time level. This might require having one
      // mfab per time level
      // TODO: Output all levels into a single file
      // TODO: Output everything into a single file
      WriteSingleLevelPlotfile(filename, groupdata.mfab, varnames, level.geom,
                               cctk_time, cctk_iteration);

  }

#include <AMReX.hxx>

#include <driver.hxx>

#include <algorithm>

    const Array4<CCTK_REAL> &vars = groupdata.mfab.array(mfi);
    for (int tl = 0; tl < groupdata.numtimelevels; ++tl)
      for (int n = 0; n < groupdata.numvars; ++n)

    for (int tl = 0; tl < int(groupdata.mfab.size()); ++tl) {
      const Array4<CCTK_REAL> &vars = groupdata.mfab.at(tl)->array(mfi);
      for (int n = 0; n < groupdata.numvars; ++n) {

            vars.ptr(imin.x, imin.y, imin.z, tl * groupdata.numvars + n);

            vars.ptr(imin.x, imin.y, imin.z, n);
      }
    }

    for (int tl = 0; tl < groupdata.numtimelevels; ++tl)

    for (int tl = 0; tl < int(groupdata.mfab.size()); ++tl) {

    }

  for (auto &restrict level : ghext->levels)

  for (auto &restrict level : ghext->levels) {

      for (int tl = groupdata.numtimelevels - 1; tl > 0; --tl)
        MultiFab::Copy(groupdata.mfab, groupdata.mfab,
                       (tl - 1) * groupdata.numvars, tl * groupdata.numvars,
                       groupdata.numvars, ghost_size);

      const int ntls = groupdata.mfab.size();
      if (ntls > 1) {
        unique_ptr<MultiFab> tmp = move(groupdata.mfab.at(ntls - 1));
        for (int tl = ntls - 1; tl > 0; --tl)
          groupdata.mfab.at(tl) = move(groupdata.mfab.at(tl - 1));
        groupdata.mfab.at(0) = move(tmp);
      }

  }

        MultiFab &mfab = level.groupdata.at(0).mfab;

        MultiFab &mfab = *level.groupdata.at(0).mfab.at(0);

  for (auto &restrict level : ghext->levels)

  for (auto &restrict level : ghext->levels) {

      // we always sync all directions
      groupdata.mfab.FillBoundary(level.geom.periodicity());

      // We always sync all directions.
      // If there is more than one time level, then we don't sync the
      // oldest.
      int ntls = groupdata.mfab.size();
      int sync_tl = ntls > 1 ? ntls - 1 : ntls;
      for (int tl = 0; tl < sync_tl; ++tl)
        groupdata.mfab.at(tl)->FillBoundary(level.geom.periodicity());

  }

  for (int k = 0; k < cctk_lsh[2]; ++k)
    for (int j = 0; j < cctk_lsh[1]; ++j)

  for (int k = 0; k < cctk_lsh[2]; ++k) {
    for (int j = 0; j < cctk_lsh[1]; ++j) {

    }
  }

  for (int k = 0; k < cctk_lsh[2]; ++k)
    for (int j = 0; j < cctk_lsh[1]; ++j)

  for (int k = 0; k < cctk_lsh[2]; ++k) {
    for (int j = 0; j < cctk_lsh[1]; ++j) {

    }
  }

  for (int k = 0; k < cctk_lsh[2]; ++k)
    for (int j = 0; j < cctk_lsh[1]; ++j)

  for (int k = 0; k < cctk_lsh[2]; ++k) {
    for (int j = 0; j < cctk_lsh[1]; ++j) {

    }
  }