EXG-GWTGWT

Structure of Blocks

FOR EACH SIMULATION

    BEGIN OPTIONS
      GWFMODELNAME1 <gwfmodelname1>
      GWFMODELNAME2 <gwfmodelname2>
      [AUXILIARY <auxiliary(naux)>]
      [BOUNDNAMES]
      [PRINT_INPUT]
      [PRINT_FLOWS]
      [SAVE_FLOWS]
      [ADV_SCHEME <adv_scheme>]
      [DSP_XT3D_OFF]
      [DSP_XT3D_RHS]
      [MVT6 FILEIN <mvt6_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

  • gwfmodelname1 keyword to specify name of first corresponding GWF Model. In the simulation name file, the GWT6-GWT6 entry contains names for GWT Models (exgmnamea and exgmnameb). The GWT Model with the name exgmnamea must correspond to the GWF Model with the name gwfmodelname1.

  • gwfmodelname2 keyword to specify name of second corresponding GWF Model. In the simulation name file, the GWT6-GWT6 entry contains names for GWT Models (exgmnamea and exgmnameb). The GWT Model with the name exgmnameb must correspond to the GWF Model with the name gwfmodelname2.

  • 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 GWT-GWT 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” is provided in the description of the DISU Package. ANGLDEGX must be specified if dispersion is simulated in the connected GWT models.

  • BOUNDNAMES keyword to indicate that boundary names may be provided with the list of GWT 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.

  • adv_scheme scheme used to solve the advection term. Can be upstream, central, or TVD. If not specified, upstream weighting is the default weighting scheme.

  • DSP_XT3D_OFF deactivate the xt3d method for the dispersive flux and use the faster and less accurate approximation for this exchange.

  • DSP_XT3D_RHS add xt3d dispersion terms to right-hand side, when possible, for this exchange.

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

  • MVT6 keyword to specify that record corresponds to a transport mover file.

  • mvt6_filename is the file name of the transport mover input file to apply to this exchange. Information for the transport mover are provided in the file provided with these keywords.

  • 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:gwt-obstypetable} lists observation type(s) supported by the GWT-GWT package.

Block: DIMENSIONS

  • nexg keyword and integer value specifying the number of GWT-GWT 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 GWTGWT 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 GWT 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
    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