GWF-DISV

Structure of Blocks

FOR EACH SIMULATION

    BEGIN OPTIONS
      [LENGTH_UNITS <length_units>]
      [NOGRB]
      [XORIGIN <xorigin>]
      [YORIGIN <yorigin>]
      [ANGROT <angrot>]
      [EXPORT_ARRAY_ASCII]
    END OPTIONS
    BEGIN DIMENSIONS
      NLAY <nlay>
      NCPL <ncpl>
      NVERT <nvert>
    END DIMENSIONS
    BEGIN GRIDDATA
      TOP
            <top(ncpl)> -- READARRAY
      BOTM [LAYERED]
            <botm(ncpl, nlay)> -- READARRAY
      [IDOMAIN [LAYERED]
            <idomain(ncpl, nlay)> -- READARRAY]
    END GRIDDATA
    BEGIN VERTICES
      <iv> <xv> <yv>
      <iv> <xv> <yv>
      ...
    END VERTICES
    BEGIN CELL2D
      <icell2d> <xc> <yc> <ncvert> <icvert(ncvert)>
      <icell2d> <xc> <yc> <ncvert> <icvert(ncvert)>
      ...
    END CELL2D

Explanation of Variables

Block: OPTIONS

  • length_units is the length units used for this model. Values can be “FEET”, “METERS”, or “CENTIMETERS”. If not specified, the default is “UNKNOWN”.

  • NOGRB keyword to deactivate writing of the binary grid file.

  • xorigin x-position of the origin used for model grid vertices. This value should be provided in a real-world coordinate system. A default value of zero is assigned if not specified. The value for XORIGIN does not affect the model simulation, but it is written to the binary grid file so that postprocessors can locate the grid in space.

  • yorigin y-position of the origin used for model grid vertices. This value should be provided in a real-world coordinate system. If not specified, then a default value equal to zero is used. The value for YORIGIN does not affect the model simulation, but it is written to the binary grid file so that postprocessors can locate the grid in space.

  • angrot counter-clockwise rotation angle (in degrees) of the model grid coordinate system relative to a real-world coordinate system. If not specified, then a default value of 0.0 is assigned. The value for ANGROT does not affect the model simulation, but it is written to the binary grid file so that postprocessors can locate the grid in space.

  • EXPORT_ARRAY_ASCII keyword that specifies input griddata arrays should be written to layered ascii output files.

Block: DIMENSIONS

  • nlay is the number of layers in the model grid.

  • ncpl is the number of cells per layer. This is a constant value for the grid and it applies to all layers.

  • nvert is the total number of (x, y) vertex pairs used to characterize the horizontal configuration of the model grid.

Block: GRIDDATA

  • top is the top elevation for each cell in the top model layer.

  • botm is the bottom elevation for each cell.

  • idomain is an optional array that characterizes the existence status of a cell. If the IDOMAIN array is not specified, then all model cells exist within the solution. If the IDOMAIN value for a cell is 0, the cell does not exist in the simulation. Input and output values will be read and written for the cell, but internal to the program, the cell is excluded from the solution. If the IDOMAIN value for a cell is 1 or greater, the cell exists in the simulation. If the IDOMAIN value for a cell is -1, the cell does not exist in the simulation. Furthermore, the first existing cell above will be connected to the first existing cell below. This type of cell is referred to as a “vertical pass through” cell.

Block: VERTICES

  • iv is the vertex number. Records in the VERTICES block must be listed in consecutive order from 1 to NVERT.

  • xv is the x-coordinate for the vertex.

  • yv is the y-coordinate for the vertex.

Block: CELL2D

  • icell2d is the CELL2D number. Records in the CELL2D block must be listed in consecutive order from the first to the last.

  • xc is the x-coordinate for the cell center.

  • yc is the y-coordinate for the cell center.

  • ncvert is the number of vertices required to define the cell. There may be a different number of vertices for each cell.

  • icvert is an array of integer values containing vertex numbers (in the VERTICES block) used to define the cell. Vertices must be listed in clockwise order. Cells that are connected must share vertices.

Example Input File

    #The OPTIONS block is optional
    BEGIN OPTIONS
      LENGTH_UNITS METERS
    END OPTIONS
    
    #The DIMENSIONS block is required
    BEGIN DIMENSIONS
      NCPL 4
      NLAY 3
      NVERT 9
    END DIMENSIONS
    
    #The GRIDDATA block is required
    BEGIN GRIDDATA
      TOP
        CONSTANT 3.0
      BOTM LAYERED
        CONSTANT 2.0
        CONSTANT 1.0
        CONSTANT 0.0
      IDOMAIN LAYERED
        INTERNAL FACTOR 1
          1 1 1 0
        CONSTANT 1
        CONSTANT 1
    END GRIDDATA
    
    #The VERTICES block is required
    BEGIN VERTICES
      1 0. 1.
      2 .5 1.
      3 1. 1.
      4 0 .5
      5 .5 .5
      6 1. .5
      7 0. 0.
      8 .5 0.
      9 1. 0.
    END VERTICES
    
    BEGIN CELL2D
      1 .25 .75 4 1 2 5 4
      2 .75 .75 4 2 3 6 5
      3 .25 .25 4 4 5 8 7
      4 .75 .25 4 5 6 9 8
    END CELL2D