EXG-GWFGWF

Structure of Blocks

FOR EACH SIMULATION

    BEGIN OPTIONS
      [AUXILIARY <auxiliary(naux)>]
      [BOUNDNAMES]
      [PRINT_INPUT]
      [PRINT_FLOWS]
      [SAVE_FLOWS]
      [CELL_AVERAGING <cell_averaging>]
      [VARIABLECV [DEWATERED]]
      [NEWTON]
      [XT3D]
      [GNC6 FILEIN <gnc6_filename>]
      [MVR6 FILEIN <mvr6_filename>]
      [OBS6 FILEIN <obs6_filename>]
    END OPTIONS
    BEGIN DIMENSIONS
      NEXG <nexg>
    END DIMENSIONS
    BEGIN EXCHANGEDATA
      <cellidm1> <cellidm2> <ihc> <cl1> <cl2> <hwva> [<aux(naux)>] [<boundname>]
      <cellidm1> <cellidm2> <ihc> <cl1> <cl2> <hwva> [<aux(naux)>] [<boundname>]
      ...
    END EXCHANGEDATA

Explanation of Variables

Block: OPTIONS

  • auxiliary an array of auxiliary variable names. There is no limit on the number of auxiliary variables that can be provided. Most auxiliary variables will not be used by the GWF-GWF Exchange, but they will be available for use by other parts of the program. If an auxiliary variable with the name “ANGLDEGX” is found, then this information will be used as the angle (provided in degrees) between the connection face normal and the x axis, where a value of zero indicates that a normal vector points directly along the positive x axis. The connection face normal is a normal vector on the cell face shared between the cell in model 1 and the cell in model 2 pointing away from the model 1 cell. Additional information on “ANGLDEGX” and when it is required is provided in the description of the DISU Package. If an auxiliary variable with the name “CDIST” is found, then this information will be used in the calculation of specific discharge within model cells connected by the exchange. For a horizontal connection, CDIST should be specified as the horizontal distance between the cell centers, and should not include the vertical component. For vertical connections, CDIST should be specified as the difference in elevation between the two cell centers. Both ANGLDEGX and CDIST are required if specific discharge is calculated for either of the groundwater models.

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

  • PRINT_INPUT keyword to indicate that the list of exchange entries will be echoed to the listing file immediately after it is read.

  • PRINT_FLOWS keyword to indicate that the list of exchange flow rates will be printed to the listing file for every stress period in which “SAVE BUDGET” is specified in Output Control.

  • SAVE_FLOWS keyword to indicate that cell-by-cell flow terms will be written to the budget file for each model provided that the Output Control for the models are set up with the “BUDGET SAVE FILE” option.

  • cell_averaging is a keyword and text keyword to indicate the method that will be used for calculating the conductance for horizontal cell connections. The text value for CELL_AVERAGING can be “HARMONIC”, “LOGARITHMIC”, or “AMT-LMK”, which means “arithmetic-mean thickness and logarithmic-mean hydraulic conductivity”. If the user does not specify a value for CELL_AVERAGING, then the harmonic-mean method will be used.

  • VARIABLECV keyword to indicate that the vertical conductance will be calculated using the saturated thickness and properties of the overlying cell and the thickness and properties of the underlying cell. If the DEWATERED keyword is also specified, then the vertical conductance is calculated using only the saturated thickness and properties of the overlying cell if the head in the underlying cell is below its top. If these keywords are not specified, then the default condition is to calculate the vertical conductance at the start of the simulation using the initial head and the cell properties. The vertical conductance remains constant for the entire simulation.

  • DEWATERED If the DEWATERED keyword is specified, then the vertical conductance is calculated using only the saturated thickness and properties of the overlying cell if the head in the underlying cell is below its top.

  • NEWTON keyword that activates the Newton-Raphson formulation for groundwater flow between connected, convertible groundwater cells. Cells will not dry when this option is used.

  • XT3D keyword that activates the XT3D formulation between the cells connected with this GWF-GWF Exchange.

  • FILEIN keyword to specify that an input filename is expected next.

  • GNC6 keyword to specify that record corresponds to a ghost-node correction file.

  • gnc6_filename is the file name for ghost node correction input file. Information for the ghost nodes are provided in the file provided with these keywords. The format for specifying the ghost nodes is the same as described for the GNC Package of the GWF Model. This includes specifying OPTIONS, DIMENSIONS, and GNCDATA blocks. The order of the ghost nodes must follow the same order as the order of the cells in the EXCHANGEDATA block. For the GNCDATA, noden and all of the nodej values are assumed to be located in model 1, and nodem is assumed to be in model 2.

  • MVR6 keyword to specify that record corresponds to a mover file.

  • mvr6_filename is the file name of the water mover input file to apply to this exchange. Information for the water mover are provided in the file provided with these keywords. The format for specifying the water mover information is the same as described for the Water Mover (MVR) Package of the GWF Model, with two exceptions. First, in the PACKAGES block, the model name must be included as a separate string before each package. Second, the appropriate model name must be included before package name 1 and package name 2 in the BEGIN PERIOD block. This allows providers and receivers to be located in both models listed as part of this exchange.

  • OBS6 keyword to specify that record corresponds to an observations file.

  • obs6_filename is the file name of the observations input file for this exchange. See the “Observation utility” section for instructions for preparing observation input files. Table ref{table:gwf-obstypetable} lists observation type(s) supported by the GWF-GWF package.

