Degenerate parabolic variational inequalities with convection are solved bymeans of a combined relaxation method and method of characteristics. Themathematical problem is motivated by Richard's equation, modelling theunsaturated – saturated flow in porous media. By means of the relaxationmethod we control the degeneracy. The dominance of the convection iscontrolled by the method of characteristics.