Dennis Gläser authored Jul 05, 2016 and Timo Koch committed Nov 23, 2016

the boundary treatment is now performed after all the interior
fluxes have been calculated. This corresponds more to the previous
implementation and is assumend to be more correct for mixed BCs.