A numerical model for the movement of H2O, H218O, and 2HHO in the unsaturated zone