Describir: High-resolution shock-capturing numerical simulations of three-phase immiscible fluids from the unsaturated to the saturated zone