A flexible numerical approach to solving a generalized Richard's equation problem and some applications