A mixed MoL–TMoL for the numerical solution of the 2D Richards’ equation in layered soils