Write a multi level sliced plotfile

[dalarlla]

I am interested in being able to write effectively 2d slices of my 3d multilevel domain. I have seen some previous discussion of this and have tried to follow that approach but I am having issues.

Here is my implementation:

`void MyLevel::sliceWriter(const Output* out, int slice_dir, Real slice_coord) {

const char* dirs = "XYZ";
std::string dir(1, dirs[slice_dir]);

std::string prefix = domPtr->basePath.string()
    + "/" + out->name
    + "/pltslice_" + dir + "_"
    + Blast::formatString("%g", slice_coord) + "_";

Vector<std::unique_ptr<MultiFab>>  data_vec;
Vector<Geometry>                   geom_vec;
Vector<int>                        level_steps;
Vector<IntVect>                    ref_ratios;

Real time = state[Conserved_Type].curTime();
int  step = parent->levelSteps(0);
int  finest = parent->finestLevel();

for (int lev = 0; lev <= finest; lev++) {

    MyLevel& lev_slice = getLevel(lev);
    const Geometry& geom_full = lev_slice.geom;

    auto selected = createSelectedVariableMultiFab(out->fields, lev);

    auto sliced = createSlicedMultiFab(*selected, slice_dir, slice_coord, geom_full);

    RealBox real_slice = geom_full.ProbDomain();

    Real dx = geom_full.CellSize(slice_dir);
    real_slice.setLo(slice_dir, slice_coord - 0.5 * dx);
    real_slice.setHi(slice_dir, slice_coord + 0.5 * dx);

    Box domain_of_slice = geom_full.Domain();
    domain_of_slice.setRange(slice_dir, sliced->boxArray()[0].smallEnd(slice_dir));

    amrex::Geometry geom_of_slice(domain_of_slice, &real_slice);

    geom_vec.emplace_back(geom_of_slice);
    data_vec.emplace_back(std::move(sliced));

    if (lev < finest) {

        ref_ratios.push_back(parent->refRatio(lev));
        level_steps.push_back(step);
    }
}

std::string fname = amrex::Concatenate(prefix, step);

amrex::WriteMultiLevelPlotfile(
    fname,
    data_vec.size(),
    amrex::GetVecOfConstPtrs(data_vec),
    out->fields,
    geom_vec,
    time,
    level_steps,
    ref_ratios);

}

std::unique_ptramrex::MultiFab MyLevel::createSlicedMultiFab(
const amrex::MultiFab& mf, int slice_dir, Real slice_coord, Geometry levG) const {

auto sliced_mf = amrex::get_slice_data(slice_dir, slice_coord, mf, levG, 0, mf.nComp());

return sliced_mf;

}

std::unique_ptramrex::MultiFab MyLevel::createSelectedVariableMultiFab(
const std::vectorstd::string& var_names, int lev) const {

MyLevel& amr_level = static_cast<MyLevel&>(parent->getLevel(lev));

const amrex::MultiFab& ref_mf = amr_level.get_new_data(Conserved_Type);

int ncomp = var_names.size();
auto output_mf = std::make_unique<amrex::MultiFab>(ref_mf.boxArray(), ref_mf.DistributionMap(), ncomp, 0);

for (int i = 0; i < ncomp; ++i) {

    StateType stateType = varMap[var_names[i]].first;
    int comp = varMap[var_names[i]].second;

    const amrex::MultiFab& source_mf = amr_level.get_new_data(stateType);

    amrex::MultiFab::Copy(*output_mf, source_mf, comp, i, 1, 0);
}

return output_mf;

}
`

The main issue is that when I open the plotfile in paraview all the levels appear at once and overlap and clip into each other. Is there an obvious error in my method?