-
Notifications
You must be signed in to change notification settings - Fork 122
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
feat(swf-pcp): add precipitation package for CHF and OLF models (#2160)
* feat(swf-pcp): add precipitation package for CHF and OLF models * fprettify * update makefile * update io guide * apply precip to maximum top width of the channel ensure cross sections are defined by at least 2 points ensure precip is greater than or equal to zero
- Loading branch information
1 parent
2fb2872
commit 3f6af23
Showing
37 changed files
with
2,992 additions
and
156 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,47 @@ | ||
| Input to the Precipitation (PCP) Package is read from the file that has type ``PCP6'' in the Name File. Any number of PCP Packages can be specified for a single surface water flow model. | ||
|
|
||
| \vspace{5mm} | ||
| \subsubsection{Structure of Blocks} | ||
| \vspace{5mm} | ||
|
|
||
| \noindent \textit{FOR EACH SIMULATION} | ||
| \lstinputlisting[style=blockdefinition]{./mf6ivar/tex/chf-pcp-options.dat} | ||
| \lstinputlisting[style=blockdefinition]{./mf6ivar/tex/chf-pcp-dimensions.dat} | ||
| \vspace{5mm} | ||
| \noindent \textit{FOR ANY STRESS PERIOD} | ||
| \lstinputlisting[style=blockdefinition]{./mf6ivar/tex/chf-pcp-period.dat} | ||
| \packageperioddescription | ||
|
|
||
| \vspace{5mm} | ||
| \subsubsection{Explanation of Variables} | ||
| \begin{description} | ||
| \input{./mf6ivar/tex/chf-pcp-desc.tex} | ||
| \end{description} | ||
|
|
||
| \vspace{5mm} | ||
| \subsubsection{Example Input File} | ||
| \lstinputlisting[style=inputfile]{./mf6ivar/examples/chf-pcp-example.dat} | ||
|
|
||
| %\vspace{5mm} | ||
| %\subsubsection{Available observation types} | ||
| %Well Package observations include the simulated well rates (\texttt{wel}), the well discharge that is available for the MVR package (\texttt{to-mvr}), and the reduction in the specified \texttt{q} when the \texttt{AUTO\_FLOW\_REDUCE} option is enabled. The data required for each WEL Package observation type is defined in table~\ref{table:gwf-welobstype}. The sum of \texttt{wel} and \texttt{to-mvr} is equal to the simulated well discharge rate, which may be less than the specified \texttt{q} if the \texttt{AUTO\_FLOW\_REDUCE} option is enabled. The \texttt{DNODATA} value is returned if the \texttt{wel-reduction} observation is specified but the \texttt{AUTO\_FLOW\_REDUCE} option is not enabled. Negative and positive values for an observation represent a loss from and gain to the GWF model, respectively. | ||
|
|
||
| %\begin{longtable}{p{2cm} p{2.75cm} p{2cm} p{1.25cm} p{7cm}} | ||
| %\caption{Available WEL Package observation types} \tabularnewline | ||
|
|
||
| %\hline | ||
| %\hline | ||
| %\textbf{Stress Package} & \textbf{Observation type} & \textbf{ID} & \textbf{ID2} & \textbf{Description} \\ | ||
| %\hline | ||
| %\endhead | ||
|
|
||
| %\hline | ||
| %\endfoot | ||
|
|
||
| %\input{../Common/gwf-welobs.tex} | ||
| %\label{table:gwf-welobstype} | ||
| %\end{longtable} | ||
|
|
||
| %\vspace{5mm} | ||
| %\subsubsection{Example Observation Input File} | ||
| %\lstinputlisting[style=inputfile]{./mf6ivar/examples/gwf-wel-example-obs.dat} |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,211 @@ | ||
| # --------------------- swf pcp options --------------------- | ||
| # flopy multi-package | ||
|
|
||
| block options | ||
| name auxiliary | ||
| type string | ||
| shape (naux) | ||
| reader urword | ||
| optional true | ||
| longname keyword to specify aux variables | ||
| description REPLACE auxnames {'{#1}': 'Surface Water Flow'} | ||
|
|
||
| block options | ||
| name auxmultname | ||
| type string | ||
| shape | ||
| reader urword | ||
| optional true | ||
| longname name of auxiliary variable for multiplier | ||
| description REPLACE auxmultname {'{#1}': 'precipitation'} | ||
|
|
||
| block options | ||
| name boundnames | ||
| type keyword | ||
| shape | ||
| reader urword | ||
| optional true | ||
| longname | ||
| description REPLACE boundnames {'{#1}': 'precipitation'} | ||
|
|
||
| block options | ||
| name print_input | ||
| type keyword | ||
| reader urword | ||
| optional true | ||
| longname print input to listing file | ||
| description REPLACE print_input {'{#1}': 'precipitation'} | ||
| mf6internal iprpak | ||
|
|
||
| block options | ||
| name print_flows | ||
| type keyword | ||
| reader urword | ||
| optional true | ||
| longname print precipitation rates to listing file | ||
| description REPLACE print_flows {'{#1}': 'precipitation'} | ||
| mf6internal iprflow | ||
|
|
||
| block options | ||
| name save_flows | ||
| type keyword | ||
| reader urword | ||
| optional true | ||
| longname save precipitation to budget file | ||
| description REPLACE save_flows {'{#1}': 'precipitation'} | ||
| mf6internal ipakcb | ||
|
|
||
| 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 obs_filerecord | ||
| type record obs6 filein obs6_filename | ||
| shape | ||
| reader urword | ||
| tagged true | ||
| optional true | ||
| longname | ||
| description | ||
|
|
||
| block options | ||
| name obs6 | ||
| type keyword | ||
| shape | ||
| in_record true | ||
| reader urword | ||
| tagged true | ||
| optional false | ||
| longname obs keyword | ||
| description keyword to specify that record corresponds to an observations file. | ||
|
|
||
| block options | ||
| name obs6_filename | ||
| type string | ||
| preserve_case true | ||
| in_record true | ||
| tagged false | ||
| reader urword | ||
| optional false | ||
| longname obs6 input filename | ||
| description REPLACE obs6_filename {'{#1}': 'Precipitation'} | ||
|
|
||
| # --------------------- swf pcp dimensions --------------------- | ||
|
|
||
| block dimensions | ||
| name maxbound | ||
| type integer | ||
| reader urword | ||
| optional false | ||
| longname maximum number of precipitation cells | ||
| description REPLACE maxbound {'{#1}': 'precipitation cells'} | ||
|
|
||
|
|
||
| # --------------------- swf pcp 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 stress_period_data | ||
| type recarray cellid precipitation aux boundname | ||
| shape (maxbound) | ||
| reader urword | ||
| longname | ||
| description | ||
| mf6internal spd | ||
|
|
||
| block period | ||
| name cellid | ||
| type integer | ||
| shape (ncelldim) | ||
| tagged false | ||
| in_record true | ||
| reader urword | ||
| longname cell identifier | ||
| description REPLACE cellid {} | ||
|
|
||
| block period | ||
| name precipitation | ||
| type double precision | ||
| shape | ||
| tagged false | ||
| in_record true | ||
| reader urword | ||
| time_series true | ||
| longname precipitation rate | ||
| description is the precipitation flux rate ($LT^{-1}$). This rate is multiplied inside the program by the surface area of the cell to calculate the volumetric precipitation rate. 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 aux | ||
| type double precision | ||
| in_record true | ||
| tagged false | ||
| shape (naux) | ||
| reader urword | ||
| optional true | ||
| time_series true | ||
| longname auxiliary variables | ||
| description REPLACE aux {'{#1}': 'precipitation'} | ||
| mf6internal auxvar | ||
|
|
||
| block period | ||
| name boundname | ||
| type string | ||
| shape | ||
| tagged false | ||
| in_record true | ||
| reader urword | ||
| optional true | ||
| longname precipitation name | ||
| description REPLACE boundname {'{#1}': 'precipitation'} |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.