### GWF-UZF #### Structure of Blocks _FOR EACH SIMULATION_ ``` BEGIN OPTIONS [AUXILIARY ] [AUXMULTNAME ] [BOUNDNAMES] [PRINT_INPUT] [PRINT_FLOWS] [SAVE_FLOWS] [WATER_CONTENT FILEOUT ] [BUDGET FILEOUT ] [BUDGETCSV FILEOUT ] [PACKAGE_CONVERGENCE FILEOUT ] [TS6 FILEIN ] [OBS6 FILEIN ] [MOVER] [SIMULATE_ET] [LINEAR_GWET] [SQUARE_GWET] [UNSAT_ETWC] [UNSAT_ETAE] END OPTIONS ``` ``` BEGIN DIMENSIONS NUZFCELLS NTRAILWAVES NWAVESETS END DIMENSIONS ``` ``` BEGIN PACKAGEDATA [] [] ... END PACKAGEDATA ``` _FOR ANY STRESS PERIOD_ ``` BEGIN PERIOD [] [] ... END PERIOD ``` #### Explanation of Variables ##### Block: OPTIONS * `auxiliary` defines an array of one or more auxiliary variable names. There is no limit on the number of auxiliary variables that can be provided on this line; however, lists of information provided in subsequent blocks must have a column of data for each auxiliary variable name defined here. The number of auxiliary variables detected on this line determines the value for naux. Comments cannot be provided anywhere on this line as they will be interpreted as auxiliary variable names. Auxiliary variables may not be used by the package, but they will be available for use by other parts of the program. The program will terminate with an error if auxiliary variables are specified on more than one line in the options block. * `auxmultname` name of auxiliary variable to be used as multiplier of GWF cell area used by UZF cell. * `BOUNDNAMES` keyword to indicate that boundary names may be provided with the list of UZF cells. * `PRINT_INPUT` keyword to indicate that the list of UZF information will be written to the listing file immediately after it is read. * `PRINT_FLOWS` keyword to indicate that the list of UZF flow rates will be printed to the listing file for every stress period time step in which "BUDGET PRINT" is specified in Output Control. If there is no Output Control option and "PRINT_FLOWS" is specified, then flow rates are printed for the last time step of each stress period. * `SAVE_FLOWS` keyword to indicate that UZF flow terms will be written to the file specified with "BUDGET FILEOUT" in Output Control. * `WATER_CONTENT` keyword to specify that record corresponds to unsaturated zone water contents. * `wcfile` name of the binary output file to write water content information. * `BUDGET` keyword to specify that record corresponds to the budget. * `FILEOUT` keyword to specify that an output filename is expected next. * `budgetfile` name of the binary output file to write budget information. * `BUDGETCSV` keyword to specify that record corresponds to the budget CSV. * `budgetcsvfile` name of the comma-separated value (CSV) output file to write budget summary information. A budget summary record will be written to this file for each time step of the simulation. * `PACKAGE_CONVERGENCE` keyword to specify that record corresponds to the package convergence comma spaced values file. * `package_convergence_filename` name of the comma spaced values output file to write package convergence information. * `TS6` keyword to specify that record corresponds to a time-series file. * `FILEIN` keyword to specify that an input filename is expected next. * `ts6_filename` defines a time-series file defining time series that can be used to assign time-varying values. See the "Time-Variable Input" section for instructions on using the time-series capability. * `OBS6` keyword to specify that record corresponds to an observations file. * `obs6_filename` name of input file to define observations for the UZF package. See the "Observation utility" section for instructions for preparing observation input files. Table ref{table:gwf-obstypetable} lists observation type(s) supported by the UZF package. * `MOVER` keyword to indicate that this instance of the UZF Package can be used with the Water Mover (MVR) Package. When the MOVER option is specified, additional memory is allocated within the package to store the available, provided, and received water. * `SIMULATE_ET` keyword specifying that ET in the unsaturated (UZF) and saturated zones (GWF) will be simulated. ET can be simulated in the UZF cell and not the GWF cell by omitting keywords LINEAR_GWET and SQUARE_GWET. * `LINEAR_GWET` keyword specifying that groundwater ET will be simulated using the original ET formulation of MODFLOW-2005. * `SQUARE_GWET` keyword specifying that groundwater ET will be simulated by assuming a constant ET rate for groundwater levels between land surface (TOP) and land surface minus the ET extinction depth (TOP-EXTDP). Groundwater ET is smoothly reduced from the PET rate to zero over a nominal interval at TOP-EXTDP. * `UNSAT_ETWC` keyword specifying that ET in the unsaturated zone will be simulated as a function of the specified PET rate while the water content (THETA) is greater than the ET extinction water content (EXTWC). * `UNSAT_ETAE` keyword specifying that ET in the unsaturated zone will be simulated using a capillary pressure based formulation. Capillary pressure is calculated using the Brooks-Corey retention function. ##### Block: DIMENSIONS * `nuzfcells` is the number of UZF cells. More than one UZF cell can be assigned to a GWF cell; however, only one GWF cell can be assigned to a single UZF cell. If more than one UZF cell is assigned to a GWF cell, then an auxiliary variable should be used to reduce the surface area of the UZF cell with the AUXMULTNAME option. * `ntrailwaves` is the number of trailing waves. A recommended value of 7 can be used for NTRAILWAVES. This value can be increased to lower mass balance error in the unsaturated zone. * `nwavesets` is the number of wave sets. A recommended value of 40 can be used for NWAVESETS. This value can be increased if more waves are required to resolve variations in water content within the unsaturated zone. ##### Block: PACKAGEDATA * `ifno` integer value that defines the feature (UZF object) number associated with the specified PACKAGEDATA data on the line. IFNO must be greater than zero and less than or equal to NUZFCELLS. UZF information must be specified for every UZF cell or the program will terminate with an error. The program will also terminate with an error if information for a UZF cell is specified more than once. * `cellid` is the cell identifier, and depends on the type of grid that is used for the simulation. For a structured grid that uses the DIS input file, CELLID is the layer, row, and column. For a grid that uses the DISV input file, CELLID is the layer and CELL2D number. If the model uses the unstructured discretization (DISU) input file, CELLID is the node number for the cell. * `landflag` integer value set to one for land surface cells indicating that boundary conditions can be applied and data can be specified in the PERIOD block. A value of 0 specifies a non-land surface cell. * `ivertcon` integer value set to specify underlying UZF cell that receives water flowing to bottom of cell. If unsaturated zone flow reaches the water table before the cell bottom, then water is added to the GWF cell instead of flowing to the underlying UZF cell. A value of 0 indicates the UZF cell is not connected to an underlying UZF cell. * `surfdep` is the surface depression depth of the UZF cell. * `vks` is the saturated vertical hydraulic conductivity of the UZF cell. This value is used with the Brooks-Corey function and the simulated water content to calculate the partially saturated hydraulic conductivity. * `thtr` is the residual (irreducible) water content of the UZF cell. This residual water is not available to plants and will not drain into underlying aquifer cells. * `thts` is the saturated water content of the UZF cell. The values for saturated and residual water content should be set in a manner that is consistent with the specific yield value specified in the Storage Package. The saturated water content must be greater than the residual content. * `thti` is the initial water content of the UZF cell. The value must be greater than or equal to the residual water content and less than or equal to the saturated water content. * `eps` is the exponent used in the Brooks-Corey function. The Brooks-Corey function is used by UZF to calculated hydraulic conductivity under partially saturated conditions as a function of water content and the user-specified saturated hydraulic conductivity. * `boundname` name of the UZF cell cell. BOUNDNAME is an ASCII character variable that can contain as many as 40 characters. If BOUNDNAME contains spaces in it, then the entire name must be enclosed within single quotes. ##### Block: PERIOD * `iper` integer value specifying the starting stress period number for which the data specified in the PERIOD block apply. IPER must be less than or equal to NPER in the TDIS Package and greater than zero. The IPER value assigned to a stress period block must be greater than the IPER value assigned for the previous PERIOD block. The information specified in the PERIOD block will continue to apply for all subsequent stress periods, unless the program encounters another PERIOD block. * `ifno` integer value that defines the feature (UZF object) number associated with the specified PERIOD data on the line. * `finf` real or character value that defines the applied infiltration rate of the UZF cell (LT-1). 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. * `pet` real or character value that defines the potential evapotranspiration rate of the UZF cell and specified GWF cell. Evapotranspiration is first removed from the unsaturated zone and any remaining potential evapotranspiration is applied to the saturated zone. If IVERTCON is greater than zero then residual potential evapotranspiration not satisfied in the UZF cell is applied to the underlying UZF and GWF cells. PET is always specified, but is only used if SIMULATE_ET is specified in the OPTIONS block. 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. * `extdp` real or character value that defines the evapotranspiration extinction depth of the UZF cell. If IVERTCON is greater than zero and EXTDP extends below the GWF cell bottom then remaining potential evapotranspiration is applied to the underlying UZF and GWF cells. EXTDP is always specified, but is only used if SIMULATE_ET is specified in the OPTIONS block. 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. * `extwc` real or character value that defines the evapotranspiration extinction water content of the UZF cell. EXTWC is always specified, but is only used if SIMULATE_ET and UNSAT_ETWC are specified in the OPTIONS block. The evapotranspiration rate from the unsaturated zone will be set to zero when the calculated water content is at or less than this value. The value for EXTWC cannot be less than the residual water content, and if it is specified as being less than the residual water content it is set to the residual water content. 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. * `ha` real or character value that defines the air entry potential (head) of the UZF cell. HA is always specified, but is only used if SIMULATE_ET and UNSAT_ETAE are specified in the OPTIONS block. 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. * `hroot` real or character value that defines the root potential (head) of the UZF cell. HROOT is always specified, but is only used if SIMULATE_ET and UNSAT_ETAE are specified in the OPTIONS block. 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. * `rootact` real or character value that defines the root activity function of the UZF cell. ROOTACT is the length of roots in a given volume of soil divided by that volume. Values range from 0 to about 3 cm-2, depending on the plant community and its stage of development. ROOTACT is always specified, but is only used if SIMULATE_ET and UNSAT_ETAE are specified in the OPTIONS block. 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. * `aux` represents the values of the auxiliary variables for each UZF. The values of auxiliary variables must be present for each UZF. The values must be specified in the order of the auxiliary variables specified in the OPTIONS block. If the package supports time series and 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. #### Example Input File ``` BEGIN OPTIONS OBS6 UZF.obs SIMULATE_ET UNSAT_ETWC LINEAR_GWET END OPTIONS BEGIN DIMENSIONS NUZFCELLS 10 NTRAILWAVES 7 NWAVESETS 40 END DIMENSIONS BEGIN PACKAGEDATA 1 1 1 1 1 1.0 1.0 0.05 0.35 0.1 4.0 2 1 1 2 1 1.0 1.0 0.05 0.35 0.1 4.0 3 1 1 3 1 1.0 1.0 0.05 0.35 0.1 4.0 4 1 1 4 1 1.0 1.0 0.05 0.35 0.1 4.0 5 1 1 5 1 1.0 1.0 0.05 0.35 0.1 4.0 6 1 1 6 1 1.0 1.0 0.05 0.35 0.1 4.0 7 1 1 7 1 1.0 1.0 0.05 0.35 0.1 4.0 8 1 1 8 1 1.0 1.0 0.05 0.35 0.1 4.0 9 1 1 9 1 1.0 1.0 0.05 0.35 0.1 4.0 10 1 1 10 1 1.0 1.0 0.05 0.35 0.1 4.0 END PACKAGEDATA BEGIN PERIOD 1 2 0.00005 0.00002 2.0 0.10 3 0.00008 0.00002 2.0 0.10 4 0.00009 0.00002 2.0 0.10 5 0.0001 0.00002 2.0 0.10 6 0.0001 0.00002 2.0 0.10 7 0.00009 0.00002 2.0 0.10 8 0.00008 0.00002 2.0 0.10 9 0.00005 0.00002 2.0 0.10 END PERIOD BEGIN PERIOD 2 2 0.00009 0.00003 2.0 0.10 3 0.0001 0.00003 2.0 0.10 4 0.0001 0.00003 2.0 0.10 5 0.00015 0.00003 2.0 0.10 6 0.00015 0.00003 2.0 0.10 7 0.0001 0.00003 2.0 0.10 8 0.0001 0.00003 2.0 0.10 9 0.00009 0.00003 2.0 0.10 END PERIOD ``` #### Available Observation Types | Stress Package | Observation Type | ID1 | ID2 | Description | |----------------|------------------|-----|-----|-------------| | UZF | uzf-gwrch | ifno or boundname | -- | Simulated recharge to the aquifer calculated by the UZF package for a UZF cell or a group of UZF cells. | | UZF | uzf-gwd | ifno or boundname | -- | Simulated groundwater discharge to the land surface calculated by the UZF package for a UZF cell or a group of UZF cells. | | UZF | uzf-gwd-to-mvr | ifno or boundname | -- | Simulated groundwater discharge to the land surface calculated by the UZF package that is available to the MVR package for a UZF cell or a group of UZF cells. | | UZF | uzf-gwet | ifno or boundname | -- | Simulated groundwater evapotranspiration calculated by the UZF package for a UZF cell or a group of UZF cells. | | UZF | infiltration | ifno or boundname | -- | Specified infiltration rate applied to a UZF package for a UZF cell or a group of UZF cells with landflag values not equal to zero. | | UZF | from-mvr | ifno or boundname | -- | Inflow into a UZF cell from the MVR package for a UZF cell or a group of UZF cells. | | UZF | rej-inf | ifno or boundname | -- | Simulated rejected infiltration calculated by the UZF package for a UZF cell or a group of UZF cells. | | UZF | rej-inf-to-mvr | ifno or boundname | -- | Simulated rejected infiltration calculated by the UZF package that is available to the MVR package for a UZF cell or a group of UZF cells. | | UZF | uzet | ifno or boundname | -- | Simulated unsaturated evapotranspiration calculated by the UZF package for a UZF cell or a group of UZF cells. | | UZF | storage | ifno or boundname | -- | Simulated storage flow rate for a UZF package cell or a group of UZF cells. | | UZF | net-infiltration | ifno or boundname | -- | Simulated net infiltration rate for a UZF package cell or a group of UZF cells. | | UZF | water-content | ifno or boundname | depth | Unsaturated-zone water content at a user-specified depth (ID2) relative to the top of GWF cellid for a UZF cell. The user-specified depth must be greater than or equal to zero and less than the thickness of GWF cellid (TOP - BOT). If boundname is specified, boundname must be unique for each UZF cell. | #### Example Observation Input File ``` BEGIN CONTINUOUS FILEOUT my_model.obs.uzf.csv id26_infil infiltration 26 id126_infil infiltration 126 id26_dpth=20 water-content 26 20.0 id126_dpth=51 water-content 126 1.0 #depth is below celtop id126_rch uzf-gwrch 126 END CONTINUOUS BEGIN CONTINUOUS FILEOUT my_model.uzf.budget.uzf.csv sinf infiltration uzfcells frommvr from-mvr uzfcells rejinf rej-inf uzfcells rejinftomvr rej-inf-to-mvr uzfcells uzet uzet uzfcells storage storage uzfcells net-inf net-infiltration uzfcells END CONTINUOUS BEGIN CONTINUOUS FILEOUT my_model.uzf.budget.gwf.csv gwrch uzf-gwrch uzfcells gwd uzf-gwd uzfcells gwdtomvr uzf-gwd-to-mvr uzfcells gwet uzf-gwet uzfcells END CONTINUOUS ```