VanLeerPhi2 uses boundaries conditions
In VanLeerPhi2, we make decisions based on the boundaries decision, which we cannot do, if only because boundaries condition are configurable by users.
bool openZ = disk().physic().boundaries.height == BCT_OPEN;
.....
if ( openZ && (h == 1 || h == ni-2) ) {
lostByDisk += varq*volume[l];
}
By the way, shouldn't the mass lost at boundaries be computed by boundaries conditions ?