~ ~ ~ PURPOSE ~ ~ ~ this subroutine calculates solute advection, dispersion, chemical reactions, and sorption for groundwater solute transport. Called from gwflow_lateral once per flow time step. Contains its own sub-timestep loop (num_ts_transport steps). Uses cell_con%latl and cell_con%sat populated during the Darcy flow calculation in the same time step.