A discrete BGK approximation for strongly degenerate parabolic problems with boundary conditions