How do I set up a heat transport model in MODFLOW 6 using the GWE (energy) package?
A GWE model couples to a GWF flow simulation and adds packages specific to thermal energy. In FloPy you instantiate ModflowGwe with an EST package for thermal storage, a CND package for conduction and mechanical dispersion, and boundary conditions through CTP or the advanced packages (SFE for streams, LKE for lakes, MWE for multi-aquifer wells). A GwfGweExchange block wires the two models together so heads from GWF drive temperature in GWE.