GWE-LKE

Structure of Blocks

FOR EACH SIMULATION

    BEGIN OPTIONS
      [FLOW_PACKAGE_NAME <flow_package_name>]
      [AUXILIARY <auxiliary>]
      [FLOW_PACKAGE_AUXILIARY_NAME <flow_package_auxiliary_name>]
      [BOUNDNAMES]
      [PRINT_INPUT]
      [PRINT_TEMPERATURE]
      [PRINT_FLOWS]
      [SAVE_FLOWS]
      [TEMPERATURE FILEOUT <tempfile>]
      [BUDGET FILEOUT <budgetfile>]
      [BUDGETCSV FILEOUT <budgetcsvfile>]
      [TS6 FILEIN <ts6_filename>]
      [OBS6 FILEIN <obs6_filename>]
    END OPTIONS
    BEGIN PACKAGEDATA
      <lakeno> <strt> <ktf> <rbthcnd> [<aux(naux)>] [<boundname>]
      <lakeno> <strt> <ktf> <rbthcnd> [<aux(naux)>] [<boundname>]
      ...
    END PACKAGEDATA

FOR ANY STRESS PERIOD

    BEGIN PERIOD <iper>
      <lakeno> <laksetting>
      <lakeno> <laksetting>
      ...
    END PERIOD

Explanation of Variables

Block: OPTIONS

  • flow_package_name keyword to specify the name of the corresponding flow package. If not specified, then the corresponding flow package must have the same name as this advanced transport package (the name associated with this package in the GWE name file).

  • 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.

  • flow_package_auxiliary_name keyword to specify the name of an auxiliary variable in the corresponding flow package. If specified, then the simulated temperatures from this advanced transport package will be copied into the auxiliary variable specified with this name. Note that the flow package must have an auxiliary variable with this name or the program will terminate with an error. If the flows for this advanced transport package are read from a file, then this option will have no effect.

  • BOUNDNAMES keyword to indicate that boundary names may be provided with the list of lake cells.

  • PRINT_INPUT keyword to indicate that the list of lake information will be written to the listing file immediately after it is read.

  • PRINT_TEMPERATURE keyword to indicate that the list of lake temperature will be printed to the listing file for every stress period in which “TEMPERATURE PRINT” is specified in Output Control. If there is no Output Control option and PRINT_TEMPERATURE is specified, then temperature are printed for the last time step of each stress period.

  • PRINT_FLOWS keyword to indicate that the list of lake 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 lake flow terms will be written to the file specified with “BUDGET FILEOUT” in Output Control.

  • TEMPERATURE keyword to specify that record corresponds to temperature.

  • tempfile name of the binary output file to write temperature 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.

  • 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 LKE package. See the “Observation utility” section for instructions for preparing observation input files. Tables ref{table:gwf-obstypetable} and ref{table:gwt-obstypetable} lists observation type(s) supported by the LKE package.

