Provide implementation of the McKee, Wall and Wilson ADI with 2nd order upwinding.