MODFLOW-ORG · emorway-usgs · Aug 22, 2025 · Mar 6, 2025 · Mar 13, 2025 · Apr 24, 2025
diff --git a/autotest/test_gwe_sfe_abc01.py b/autotest/test_gwe_sfe_abc01.py
diff --git a/autotest/test_gwe_sfe_abc_lhf00.py b/autotest/test_gwe_sfe_abc_lhf00.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,6 @@
 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 & shf & rno or boundname & -- & Amount of sensible heat exchanged with the atmosphere. \\
+SFE & swr & rno or boundname & -- & Amount of shortwave radiation on stream reach.
diff --git a/doc/mf6io/gwe/gwe.tex b/doc/mf6io/gwe/gwe.tex
@@ -139,6 +139,14 @@ \subsection{Multi-Aquifer Well Energy Transport (MWE) Package}
 \subsection{Unsaturated-Zone Energy Transport (UZE) Package}
 \input{gwe/uze}
 
+\newpage
+\subsection{Sensible Heat Flux (SHF) Package}
+\input{gwe/shf}
+
+\newpage
+\subsection{Shortwave Radiation (SWR) Package}
+\input{gwe/swr}
+
 \newpage
 \subsection{Flow Model Interface (FMI) Package}
 \input{gwe/fmi}

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
@@ -254,6 +254,37 @@ optional false
 longname obs6 input filename
 description REPLACE obs6_filename {'{#1}': 'SFE'}
 
+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 ---------------------
 

