"AMReX: Domain size mismatch when using BoxArray and MFIter in 2D simulation"

[Afraid31]

I would like to set up a 2D fluid domain where nx is not equal to ny. Here's how I define the domain:
IntVect dom_lo((0, 0)); IntVect dom_hi((31, 15)); Box domain(dom_lo, dom_hi); ba.define(domain); ba.maxSize(16);
However, I found that when I loop through the domain using the following method, the simulation appears to run over a 16×16 domain instead of the intended 32×16 one. Here's the loop:
`for (MFIter mfi(macro); mfi.isValid(); ++mfi)
{
const Box& bx = mfi.validbox();
const Array4& fi_ = fi.array(mfi);

amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j)
{
    ...
});

}
`
Why is the loop running over a 16×16 domain instead of 32×16? Is there something I'm misunderstanding in how the domain or BoxArray is defined?

[WeiqunZhang]

ba.maxSize(16) chopped into two Boxes. Each Box has 16x16 cells. https://amrex-codes.github.io/amrex/docs_html/Basics.html#boxarray might be helpful.