An upgraded version of a numerical solver is introduced to better capture the 3D interactions between surface water and groundwater. Built using open-source software, the custom solver adds new features to handle the complexity of real environments, including the representation of subsurface geology and the simulation of diverse dynamic processes, such as solute transport and heat transfer, in both domains. A test case and a full description of the novel features are provided in this paper.