I can recommend to you to use HP1. In this program, Hydrus-1D is coupled with the PHREEQC geochemical code (Parkhurst & Appelo 1999) to create a new comprehensive simulation tool, HP1 (acronym for HYDRUS1D-PHREEQC) (Jacques and Šimůnek 2005; Jacques et al. 2006; Šimůnek et al. 2006, 2008). This code contains modules simulating (1) transient water flow, (2) the transport of multiple components, (3) mixed equilibrium/kinetic biogeochemical reactions, and (4) heat transport in one-dimensional variably-saturated porous media (soils). HP1 is a significant expansion of the individual Hydrus-1D and PHREEQC programs by preserving most of their original features.