Mass conservative numerical solutions of the head‐based Richards equation