MODFLOW-ORG · emorway-usgs · Jun 3, 2026 · Jun 3, 2026 · Jun 16, 2026 · Jun 17, 2026
diff --git a/autotest/test_gwe_sfe_abc00.py b/autotest/test_gwe_sfe_abc00.py
diff --git a/autotest/test_gwe_sfe_abc01.py b/autotest/test_gwe_sfe_abc01.py
diff --git a/autotest/test_gwe_sfe_abc02.py b/autotest/test_gwe_sfe_abc02.py
diff --git a/autotest/test_gwe_sfe_abc03.py b/autotest/test_gwe_sfe_abc03.py
diff --git a/autotest/test_gwe_sfe_lhf00.py b/autotest/test_gwe_sfe_lhf00.py
diff --git a/autotest/test_gwe_sfe_lwr00.py b/autotest/test_gwe_sfe_lwr00.py
diff --git a/autotest/test_gwe_sfe_shf00.py b/autotest/test_gwe_sfe_shf00.py
diff --git a/autotest/test_gwe_sfe_shf01.py b/autotest/test_gwe_sfe_shf01.py
diff --git a/autotest/test_gwe_sfe_swr00.py b/autotest/test_gwe_sfe_swr00.py
diff --git a/doc/Common/gwe-sfeobs.tex b/doc/Common/gwe-sfeobs.tex
@@ -14,4 +14,9 @@
 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. \\
-SFE & strmbd-cond & rno or boundname & -- & Amount of heat conductively exchanged with the streambed material. 
+SFE & strmbd-cond & rno or boundname & -- & Amount of heat conductively exchanged with the streambed material. \\
+SFE & swr & rno or boundname & -- & Amount of energy entering the stream reach from shortwave radiation. \\
+SFE & lwr & rno or boundname & -- & Net amount of longwave radiation energy exchanged with the stream reach. \\
+SFE & lhf & rno or boundname & -- & Amount of energy lost from the stream reach owing to evaporative cooling. \\
+SFE & shf & rno or boundname & -- & Amount of sensible heat exchanged with the atmosphere. \\
+SFE & surfevap & rno & -- & Evaporation rate for a given stream reach. \\
diff --git a/doc/mf6io/gwe/abc.tex b/doc/mf6io/gwe/abc.tex
@@ -0,0 +1,23 @@
+Input to the Atmospheric Boundary Condition (ABC) Package is read from the file that is specified in the ABC6 record of the OPTIONS block in the SFE package.  There can only be one ABC utility activated for each instance of SFE.  The ABC utility cannot be used without a corresponding SFE Package.
+
+The ABC utility does not have a dimensions block; instead, dimensions for the ABC utility are set using the dimensions from the corresponding SFE Package.  For example, an SFR Package corresponding to an instance of SFE requires specification of the number of reaches (NREACHES).  ABC sets the number of reaches equal to NREACHES.  Therefore, the PERIOD block below cannot contain boundary information for a reach with an index greater than NREACHES.
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\vspace{5mm}
+
+\noindent \textit{FOR EACH SIMULATION}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-abc-options.dat}
+\vspace{5mm}
+\noindent \textit{FOR ANY STRESS PERIOD}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-abc-period.dat}
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/utl-abc-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/utl-abc-example.dat}
diff --git a/doc/mf6io/gwe/gwe.tex b/doc/mf6io/gwe/gwe.tex
@@ -129,6 +129,10 @@ \subsection{Energy Source Loading (ESL) Package}
 \subsection{Streamflow Energy Transport (SFE) Package}
 \input{gwe/sfe}
 
+\newpage
+\subsection{Atmospheric Boundary Conditions (ABC) Package}
+\input{gwe/abc}
+
 \newpage
 \subsection{Lake Energy Transport (LKE) Package}
 \input{gwe/lke}

diff --git a/doc/mf6io/gwe/shf.tex b/doc/mf6io/gwe/shf.tex
@@ -0,0 +1,21 @@
+Input to the Sensible Heat Flux (SHF) Package is read from the file that is specified in the SHF6 record of the OPTIONS block in the SFE package.
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\vspace{5mm}
+
+\noindent \textit{FOR EACH SIMULATION}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-shf-options.dat}
+\vspace{5mm}
+\noindent \textit{FOR ANY STRESS PERIOD}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-shf-period.dat}
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/utl-shf-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/utl-shf-example.dat}
diff --git a/doc/mf6io/gwe/swr.tex b/doc/mf6io/gwe/swr.tex
@@ -0,0 +1,21 @@
+Input to the Shortwave Radiation (SWR) Package is read from the file that is specified in the SWR6 record of the OPTIONS block in the SFE package.
+
+\vspace{5mm}
+\subsubsection{Structure of Blocks}
+\vspace{5mm}
+
+\noindent \textit{FOR EACH SIMULATION}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-swr-options.dat}
+\vspace{5mm}
+\noindent \textit{FOR ANY STRESS PERIOD}
+\lstinputlisting[style=blockdefinition]{./mf6ivar/tex/utl-swr-period.dat}
+
+\vspace{5mm}
+\subsubsection{Explanation of Variables}
+\begin{description}
+\input{./mf6ivar/tex/utl-swr-desc.tex}
+\end{description}
+
+\vspace{5mm}
+\subsubsection{Example Input File}
+\lstinputlisting[style=inputfile]{./mf6ivar/examples/utl-swr-example.dat}
diff --git a/doc/mf6io/mf6ivar/dfn/gwe-sfe.dfn b/doc/mf6io/mf6ivar/dfn/gwe-sfe.dfn
@@ -255,6 +255,37 @@ optional false
 longname obs6 input filename
 description REPLACE obs6_filename {'{#1}': 'SFE', '{#2}': '\\ref{table:gwe-obstypetable}'}
 
+block options
+name abc_filerecord
+type record abc6 filein abc6_filename
+shape
+reader urword
+tagged true
+optional true
+longname
+description
+
+block options
+name abc6
+type keyword
+shape
+in_record true
+reader urword
+tagged true
+optional false
+longname abc keyword
+description keyword to specify that record corresponds to the atmospheric boundary condition (ABC) file.  The behavior of the ABC utility and a description of the input file is provided separately.
+
+block options
+name abc6_filename
+type string
+preserve_case true
+in_record true
+reader urword
+optional false
+tagged false
+longname file name of ABC information
+description defines a sensible heat flux (ABC) input file.  Records in the ABC file are used to calculate the flux of thermal energy in or out of a stream reach resulting from sensible heat exchange.
 
 # --------------------- gwe sfe packagedata ---------------------