GWF-RCHA

Structure of Blocks

FOR EACH SIMULATION

    BEGIN OPTIONS
      READASARRAYS
      [FIXED_CELL]
      [AUXILIARY <auxiliary(naux)>]
      [AUXMULTNAME <auxmultname>]
      [PRINT_INPUT]
      [PRINT_FLOWS]
      [SAVE_FLOWS]
      [TAS6 FILEIN <tas6_filename>]
      [OBS6 FILEIN <obs6_filename>]
      [EXPORT_ARRAY_NETCDF]
    END OPTIONS

FOR ANY STRESS PERIOD

    BEGIN PERIOD <iper>
      [IRCH
            <irch(ncol*nrow; ncpl)> -- READARRAY]
      RECHARGE
            <recharge(ncol*nrow; ncpl)> -- READARRAY
      [AUX
            <aux(ncol*nrow; ncpl)> -- READARRAY]
    END PERIOD

Explanation of Variables

Block: OPTIONS

  • READASARRAYS indicates that array-based input will be used for the Recharge Package. This keyword must be specified to use array-based input. When READASARRAYS is specified, values must be provided for every cell within a model layer, even those cells that have an IDOMAIN value less than one. Values assigned to cells with IDOMAIN values less than one are not used and have no effect on simulation results.

  • FIXED_CELL indicates that recharge will not be reassigned to a cell underlying the cell specified in the list if the specified cell is inactive.

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

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

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

  • TAS6 keyword to specify that record corresponds to a time-array-series file.

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

  • tas6_filename defines a time-array-series file defining a time-array series that can be used to assign time-varying values. See the “Time-Variable Input” section for instructions on using the time-array series capability.

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

  • obs6_filename name of input file to define observations for the Recharge 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 Recharge package.

  • EXPORT_ARRAY_NETCDF keyword that specifies input griddata arrays should be written to the model output netcdf file.

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.

  • irch IRCH is the layer number that defines the layer in each vertical column where recharge is applied. If IRCH is omitted, recharge by default is applied to cells in layer 1. IRCH can only be used if READASARRAYS is specified in the OPTIONS block. If IRCH is specified, it must be specified as the first variable in the PERIOD block or MODFLOW will terminate with an error.

  • recharge is the recharge flux rate (LT-1). This rate is multiplied inside the program by the surface area of the cell to calculate the volumetric recharge rate. The recharge array may be defined by a time-array series (see the “Using Time-Array Series in a Package” section).

  • aux is an array of values for auxiliary variable aux(iaux), where iaux is a value from 1 to naux, and aux(iaux) must be listed as part of the auxiliary variables. A separate array can be specified for each auxiliary variable. If an array is not specified for an auxiliary variable, then a value of zero is assigned. If the value specified here for the auxiliary variable is the same as auxmultname, then the recharge array will be multiplied by this array.

Example Input File

    BEGIN OPTIONS
      AUXILIARY var1 var2 mymult
      READASARRAYS
      PRINT_INPUT
      PRINT_FLOWS
      SAVE_FLOWS
      AUXMULTNAME mymult
    END OPTIONS
    
    BEGIN PERIOD 1
    
      # For this model, the absence of an IRCH array causes
      # recharge to apply to model layer 1. To make recharge
      # apply to layer 2 instead, the following lines
      # (uncommented) could be used:
      # IRCH
      #   constant 2
    
      # recharge rate
      RECHARGE
        constant 0.0040
    
      # auxiliary variable (var1) array
      var1
        constant 100.
    
      # auxiliary variable (var2) array
      var2
        constant 0.
    
      # auxiliary variable (mymult) array
      # Because ``AUXMULTNAME mymult'' was specified in the
      # options block, the MYMULT array will be used to multiply
      # the values in the RECHARGE array
      MYMULT
        INTERNAL FACTOR  1.0
          0.5  1.0  1.0  0.5
          1.0  1.0  1.0  1.0
          0.5  1.0  1.0  0.5
    
    END PERIOD