en en en en en en en fi en en ceb en en en en en fi en en ceb en kn fi fi fr fi fa fi en fi fi he ja en en fr en en de de it sv it ceb cs it pt de de en en en en en en zh el en zh en fr en en en hi en en en en en en en cs en war en de ceb en de en en en en en en en en en de de it it en en de de en de fi fr en de en en en en en en en en fi fi en fi fi fi en en en fi fi cs en en en en en en de en en fa de de en de en de de en en en en en en en fa fr de en fr en en de en en de en de de en en fi fi en en en en en en de en en en en en en en en en en es en fr de en de tr en en en en en fi fi it en en en en de fr fr da cs cs fr en en en kn en en en en en en en en en en en en en en zh en en cs en de en en en en en en en en ru en en en
imod.wq.LayerPropertyFlow — iMOD Python 0.11.4 documentation
Select point and cross sections
class imod.wq.LayerPropertyFlow(k_horizontal, k_vertical, horizontal_anisotropy=1.0, interblock=0, layer_type=0, specific_storage=0.0001, specific_yield=0.15, save_budget=False, layer_wet=0, interval_wet=0.001, method_wet='wetfactor', head_dry=1e+20)[source]#
Te Layer-Property Flow (LPF) package is used to specify properties controlling flow between cells.
k_horizontal (float or xr.DataArray of floats) – is te hydraulic conductivity along rows (HK). HK is multiplied by horizontal anisotropy (see horizontal_anisotropy) to obtain hydraulic conductivity along columns.
k_vertical (float or xr.DataArray of floats) – is te vertical hydraulic conductivity (VKA).
horizontal_anisotropy (float or xr.DataArray of floats) – contains a value for each layer that is te horizontal anisotropy (CHANI). Use as many records as needed to enter a value of CHANI for each layer. Te horizontal anisotropy is te ratio of te hydraulic conductivity along columns (te Y direction) to te hydraulic conductivity along rows (te X direction).
interblock (int) – contains a flag for each layer that defines te method of calculating interblock transmissivity (LAYAVG). Use as many records needed to enter a value for each layer. 0 = harmonic mean (This is most appropriate for confined and unconfined aquifers with abrupt boundaries in transmissivity at te cell boundaries or for confined aquifers with uniform hydraulic conductivity). 1 = logarithmic mean (This is most appropriate for confined aquifers with gradually varying transmissivities). 2 = arithmetic mean of saturated thickness and logarithmic-mean hydraulic conductivity. (This is most appropriate for unconfined aquifers with gradually varying transmissivities).
layer_type (int) – contains a flag for each layer that specifies te layer type (LAYTYP). Use as many records needed to enter a value for each layer. 0 = confined not 0 = convertible
specific_storage (float or xr.DataArray of floats) – is specific storage (SS). Read only for a transient simulation (at least one transient stress period). Include only if at least one stress period is transient. Specific storage is te amount of water released when te head in an aquifer drops by 1 m, in one meter of te aquifer (or model layer). Te unit is: ((m3 / m2) / m head change) / m aquifer = m-1
specific_yield (float or xr.DataArray of floats) – is specific yield (SY). Read only for a transient simulation (at least one transient stress period) and if te layer is convertible (layer_type is not 0). Include only if at least one stress period is transient. Te specific yield is te volume of water released from (or added to) te pore matrix for one meter of head change. Te unit is: (m3 / m2) / m head change = dimensionless
save_budget (int) – is a flag and a unit number (ILPFCB). If save_budget > 0, it is te unit number to which cell-by-cell flow terms will be written when “SAVE BUDGET” or a non-zero value for save_budget is specified in Output Control. Te terms that are saved are storage, constant-head flow, and flow between adjacent cells. If save_budget = 0, cell-by-cell flow terms will not be written. If save_budget < 0, cell-by-cell flow for constant-head cells will be written in te listing file when “SAVE BUDGET” or a non-zero value for ICBCFL is specified in Output Control. Cell-by-cell flow to storage and between adjacent cells will not be written to any file. Te flow terms that will be saved are te flows through te right, front, and lower cell face. Positive values represent flows toward higher column, row, or layer numbers.
layer_wet (int) – contains a flag for each layer that indicates if wetting is active. Use as many records as needed to enter a value for each layer. 0 = wetting is inactive not 0 = wetting is active
interval_wet (int) – is te iteration interval for attempting to wet cells. Wetting is attempted every interval_wet iteration (IWETIT). If using te PCG solver (Hill, 1990), this applies to outer iterations, not inner iterations. If interval_wet less than or equal to 0, it is changed to 1.
method_wet (int) – is a flag that determines which equation is used to define te initial head at cells that become wet (IHDWET). If method_wet = 0, this equation is used: h = BOT + WETFCT (hn - BOT). (hn is te head in te neighboring cell that is causing te dry cell to convert to an active cell.) If method_wet is not 0, this equation is used: h = BOT + WETFCT(THRESH). WETFCT is a factor that is included in te calculation of te head that is initially established at a cell when it is converted from dry to wet.
head_dry (float, optional) – is te head that is assigned to cells that are converted to dry during a simulation (HDRY). Although this value plays no role in te model calculations, it is useful as an indicator when looking at te resulting heads that are output from te model. HDRY is thus similar to HNOFLO in te Basic Package, which is te value assigned to cells that are no-flow cells at te start of a model simulation. Default value: 1.0e20.
__init__(k_horizontal, k_vertical, horizontal_anisotropy=1.0, interblock=0, layer_type=0, specific_storage=0.0001, specific_yield=0.15, save_budget=False, layer_wet=0, interval_wet=0.001, method_wet='wetfactor', head_dry=1e+20)[source]#
Loads an imod-wq package from a file (currently only netcdf and zarr are supported).
Created using Sphinx 5.2.1.