Block: PACKAGEDATA

  • lakeno integer value that defines the lake number associated with the specified PACKAGEDATA data on the line. LAKENO must be greater than zero and less than or equal to NLAKES. Lake information must be specified for every lake or the program will terminate with an error. The program will also terminate with an error if information for a lake is specified more than once.

  • strt real value that defines the starting temperature for the lake.

  • ktf is the thermal conductivity of the material between the aquifer cell and the lake. The thickness of the material is defined by the variable RBTHCND.

  • rbthcnd real value that defines the thickness of the lakebed material through which conduction occurs. Must be greater than 0.

  • aux represents the values of the auxiliary variables for each lake. The values of auxiliary variables must be present for each lake. 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.

  • boundname name of the lake 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.

  • lakeno integer value that defines the lake number associated with the specified PERIOD data on the line. LAKENO must be greater than zero and less than or equal to NLAKES.

  • laksetting line of information that is parsed into a keyword and values. Keyword values that can be used to start the LAKSETTING string include: STATUS, TEMPERATURE, RAINFALL, EVAPORATION, RUNOFF, and AUXILIARY. These settings are used to assign the temperature associated with the corresponding flow terms. Temperatures cannot be specified for all flow terms. For example, the Lake Package supports a “WITHDRAWAL” flow term. If this withdrawal term is active, then water will be withdrawn from the lake at the calculated temperature of the lake.

      STATUS <status>
      TEMPERATURE <temperature>
      RAINFALL <rainfall>
      EVAPORATION <evaporation>
      RUNOFF <runoff>
      EXT-INFLOW <ext-inflow>
      AUXILIARY <auxname> <auxval> 
    
  • status keyword option to define lake status. STATUS can be ACTIVE, INACTIVE, or CONSTANT. By default, STATUS is ACTIVE, which means that temperature will be calculated for the lake. If a lake is inactive, then there will be no energy fluxes into or out of the lake and the inactive value will be written for the lake temperature. If a lake is constant, then the temperature for the lake will be fixed at the user specified value.

  • temperature real or character value that defines the temperature for the lake. The specified TEMPERATURE is only applied if the lake is a constant temperature lake. 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.

  • rainfall real or character value that defines the rainfall temperature for the lake. 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.

  • evaporation use of the EVAPORATION keyword is allowed in the LKE package; however, the specified value is not currently used in LKE calculations. Instead, the latent heat of evaporation is multiplied by the simulated evaporation rate for determining the thermal energy lost from a stream reach.

  • runoff real or character value that defines the temperature of runoff for the lake. Users are free to use whatever temperature scale they want, which might include negative temperatures. 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.

  • ext-inflow real or character value that defines the temperature of external inflow for the lake. Users are free to use whatever temperature scale they want, which might include negative temperatures. 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.

  • AUXILIARY keyword for specifying auxiliary variable.

  • auxname name for the auxiliary variable to be assigned AUXVAL. AUXNAME must match one of the auxiliary variable names defined in the OPTIONS block. If AUXNAME does not match one of the auxiliary variable names defined in the OPTIONS block the data are ignored.

  • auxval value for the auxiliary variable. 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.

Example Input File

    BEGIN OPTIONS
      AUXILIARY  aux1  aux2
      BOUNDNAMES
      PRINT_INPUT
      PRINT_TEMPERATURE
      PRINT_FLOWS
      SAVE_FLOWS
      TEMPERATURE  FILEOUT  gwe_lke_02.lke.bin
      BUDGET  FILEOUT  gwe_lke_02.lke.bud
      OBS6  FILEIN  gwe_lke_02.lke.obs
    END OPTIONS
    
    BEGIN PACKAGEDATA
    # L      STRT      aux1      aux2       bname
      1       5.0      99.0     999.0     MYLAKE1
      2       6.0      99.0     999.0     MYLAKE2
      3       7.0      99.0     999.0     MYLAKE3
    END PACKAGEDATA
    
    BEGIN PERIOD  1
      1  STATUS  ACTIVE
      2  STATUS  ACTIVE
      3  STATUS  ACTIVE
    END PERIOD  1

Available Observation Types

Stress Package Observation Type ID1 ID2 Description
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 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 iconn for lake lakeno.
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.

Example Observation Input File

    BEGIN options
      DIGITS  7
      PRINT_INPUT
    END options
    
    BEGIN continuous  FILEOUT  gwe_lke02.lke.obs.csv
    # obsname            obstype        id1        id2
      lke-1-temp         TEMPERATURE    1
      lke-1-extinflow    EXT-INFLOW     1
      lke-1-rain         RAINFALL       1
      lke-1-roff         RUNOFF         1
      lke-1-wdrl         WITHDRAWAL     1
      lke-1-stor         STORAGE        1
      lke-1-const        CONSTANT       1
      lke-1-gwe1         LKE            1          1
      lke-1-gwe2         LKE            1          2
      lke-2-gwe1         LKE            2          1
      lke-1-mylake1      LKE            MYLAKE1
      lke-1-fjf          FLOW-JA-FACE   1          2
      lke-2-fjf          FLOW-JA-FACE   2          1
      lke-3-fjf          FLOW-JA-FACE   2          3
      lke-4-fjf          FLOW-JA-FACE   3          2
      lke-5-fjf          FLOW-JA-FACE   MYLAKE1
      lke-6-fjf          FLOW-JA-FACE   MYLAKE2
      lke-7-fjf          FLOW-JA-FACE   MYLAKE3
    END continuous