An adaptive numerical method for the Richards equation with root growth