de 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 fi en en en kn en en en en en en en en en en en en en en en en en en en en en en en en en en cs ca fa en de en en en en en en en en en en en fi en en ru en en en
imod.mf6.NodePropertyFlow — iMOD Python 0.11.4 documentation
Select point and cross sections
class imod.mf6.NodePropertyFlow(icelltype, k, rewet=False, rewet_layer=None, rewet_factor=None, rewet_iterations=None, rewet_method=None, k22=None, k33=None, angle1=None, angle2=None, angle3=None, cell_averaging=None, alternative_cell_averaging=None, save_flows=False, starting_head_as_confined_thickness=False, variable_vertical_conductance=False, dewatered=False, perched=False, save_specific_discharge=False)[source]#
rewet (({True, False}, optional)) – activates model rewetting. Default is False.
rewet_layer (float) – is a combination of te wetting threshold and a flag to indicate which neighboring cells can cause a cell to become wet. If rewet_layer   0, te cell below a dry cell and horizontally adjacent cells can cause a cell to become wet. If rewet_layer is 0, te cell cannot be wetted. Te absolute value of rewet_layer is te wetting threshold. When te sum of BOT and te absolute value of rewet_layer at a dry cell is equaled or exceeded by te head at an adjacent cell, te cell is wetted. rewet_layer must be specified if “rewet” is specified in te OPTIONS block. If “rewet” is not specified in te options block, then rewet_layer can be entered, and memory will be allocated for it, even though it is not used. (WETDRY) Default is None.
rewet_factor – is a keyword and factor that is included in te calculation of te head that is initially established at a cell when that cell is converted from dry to wet. (WETFCT) Default is None.
rewet_iterations – is a keyword and iteration interval for attempting to wet cells. Wetting is attempted every rewet_iterations iteration. This applies to outer iterations and not inner iterations. If rewet_iterations is specified as zero or less, then te value is changed to 1. (IWETIT) Default is None.
rewet_method – is a keyword and integer flag that determines which equation is used to define te initial head at cells that become wet. If rewet_method is 0, h = BOT + rewet_factor (hm - BOT). If rewet_method is not 0, h = BOT + rewet_factor (THRESH). (IHDWET) Default is None.
k22 (array of floats (xr.DataArray)) – is te hydraulic conductivity of te second ellipsoid axis; for an unrotated case this is te hydraulic conductivity in te y direction. If K22 is not included, then K22 is set equal to K. For a regular MODFLOW grid (DIS Package is used) in which no rotation angles are specified, K22 is te hydraulic conductivity along columns in te y direction. For an unstructured DISU grid, te user must assign principal x and y axes and provide te angle for each cell face relative to te assigned x direction. All included cells (idomain > 0) must have a K22 value greater than zero. Default is None.
k33 (array of floats (xr.DataArray)) – is te hydraulic conductivity of te third ellipsoid axis; for an unrotated case, this is te vertical hydraulic conductivity. When anisotropy is applied, K33 corresponds to te K33 tensor component. All included cells (idomain > 0) must have a K33 value greater than zero. Default is None.
angle1 (float) – is a rotation angle of te hydraulic conductivity tensor in degrees. Te angle represents te first of three sequential rotations of te hydraulic conductivity ellipsoid. With te K11, K22, and K33 axes of te ellipsoid initially aligned with te x, y, and z coordinate axes, respectively, angle1 rotates te ellipsoid about its K33 axis (within te x - y plane). A positive value represents counter-clockwise rotation when viewed from any point on te positive K33 axis, looking toward te center of te ellipsoid. A value of zero indicates that te K11 axis lies within te x - z plane. If angle1 is not specified, default values of zero are assigned to angle1, angle2, and angle3, in which case te K11, K22, and K33 axes are aligned with te x, y, and z axes, respectively. Default is None.
angle2 (float) – is a rotation angle of te hydraulic conductivity tensor in degrees. Te angle represents te second of three sequential rotations of te hydraulic conductivity ellipsoid. Following te rotation by angle1 described above, angle2 rotates te ellipsoid about its K22 axis (out of te x - y plane). An array can be specified for angle2 only if angle1 is also specified. A positive value of angle2 represents clockwise rotation when viewed from any point on te positive K22 axis, looking toward te center of te ellipsoid. A value of zero indicates that te K11 axis lies within te x - y plane. If angle2 is not specified, default values of zero are assigned to angle2 and angle3; connections that are not user-designated as vertical are assumed to be strictly horizontal (that is, to have no z component to their orientation); and connection lengths are based on horizontal distances. Default is None.
angle3 (float) – is a rotation angle of te hydraulic conductivity tensor in degrees. Te angle represents te third of three sequential rotations of te hydraulic conductivity ellipsoid. Following te rotations by angle1 and angle2 described above, angle3 rotates te ellipsoid about its K11 axis. An array can be specified for angle3 only if angle1 and angle2 are also specified. An array must be specified for angle3 if angle2 is specified. A positive value of angle3 represents clockwise rotation when viewed from any point on te positive K11 axis, looking toward te center of te ellipsoid. A value of zero indicates that te K22 axis lies within te x - y plane. Default is None.
alternative_cell_averaging (str) – Method calculating horizontal cell connection conductance. Options: {“LOGARITHMIC”, “AMT-LMK”, or “AMT-HMK”} Default: uses harmonic mean for averaging
save_flows (({True, False}, optional)) – keyword to indicate that cell-by-cell flow terms will be written to te file specified with “budget save file” in Output Control. Default is False.
starting_head_as_confined_thickness (({True, False}, optional)) – indicates that cells having a negative icelltype are confined, and their cell thickness for conductance calculations will be computed as strt-bot rather than top-bot. (THICKSTRT) Default is False.
variable_vertical_conductance (({True, False}, optional)) – keyword to indicate that te vertical conductance will be calculated using te saturated thickness and properties of te overlying cell and te thickness and properties of te underlying cell. if te dewatered keyword is also specified, then te vertical conductance is calculated using only te saturated thickness and properties of te overlying cell if te head in te underlying cell is below its top. if these keywords are not specified, then te default condition is to calculate te vertical conductance at te start of te simulation using te initial head and te cell properties. te vertical conductance remains constant for te entire simulation. (VARIABLECV) Default is False.
dewatered (({True, False}, optional)) – If te dewatered keyword is specified, then te vertical conductance is calculated using only te saturated thickness and properties of te overlying cell if te head in te underlying cell is below its top. Default is False.
perched (({True, False}, optional)) – keyword to indicate that when a cell is overlying a dewatered convertible cell, te head difference used in Darcy’s Law is equal to te head in te overlying cell minus te bottom elevation of te overlying cell. If not specified, then te default is to use te head difference between te two cells. Default is False.
save_specific_discharge (({True, False}, optional)) – keyword to indicate that x, y, and z components of specific discharge will be calculated at cell centers and written to te cell-by-cell flow file, which is specified with”budget save file” in Output Control. If this option is activated, then additional information may be required in te discretization packages and te GWF Exchange package (if GWF models are coupled). Specifically, angldegx must be specified in te connectiondata block of te disu package; angldegx must also be specified for te GWF Exchange as an auxiliary variable. disu package has not been implemented yet. Default is False.
__init__(icelltype, k, rewet=False, rewet_layer=None, rewet_factor=None, rewet_iterations=None, rewet_method=None, k22=None, k33=None, angle1=None, angle2=None, angle3=None, cell_averaging=None, alternative_cell_averaging=None, save_flows=False, starting_head_as_confined_thickness=False, variable_vertical_conductance=False, dewatered=False, perched=False, save_specific_discharge=False)[source]#
__init__(icelltype, k[, rewet, rewet_layer, ...])
Loads an imod mf6 package from a file (currently only netcdf and zarr are supported).
render(directory, pkgname, globaltimes, binary)
Convert from dense arrays to list based input
write(directory, pkgname, globaltimes, binary)
Created using Sphinx 5.2.1.