flopy.discretization.grid module

class CachedData(data)

Bases: object

property data

property data_nocopy

update_data(data)

class Grid(grid_type=None, top=None, botm=None, idomain=None, lenuni=None, epsg=None, proj4=None, prj=None, xoff=0.0, yoff=0.0, angrot=0.0)

Bases: object

Base class for a structured or unstructured model grid

Parameters

• grid_type (enumeration) – type of model grid (‘structured’, ‘vertex’, ‘unstructured’)

• top (ndarray(float)) – top elevations of cells in topmost layer

• botm (ndarray(float)) – bottom elevations of all cells

• idomain (ndarray(int)) – ibound/idomain value for each cell

• lenuni (ndarray(int)) – model length units

• espg (str, int) – optional espg projection code

• proj4 (str) – optional proj4 projection string code

• prj (str) – optional projection file name path

• xoff (float) – x coordinate of the origin point (lower left corner of model grid) in the spatial reference coordinate system

• yoff (float) – y coordinate of the origin point (lower left corner of model grid) in the spatial reference coordinate system

• angrot (float) – rotation angle of model grid, as it is rotated around the origin point

grid_type

type of model grid (‘structured’, ‘vertex’, ‘unstructured’)

Type

enumeration

top

top elevations of cells in topmost layer

Type

ndarray(float)

botm

bottom elevations of all cells

Type

ndarray(float)

idomain

ibound/idomain value for each cell

Type

ndarray(int)

proj4

spatial reference locates the grid in a coordinate system

Type

proj4 SpatialReference

epsg

spatial reference locates the grid in a coordinate system

Type

epsg SpatialReference

lenuni

modflow lenuni parameter

Type

int

xoffset

x coordinate of the origin point in the spatial reference coordinate system

Type

float

yoffset

y coordinate of the origin point in the spatial reference coordinate system

Type

float

angrot

rotation angle of model grid, as it is rotated around the origin point

Type

float

angrot_radians

rotation angle of model grid, in radians

Type

float

xgrid

returns numpy meshgrid of x edges in reference frame defined by point_type

Type

ndarray

ygrid

returns numpy meshgrid of y edges in reference frame defined by point_type

Type

ndarray

zgrid

returns numpy meshgrid of z edges in reference frame defined by point_type

Type

ndarray

xcenters

returns x coordinate of cell centers

Type

ndarray

ycenters

returns y coordinate of cell centers

Type

ndarray

ycenters

returns z coordinate of cell centers

Type

ndarray

xyzgrid

returns the location of grid edges of all model cells. if the model grid contains spatial reference information, the grid edges are in the coordinate system provided by the spatial reference information. returns a list of three ndarrays for the x, y, and z coordinates

Type

[ndarray, ndarray, ndarray]

xyzcellcenters

returns the cell centers of all model cells in the model grid. if the model grid contains spatial reference information, the cell centers are in the coordinate system provided by the spatial reference information. otherwise the cell centers are based on a 0,0 location for the upper left corner of the model grid. returns a list of three ndarrays for the x, y, and z coordinates

Type

[ndarray, ndarray, ndarray]

get_coords(x, y)

transform point or array of points x, y from model coordinates to spatial coordinates

grid_lines : (point_type=PointType.spatialxyz) : list

returns the model grid lines in a list. each line is returned as a list containing two tuples in the format [(x1,y1), (x2,y2)] where x1,y1 and x2,y2 are the endpoints of the line.

xyvertices : (point_type) : ndarray

1D array of x and y coordinates of cell vertices for whole grid (single layer) in C-style (row-major) order (same as np.ravel())

intersect(x, y, local)

returns the row and column of the grid that the x, y point is in

Notes

Examples

property angrot

property angrot_radians

attribs_from_namfile_header(namefile)

property botm

cross_section_adjust_indicies(k, cbcnt)

Method to get adjusted indicies by layer and confining bed for PlotCrossSection plotting

Parameters

• k (int) – zero based layer number

• cbcnt (int) – confining bed counter

Returns

• tuple ((int, int, int) (adjusted layer, nodeskip layer, node)

• adjustment value based on number of confining beds and the layer)

cross_section_lay_ncpl_ncb(ncb)

