~ ~ ~ PURPOSE ~ ~ ~ this subroutine calculates heat advection and dispersion for groundwater heat transport. Called from gwflow_lateral once per flow time step, after the Darcy head update. Uses cell_con%latl and cell_con%sat populated during the flow calculation in the same time step.