MODFLOW-USGS · emorway-usgs · Jul 26, 2022 · Jul 27, 2022 · Jul 28, 2022 · Jul 28, 2022
diff --git a/autotest/test_gwe_drycell_conduction0.py b/autotest/test_gwe_drycell_conduction0.py
diff --git a/autotest/test_gwe_drycell_conduction1.py b/autotest/test_gwe_drycell_conduction1.py
diff --git a/autotest/test_gwe_drycell_conduction2.py b/autotest/test_gwe_drycell_conduction2.py
diff --git a/autotest/test_gwe_dsp.py b/autotest/test_gwe_dsp.py
diff --git a/autotest/test_gwe_lke_conduction.py b/autotest/test_gwe_lke_conduction.py
diff --git a/autotest/test_gwe_mwe_conduction.py b/autotest/test_gwe_mwe_conduction.py
diff --git a/autotest/test_gwe_sfe_strmbedcond.py b/autotest/test_gwe_sfe_strmbedcond.py
diff --git a/autotest/test_gwe_stallman.py b/autotest/test_gwe_stallman.py
diff --git a/autotest/test_gwe_uze00.py b/autotest/test_gwe_uze00.py
diff --git a/autotest/test_gwe_uze00_flux.py b/autotest/test_gwe_uze00_flux.py
diff --git a/doc/Common/gwe-cncobs.tex b/doc/Common/gwe-cncobs.tex
@@ -0,0 +1 @@
+CNC & cnc & cellid or boundname & -- & Energy flow between the groundwater system and a constant-temperature boundary or a group of cells with constant-temperature boundaries.
diff --git a/doc/Common/gwe-gweobs.tex b/doc/Common/gwe-gweobs.tex
@@ -0,0 +1 @@
+GWE-GWE & flow-ja-face & exchange number or boundname & -- & Mass flow between model 1 and model 2 for a specified exchange (which is the consecutive exchange number listed in the EXCHANGEDATA block), or the sum of these exchange flows by boundname if boundname is specified.
diff --git a/doc/Common/gwe-lkeobs.tex b/doc/Common/gwe-lkeobs.tex
@@ -0,0 +1,18 @@
+% general APT observations
+LKE & temperature & lakeno or boundname & -- & Lake temperature. If boundname is specified, boundname must be unique for each lake. \\
+LKE & flow-ja-face & lakeno or boundname & lakeno or -- & Energy flow between two lakes connected by an outlet.  If more than one outlet is used to connect the same two lakes, then the energy flow for only the first outlet can be observed.  If a boundname is specified for ID1, then the result is the total energy flow for all outlets for a lake. If a boundname is specified for ID1 then ID2 is not used.\\
+LKE & storage & lakeno or boundname & -- & Simulated energy storage flow rate for a lake or group of lakes. \\
+LKE & constant & lakeno or boundname & -- & Simulated energy constant-flow rate for a lake or group of lakes. \\
+LKE & from-mvr & lakeno or boundname & -- & Simulated energy inflow into a lake or group of lakes from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
+LKE & to-mvr & outletno or boundname & -- & Energy outflow from a lake outlet, a lake, or a group of lakes that is available for the MVR package. If boundname is not specified for ID, then the outflow available for the MVR package from a specific lake outlet is observed. In this case, ID is the outlet number, which must be between 1 and NOUTLETS. \\
+LKE & lke & lakeno or boundname & \texttt{iconn} or -- & Energy flow rate for a lake or group of lakes and its aquifer connection(s). If boundname is not specified for ID, then the simulated lake-aquifer flow rate at a specific lake connection is observed. In this case, ID2 must be specified and is the connection number \texttt{iconn} for lake \texttt{lakeno}. \\
+
+%observations specific to the lake package
+% rainfall evaporation runoff ext-inflow withdrawal outflow
+LKE & rainfall & lakeno or boundname & -- & Rainfall rate applied to a lake or group of lakes multiplied by the rainfall temperature. \\
+LKE & evaporation & lakeno or boundname & -- & Simulated evaporation rate from a lake or group of lakes multiplied by the latent heat of evaporation for determining the energy lost from a lake. \\
+LKE & runoff & lakeno or boundname & -- & Runoff rate applied to a lake or group of lakes multiplied by the runoff temperature. \\
+LKE & ext-inflow & lakeno or boundname & -- & Energy inflow into a lake or group of lakes calculated as the external inflow rate multiplied by the inflow temperature. \\
+LKE & withdrawal & lakeno or boundname & -- & Specified withdrawal rate from a lake or group of lakes multiplied by the simulated lake temperature. \\
+LKE & ext-outflow & lakeno or boundname & -- & External outflow from a lake or a group of lakes, through their outlets, to an external boundary.  If the water mover is active, the reported ext-outflow value plus the rate to mover is equal to the total outlet outflow.
+
diff --git a/doc/Common/gwe-mweobs.tex b/doc/Common/gwe-mweobs.tex
@@ -0,0 +1,13 @@
+% general APT observations
+MWE & temperature & mawno or boundname & -- & Well temperature. If boundname is specified, boundname must be unique for each well. \\
+%flowjaface not included
+MWE & storage & mawno or boundname & -- & Simulated energy storage flow rate for a well or group of wells. \\
+MWE & constant & mawno or boundname & -- & Simulated energy constant-flow rate for a well or group of wells. \\
+MWE & from-mvr & mawno or boundname & -- & Simulated energy inflow into a well or group of wells from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
+MWE & mwe & mawno or boundname & \texttt{iconn} or -- & Energy flow rate for a well or group of wells and its aquifer connection(s). If boundname is not specified for ID, then the simulated well-aquifer flow rate at a specific well connection is observed. In this case, ID2 must be specified and is the connection number \texttt{iconn} for well \texttt{mawno}. \\
+
+% observations specific to the mwe package
+MWE & rate & mawno or boundname & -- & Simulated energy flow rate for a well or group of wells. \\
+MWE & fw-rate & mawno or boundname & -- & Simulated energy flow rate for a flowing well or group of flowing wells. \\
+MWE & rate-to-mvr & well or boundname & -- & Simulated energy flow rate that is sent to the MVE Package for a well or group of wells.\\
+MWE & fw-rate-to-mvr & well or boundname & -- & Simulated energy flow rate that is sent to the MVE Package from a flowing well or group of flowing wells. \\
diff --git a/doc/Common/gwe-obs.tex b/doc/Common/gwe-obs.tex
@@ -0,0 +1,2 @@
+GWE & temperature & cellid & -- & Temperature at a specified cell. \\
+GWE & flow-ja-face & cellid & cellid & Energy flow in dimensions of energy per time between two adjacent cells.  The energy flow rate includes the contributions from both convection (advection) and conduction (including mechanical dispersion) if those packages are active
diff --git a/doc/Common/gwe-sfeobs.tex b/doc/Common/gwe-sfeobs.tex
@@ -0,0 +1,16 @@
+% general APT observations
+SFE & temperature & rno or boundname & -- & Reach temperature. If boundname is specified, boundname must be unique for each reach. \\
+SFE & flow-ja-face & rno or boundname & rno or -- & Energy flow between two reaches.  If a boundname is specified for ID1, then the result is the total energy flow for all reaches. If a boundname is specified for ID1 then ID2 is not used.\\
+SFE & storage & rno or boundname & -- & Simulated energy storage flow rate for a reach or group of reaches. \\
+SFE & constant & rno or boundname & -- & Simulated energy constant-flow rate for a reach or group of reaches. \\
+SFE & from-mvr & rno or boundname & -- & Simulated energy inflow into a reach or group of reaches from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
+SFE & to-mvr & rno or boundname & -- & Energy outflow from a reach, or a group of reaches that is available for the MVR package. If boundname is not specified for ID, then the outflow available for the MVR package from a specific reach is observed. \\
+SFE & sfe & rno or boundname & -- & Energy flow rate for a reach or group of reaches and its aquifer connection(s). \\
+
+%observations specific to the stream energy transport package
+% rainfall evaporation runoff ext-inflow withdrawal outflow
+SFE & rainfall & rno or boundname & -- & Rainfall rate applied to a reach or group of reaches multiplied by the rainfall temperature. \\
+SFE & evaporation & rno or boundname & -- & Simulated evaporation rate from a reach or group of reaches multiplied by the latent heat of vaporization for determining the amount of energy lost from a reach. \\
+SFE & runoff & rno or boundname & -- & Runoff rate applied to a reach or group of reaches multiplied by the runoff temperature. \\
+SFE & ext-inflow & rno or boundname & -- & Energy inflow into a reach or group of reaches calculated as the external inflow rate multiplied by the inflow temperature. \\
+SFE & ext-outflow & rno or boundname & -- & External outflow from a reach or group of reaches to an external boundary. If boundname is not specified for ID, then the external outflow from a specific reach is observed. In this case, ID is the reach rno.
diff --git a/doc/Common/gwe-srcobs.tex b/doc/Common/gwe-srcobs.tex
@@ -0,0 +1 @@
+SRC & src & cellid or boundname & -- & Energy source loading rate between the groundwater system and a energy source loading boundary or a group of  boundaries.
diff --git a/doc/Common/gwe-uzeobs.tex b/doc/Common/gwe-uzeobs.tex
@@ -0,0 +1,14 @@
+% general APT observations
+UZE & temperature & uzeno or boundname & -- & uze cell temperature. If boundname is specified, boundname must be unique for each uze cell. \\
+UZE & flow-ja-face & uzeno or boundname & uzeno or -- & Energy flow between two uze cells.  If a boundname is specified for ID1, then the result is the total energy flow for all uze cells. If a boundname is specified for ID1 then ID2 is not used.\\
+UZE & storage & uzeno or boundname & -- & Simulated energy storage flow rate for a uze cell or group of uze cells. \\
+UZE & constant & uzeno or boundname & -- & Simulated energy constant-flow rate for a uze cell or a group of uze cells. \\
+UZE & from-mvr & uzeno or boundname & -- & Simulated energy inflow into a uze cell or group of uze cells from the MVE package. Energy inflow is calculated as the product of provider temperature and the mover flow rate. \\
+UZE & uze & uzeno or boundname & -- & Energy flow rate for a uze cell or group of uze cells and its aquifer connection(s). \\
+
+%observations specific to the uze package
+% infiltration rej-inf uzet rej-inf-to-mvr
+UZE & infiltration & uzeno or boundname & -- & Infiltration rate applied to a uze cell or group of uze cells multiplied by the infiltration temperature. \\
+UZE & rej-inf & uzeno or boundname & -- & Rejected infiltration rate applied to a uze cell or group of uze cells multiplied by the infiltration temperature. \\
+UZE & uzet & uzeno or boundname & -- & Unsaturated zone evapotranspiration rate applied to a uze cell or group of uze cells multiplied by the uze cell temperature. \\
+UZE & rej-inf-to-mvr & uzeno or boundname & -- & Rejected infiltration rate applied to a uze cell or group of uze cells multiplied by the infiltration temperature that is sent to the mover package. \\
diff --git a/doc/MODFLOW6References.bib b/doc/MODFLOW6References.bib
@@ -8,6 +8,28 @@
 
 
 
