- class imod.wq.BasicTransport(icbund, starting_concentration, porosity=0.35, n_species=1, inactive_concentration=1e+30, minimum_active_thickness=0.01)#
Handles basic tasks that are required by the entire transport model. Among these tasks are definition of the problem, specification of the boundary and initial conditions, determination of the stepsize, preparation of mass balance information, and printout of the simulation results.
icbund (xr.DataArray of int) – is an integer array specifying the boundary condition type (inactive, constant-concentration, or active) for every model cell. For multi-species simulation, ICBUND defines the boundary condition type shared by all species. Note that different species are allowed to have different constant-concentration conditions through an option in the Source and Sink Mixing Package. ICBUND=0, the cell is an inactive concentration cell for all species. Note that no-flow or “dry” cells are automatically converted into inactive concentration cells. Furthermore, active cells in terms of flow can be treated as inactive concentration cells to minimize the area needed for transport simulation, as long as the solute transport is insignificant near those cells. ICBUND<0, the cell is a constant-concentration cell for all species. The starting concentration of each species remains the same at the cell throughout the simulation. (To define different constantconcentration conditions for different species at the same cell location, refer to the Sink/Source Mixing Package.) Also note that unless explicitly defined as a constant-concentration cell, a constant-head cell in the flow model is not treated as a constantconcentration cell. If ICBUND>0, the cell is an active (variable) concentration cell where the concentration value will be calculated.
starting_concentration (float or xr.DataArray of floats) – is the starting concentration (initial condition) at the beginning of the simulation (unit: ML-3) (SCONC). For multispecies simulation, the starting concentration must be specified for all species, one species at a time.
porosity (float, optional) – is the “effective” porosity of the porous medium in a single porosity system (PRSITY). Default value is 0.35.
n_species (int, optional) – is the total number of chemical species included in the current simulation (NCOMP). For single-species simulation, set n_species = 1. Default value is 1.
inactive_concentration (float, optional) – is the value for indicating an inactive concentration cell (ICBUND=0) (CINACT). Even if it is not anticipated to have inactive cells in the model, a value for inactive_concentration still must be submitted. Default value is 1.0e30
minimum_active_thickness (float, optional) – is the minimum saturated thickness in a cell (THKMIN), expressed as the decimal fraction of the model layer thickness, below which the cell is considered inactive. Default value is 0.01 (i.e., 1% of the model layer thickness).
- __init__(icbund, starting_concentration, porosity=0.35, n_species=1, inactive_concentration=1e+30, minimum_active_thickness=0.01)#
__init__(icbund, starting_concentration[, ...])
Loads an imod-wq package from a file (currently only netcdf and zarr are supported).
write_netcdf(directory, pkgname[, ...])
Write to netcdf.