Get PlotCrossSection compatible layers, ncpl, and ncb variables

Parameters

ncb (int) – number of confining beds

Returns

tuple

Return type

(int, int, int) layers, ncpl, ncb

cross_section_nodeskip(nlay, xypts)

Get a nodeskip list for PlotCrossSection. This is a correction for UnstructuredGridPlotting

Parameters

• nlay (int) – nlay is nlay + ncb

• xypts (dict) – dictionary of node number and xyvertices of a cross-section

Returns

list

Return type

n-dimensional list of nodes to not plot for each layer

cross_section_set_contour_arrays(plotarray, xcenters, head, elev, projpts)

Method to set countour array centers for rare instances where matplotlib contouring is prefered over trimesh plotting

Parameters

• plotarray (np.ndarray) – array of data for contouring

• xcenters (np.ndarray) – xcenters array

• zcenters (np.ndarray) – zcenters array

• head (np.ndarray) – head array to adjust cell centers location

• elev (np.ndarray) – cell elevation array

• projpts (dict) – dictionary of projected cross sectional vertices

Returns

• tuple ((np.ndarray, np.ndarray, np.ndarray, bool))

• plotarray, xcenter array, ycenter array, and a boolean flag

• for contouring

property cross_section_vertices

property epsg

property extent

classmethod from_binary_grid_file(file_path, verbose=False)

get_coords(x, y)

Given x and y array-like values, apply rotation, scale and offset, to convert them from model coordinates to real-world coordinates.

get_local_coords(x, y)

Given x and y array-like values, apply rotation, scale and offset, to convert them from real-world coordinates to model coordinates.

get_number_plottable_layers(a)

get_plottable_layer_array(plotarray, layer)

get_plottable_layer_shape(layer=None)

Determine the shape that is required in order to plot a 2d array for this grid. For a regular MODFLOW grid, this is (nrow, ncol). For a vertex grid, this is (ncpl,) and for an unstructured grid this is (ncpl[layer],).

Parameters

layer (int) – Has no effect unless grid changes by layer

Returns

shape – required shape of array to plot for a layer

Return type

tuple

get_xcellcenters_for_layer(layer)

get_xvertices_for_layer(layer)

get_ycellcenters_for_layer(layer)

get_yvertices_for_layer(layer)

property grid_lines

property grid_type

property idomain

intersect(x, y, local=False, forgive=False)

property is_complete

property is_valid

property iverts

property laycbd

property lenuni

load_coord_info(namefile=None, reffile='usgs.model.reference')

Attempts to load spatial reference information from the following files (in order): 1) usgs.model.reference 2) NAM file (header comment) 3) defaults

property map_polygons

property ncpl

property nnodes

property nvert

property prj

property proj4

read_usgs_model_reference_file(reffile='usgs.model.reference')

read spatial reference info from the usgs.model.reference file https://water.usgs.gov/ogw/policy/gw-model/modelers-setup.html

remove_confining_beds(array)

Method to remove confining bed layers from an array

Parameters

array (np.ndarray) – array to remove quasi3d confining bed data from. Shape of axis 0 should be (self.lay + ncb) to remove beds

Return type

np.ndarray

saturated_thick(array, mask=None)

Get the saturated thickness for a structured, vertex, or unstructured grid. If the optional array is passed then thickness is returned relative to array values (saturated thickness). Returned values ranges from zero to cell thickness if optional array is passed.

Parameters

• array (ndarray) – array of elevations that will be used to adjust the cell thickness

• mask (float, list, tuple, ndarray) – array values to replace with a nan value.

Returns

thick

Return type

calculated saturated thickness

set_coord_info(xoff=None, yoff=None, angrot=None, epsg=None, proj4=None, merge_coord_info=True)

property shape

property thick

Get the cell thickness for a structured, vertex, or unstructured grid.

Returns

thick

Return type

calculated thickness

property top

property top_botm

property units

property verts

write_shapefile(filename='grid.shp', epsg=None, prj=None)

Write a shapefile of the grid with just the row and column attributes.

property xcellcenters

property xoffset

property xvertices

property xyzcellcenters

property xyzextent

property xyzvertices

property ycellcenters

property yoffset

property yvertices

property zcellcenters

property zvertices