mf6.open_cbc does not return the right lower face face in a 2D(layer, x) model?
Checking the result of a 2D example, I'm getting 0 values in the lower face flow, while I have only one row defined.
I didn't test the _dis_indices
function for 2D arrays, it might not be returning the right numbers here ...