A matlab code was developed to calculate pressure heads and future concentrations given known boundary conditions for pressure heads and concentrations.