diff --git a/doc/mf6io/mf6ivar/dfn/utl-abc.dfn b/doc/mf6io/mf6ivar/dfn/utl-abc.dfn
@@ -0,0 +1,215 @@
+# --------------------- gwe abc options ---------------------
+# flopy subpackage abc_filerecord abc abc_perioddata perioddata
+# flopy parent_name_type parent_package MFPackage
+
+block options
+name print_input
+type keyword
+reader urword
+optional true
+longname print input to listing file
+description REPLACE print_input {'{#1}': 'atmospheric boundary condition'}
+
+block options
+name ts_filerecord
+type record ts6 filein ts6_filename
+shape
+reader urword
+tagged true
+optional true
+longname
+description
+
+block options
+name ts6
+type keyword
+shape
+in_record true
+reader urword
+tagged true
+optional false
+longname head keyword
+description keyword to specify that record corresponds to a time-series file.
+
+block options
+name filein
+type keyword
+shape
+in_record true
+reader urword
+tagged true
+optional false
+longname file keyword
+description keyword to specify that an input filename is expected next.
+
+block options
+name ts6_filename
+type string
+preserve_case true
+in_record true
+reader urword
+optional false
+tagged false
+longname file name of time series information
+description REPLACE timeseriesfile {}
+
+block options
+name density_air
+type double precision
+reader urword
+optional true
+longname density of water
+description density of the atmospheric air used by calculations related to sensible heat flux.  This value is set to 1.225 kg/m3 if no overriding value is specified.  A user-specified value should be provided for models that use units other than kilograms and meters or if it is necessary to use a value other than the default.
+default_value 1.225
+mf6internal rhoa
+
+block options
+name heat_capacity_air
+type double precision
+reader urword
+optional true
+longname heat capacity of air
+description mass-based heat capacity of the air used by calculations related to sensible heat flux.  This value is set to 717.0 J/kg/C if no overriding value is specified.  A user-specified value should be provided for models that use units other than joules, kilograms, and degrees celsius or if it is necessary to use a value other than the default.
+default_value 717.0
+mf6internal cpa
+
+block options
+name drag_coefficient
+type double precision
+reader urword
+optional true
+longname drag coefficient
+description drag coefficient used by calculations related to sensible heat flux.  This value is set to 0.002 (unitless) if no overriding value is specified.  A user-specified value should be provided for models that require a value other than the default.
+default_value 0.0002
+mf6internal cd
+
+block options
+name wind_func_slope
+type double precision
+reader urword
+optional true
+longname wind function slope
+description wind function slope used by calculations related to the latent heat flux.  This value is set to 1.383E-08 (1/mbar) if no overriding value is specified.  A user-specified value should be provided for models that require a value other than the default.
+default_value 1.383E-08
+mf6internal wf_slope
+
+block options
+name wind_func_int
+type double precision
+reader urword
+optional true
+longname wind function intercept
+description wind function intercept used by calculations related to the latent heat flux.  This value is set to 3.445E-09 (m/s) if no overriding value is specified.  A user-specified value should be provided for models that require a value other than the default.
+default_value 3.445E-09
+mf6internal wf_int
+
+# --------------------- gwe abc period ---------------------
+
+block period
+name iper
+type integer
+block_variable True
+in_record true
+tagged false
+shape
+valid
+reader urword
+optional false
+longname stress period number
+description REPLACE iper {}
+
+block period
+name reachperioddata
+type recarray rno abcsetting
+shape
+reader urword
+longname
+description
+
+block period
+name rno
+type integer
+shape
+tagged false
+in_record true
+reader urword
+longname reach number for this entry
+description integer value that defines the reach number associated with the specified PERIOD data on the line. RNO must be greater than zero and less than or equal to NREACHES.
+numeric_index true
+
+block period
+name abcsetting
+type keystring wspd tatm solr shd swrefl rh
+auxiliaryrecord
+shape
+tagged false
+in_record true
+reader urword
+longname
+description line of information that is parsed into a parameter keyword and values.  Property name keywords that can be used to start the ABCSETTING string include: WSPD, TATM, SOLR, SHD, SWREFL, and RH.  These settings are used to assign the wind speed, temperature of the atmosphere, shortwave solar radiation, channel shade, shortwave reflectance of the water surface, and relative humidity, respectively.  
+
+block period
+name wspd
+type string
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname wind speed (L/T)
+description real or character value that defines the wind speed above of each reach. If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
+
+block period
+name tatm
+type double precision
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname air temperature
+description is the temperature of the atmosphere above the stream reach. For an unrotated case this is the hydraulic conductivity in the y direction.  If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
+
+block period
+name solr
+type string
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname shortwave solar radiation (E/T/L^2)
+description real or character value that defines the solar radiation above of each reach. If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
+
+block period
+name shd
+type double precision
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname riparian shade
+description is the shaded fraction of the wetted surface of the stream reach. If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
+
+block period
+name swrefl
+type double precision
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname reflectance of shortwave radiation of water surface
+description the fraction of solar radiation reflected by water surface. If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
+
+block period
+name rh
+type double precision
+shape
+tagged true
+in_record true
+reader urword
+time_series true
+longname relative humidity of the atmosphere
+description the relative humidity of the atmosphere. If the Options block includes a TIMESERIESFILE entry (see the ``Time-Variable Input'' section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
diff --git a/doc/mf6io/mf6ivar/examples/gwe-sfe-example-obs.dat b/doc/mf6io/mf6ivar/examples/gwe-sfe-example-obs.dat
@@ -22,4 +22,5 @@ BEGIN continuous  FILEOUT  gwe_sfe02.sfe.obs.csv
   sfe-5-fjf          FLOW-JA-FACE   MYREACH1
   sfe-6-fjf          FLOW-JA-FACE   MYREACH2
   sfe-7-fjf          FLOW-JA-FACE   MYREACH3
+  sfe-1-shf          SHF            1
 END continuous
diff --git a/doc/mf6io/mf6ivar/examples/utl-shf-example.dat b/doc/mf6io/mf6ivar/examples/utl-shf-example.dat
@@ -0,0 +1,13 @@
+BEGIN OPTIONS
+  PRINT_INPUT
+  DENSITY_AIR  1.225
+  HEAT_CAPACITY_AIR  717.0
+  DRAG_COEFFICIENT  0.002
+END OPTIONS
+
+BEGIN PERIOD  1
+  1  WSPD    5.00
+  1  TATM   12.97
+  2  WSPD    4.99
+  2  TATM   12.95
+END PERIOD  1
diff --git a/doc/mf6io/mf6ivar/examples/utl-swr-example.dat b/doc/mf6io/mf6ivar/examples/utl-swr-example.dat
@@ -0,0 +1,12 @@
+BEGIN OPTIONS
+  PRINT_INPUT
+END OPTIONS
+
+BEGIN PERIOD  1
+  1  SOLR    800.0
+  1  SHD       0.1
+  1  SWREFL    0.03
+  2  SOLR    800.0
+  2  SHD       0.25
+  2  SWREFL    0.03
+END PERIOD  1