Block: DIMENSIONS

  • nexg keyword and integer value specifying the number of GWF-GWF exchanges.

Block: EXCHANGEDATA

  • cellidm1 is the cellid of the cell in model 1 as specified in the simulation name file. For a structured grid that uses the DIS input file, CELLIDM1 is the layer, row, and column numbers of the cell. For a grid that uses the DISV input file, CELLIDM1 is the layer number and CELL2D number for the two cells. If the model uses the unstructured discretization (DISU) input file, then CELLIDM1 is the node number for the cell.

  • cellidm2 is the cellid of the cell in model 2 as specified in the simulation name file. For a structured grid that uses the DIS input file, CELLIDM2 is the layer, row, and column numbers of the cell. For a grid that uses the DISV input file, CELLIDM2 is the layer number and CELL2D number for the two cells. If the model uses the unstructured discretization (DISU) input file, then CELLIDM2 is the node number for the cell.

  • ihc is an integer flag indicating the direction between node n and all of its m connections. If IHC = 0 then the connection is vertical. If IHC = 1 then the connection is horizontal. If IHC = 2 then the connection is horizontal for a vertically staggered grid.

  • cl1 is the distance between the center of cell 1 and the its shared face with cell 2.

  • cl2 is the distance between the center of cell 2 and the its shared face with cell 1.

  • hwva is the horizontal width of the flow connection between cell 1 and cell 2 if IHC > 0, or it is the area perpendicular to flow of the vertical connection between cell 1 and cell 2 if IHC = 0.

  • aux represents the values of the auxiliary variables for each GWFGWF Exchange. The values of auxiliary variables must be present for each exchange. The values must be specified in the order of the auxiliary variables specified in the OPTIONS block.

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

Example Input File

    BEGIN OPTIONS
      PRINT_INPUT
      PRINT_FLOWS
      SAVE_FLOWS
      AUXILIARY testaux
      GNC6 FILEIN simulation.gnc
      MVR6 FILEIN simulation.mvr
    END OPTIONS
    
    BEGIN DIMENSIONS
      NEXG 36
    END DIMENSIONS
    
      # nodem1 nodem2  ihc  cl1   cl2    fahl  testaux
    BEGIN EXCHANGEDATA
    #
    #       left side
            16      1    1  50. 16.67 33.33   100.99
            16     10    1  50. 16.67 33.33   100.99
            16     19    1  50. 16.67 33.33   100.99
            23     28    1  50. 16.67 33.33   100.99
            23     37    1  50. 16.67 33.33   100.99
            23     46    1  50. 16.67 33.33   100.99
            30     55    1  50. 16.67 33.33   100.99
            30     64    1  50. 16.67 33.33   100.99
            30     73    1  50. 16.67 33.33   100.99
    #
    #       right side
            20      9    1  50. 16.67 33.33   100.99
            20     18    1  50. 16.67 33.33   100.99
            20     27    1  50. 16.67 33.33   100.99
            27     36    1  50. 16.67 33.33   100.99
            27     45    1  50. 16.67 33.33   100.99
            27     54    1  50. 16.67 33.33   100.99
            34     63    1  50. 16.67 33.33   100.99
            34     72    1  50. 16.67 33.33   100.99
            34     81    1  50. 16.67 33.33   100.99
    #
    #       back
            10      1    1  50. 17.67 33.33   100.99
            10      2    1  50. 17.67 33.33   100.99
            10      3    1  50. 17.67 33.33   100.99
            11      4    1  50. 17.67 33.33   100.99
            11      5    1  50. 17.67 33.33   100.99
            11      6    1  50. 17.67 33.33   100.99
            12      7    1  50. 17.67 33.33   100.99
            12      8    1  50. 17.67 33.33   100.99
            12      9    1  50. 17.67 33.33   100.99
    #
    #       front
            38     73    1  50. 17.67 33.33   100.99
            38     74    1  50. 17.67 33.33   100.99
            38     75    1  50. 17.67 33.33   100.99
            39     76    1  50. 17.67 33.33   100.99
            39     77    1  50. 17.67 33.33   100.99
            39     78    1  50. 17.67 33.33   100.99
            40     79    1  50. 17.67 33.33   100.99
            40     80    1  50. 17.67 33.33   100.99
            40     81    1  50. 17.67 33.33   100.99
    END EXCHANGEDATA

Example Observation Input File

    BEGIN OPTIONS
      DIGITS 10
      PRINT_INPUT
    END OPTIONS
    
    # Block defining continuous observations
    BEGIN CONTINUOUS FILEOUT simulation.obs.csv
    # obsname        obstype        id or boundname
      exg1           flow-ja-face   1
      left-face      flow-ja-face   bnameleft
      right-face     flow-ja-face   bnameright
    END CONTINUOUS