+@article{mazheng2010,
+  author = {Ma, Rui and Zheng, Chunmiao},
+  title = {Effects of Density and Viscosity in Modeling Heat as a Groundwater Tracer},
+  journal = {Groundwater},
+  year = {2010},
+  volume = {48},
+  number = {3},
+  pages = {380--389},
+  doi = {10.1111/j.1745-6584.2009.00660.x},
+  }
+
+@article{hechtmendez,
+  author = {Hecht-Mendez, J. and Molina-Giraldo, N. and Blum, P. and Bayer, P.},
+  title = {Evaluating MT3DMS for Heat Transport Simulation of Closed Geothermal Systems},
+  journal = {Groundwater},
+  year = {2010},
+  volume = {48},
+  number = {5},
+  pages = {741--756},
+  doi = {10.1111/j.1745-6584.2010.00678.x},
+  }
+
 @article{morway2021,
   author = {Morway, Eric D. and Langevin, Christian D. and Hughes, Joseph D.},
   title = {Use of the {MODFLOW 6} Water Mover Package to Represent Natural and Managed Hydrologic Connections},

diff --git a/doc/mf6io/body.tex b/doc/mf6io/body.tex
@@ -41,6 +41,11 @@
 \SECTION{Groundwater Transport (GWT) Model Input}
 \input{gwt/gwt.tex}
 
+%GWE Model Input Instructions
+\newpage
+\SECTION{Groundwater Energy Transport (GWE) Model Input}
+\input{gwe/gwe.tex}
+
 %Sparse Matrix Solution (IMS)
 \newpage
 \SECTION{Iterative Model Solution}

diff --git a/doc/mf6io/gwe/adv.tex b/doc/mf6io/gwe/adv.tex
@@ -0,0 +1,16 @@
+Advection (ADV) Package information is read from the file that is specified by ``ADV6'' as the file type.  Only one ADV Package can be specified for a GWE model. 
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-adv-options.dat}
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/gwe-adv-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-adv-example.dat}
+
diff --git a/doc/mf6io/gwe/cnc.tex b/doc/mf6io/gwe/cnc.tex
@@ -0,0 +1,47 @@
+Constant Temperature (CNC) Package information is read from the file that is specified by ``CNC6'' as the file type.  Any number of CNC Packages can be specified for a single GWE model, but the same cell cannot be designated as a constant temperature by more than one CNC entry. 
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\vspace{5mm}
+
+\noindent \textit{FOR EACH SIMULATION}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-cnc-options.dat}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-cnc-dimensions.dat}
+\vspace{5mm}
+\noindent \textit{FOR ANY STRESS PERIOD}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-cnc-period.dat}
+\packageperioddescription
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/gwe-cnc-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-cnc-example.dat}
+
+\vspace{5mm}
+\subsubsection{Available observation types}
+CNC Package observations are limited to the simulated constant temperature energy flow rate (\texttt{cnc}). The data required for the CNC Package observation type is defined in table~\ref{table:gwe-cncobstype}. Negative and positive values for an observation represent a loss from and gain to the GWE model, respectively.
+
+\begin{longtable}{p{2cm} p{2.75cm} p{2cm} p{1.25cm} p{7cm}}
+\caption{Available CNC Package observation types} \tabularnewline
+
+\hline
+\hline
+\textbf{Model} & \textbf{Observation type} & \textbf{ID} & \textbf{ID2} & \textbf{Description} \\
+\hline
+\endhead
+
+\hline
+\endfoot
+
+\input{../Common/gwe-cncobs.tex}
+\label{table:gwe-cncobstype}
+\end{longtable}
+
+\vspace{5mm}
+\subsubsection{Example Observation Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-cnc-example-obs.dat}
diff --git a/doc/mf6io/gwe/dsp.tex b/doc/mf6io/gwe/dsp.tex
@@ -0,0 +1,17 @@
+Dispersion (DSP) Package information is read from the file that is specified by ``DSP6'' as the file type.  Only one DSP Package can be specified for a GWE model.  The DSP Package is based on the mathematical formulation presented for the XT3D option of the NPF Package available to represent full three-dimensional anisotropy in groundwater flow.  XT3D can be computationally expensive and can be turned off to use a simplified and approximate form of the dispersion equations.  For most problems, however, XT3D will be required to accurately represent dispersion.
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-dsp-options.dat}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/gwe-dsp-griddata.dat}
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/gwe-dsp-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwe-dsp-example.dat}
+