From ac2626cc6a7d4cb09ab72054b453544a933e0f0e Mon Sep 17 00:00:00 2001 From: wpbonelli Date: Fri, 30 Oct 2026 14:52:42 -0700 Subject: [PATCH] chore(mf6): regenerate mf6 modules For the MF6.8.2 release Co-Authored-By: Claude Opus 5.5 --- flopy/mf6/modflow/mfgweest.py | 5 + flopy/mf6/modflow/mfgwelke.py | 35 ++++--- flopy/mf6/modflow/mfgwemwe.py | 23 ++++- flopy/mf6/modflow/mfgwesfe.py | 35 ++++--- flopy/mf6/modflow/mfgweuze.py | 27 ++++-- flopy/mf6/modflow/mfgwfdisu.py | 26 +++++- flopy/mf6/modflow/mfgwflak.py | 165 ++++++++++++++++++++++++--------- flopy/mf6/modflow/mfgwfmaw.py | 19 ++++ flopy/mf6/modflow/mfgwfmvr.py | 2 +- flopy/mf6/modflow/mfgwfsfr.py | 90 +++++++++++++----- flopy/mf6/modflow/mfgwfuzf.py | 38 +++++--- flopy/mf6/modflow/mfgwtfmi.py | 12 +-- flopy/mf6/modflow/mfgwtlkt.py | 34 ++++--- flopy/mf6/modflow/mfgwtmwt.py | 22 ++++- flopy/mf6/modflow/mfgwtsft.py | 34 ++++--- flopy/mf6/modflow/mfgwtuzt.py | 26 ++++-- flopy/mf6/modflow/mfmvr.py | 2 +- flopy/mf6/modflow/mfnam.py | 1 - flopy/mf6/modflow/mfprtprp.py | 12 ++- flopy/mf6/modflow/mfutlncf.py | 39 +++++--- flopy/mf6/modflow/mfutlspc.py | 10 +- flopy/mf6/modflow/mfutltas.py | 2 +- flopy/mf6/modflow/mfutlts.py | 24 +++-- flopy/mf6/modflow/mfutltvk.py | 15 ++- flopy/mf6/modflow/mfutltvs.py | 10 +- 25 files changed, 501 insertions(+), 207 deletions(-) diff --git a/flopy/mf6/modflow/mfgweest.py b/flopy/mf6/modflow/mfgweest.py index b7e021a94c..82917569bc 100644 --- a/flopy/mf6/modflow/mfgweest.py +++ b/flopy/mf6/modflow/mfgweest.py @@ -159,6 +159,7 @@ class ModflowGweest(MFPackage): "shape (nodes)", "reader readarray", "layered true", + "netcdf true", ], [ "block griddata", @@ -167,6 +168,7 @@ class ModflowGweest(MFPackage): "shape (nodes)", "reader readarray", "layered true", + "netcdf true", "optional true", ], [ @@ -176,6 +178,7 @@ class ModflowGweest(MFPackage): "shape (nodes)", "reader readarray", "layered true", + "netcdf true", "optional true", ], [ @@ -185,6 +188,7 @@ class ModflowGweest(MFPackage): "shape (nodes)", "reader readarray", "layered true", + "netcdf true", "mf6internal cps", ], [ @@ -194,6 +198,7 @@ class ModflowGweest(MFPackage): "shape (nodes)", "reader readarray", "layered true", + "netcdf true", "mf6internal rhos", ], ] diff --git a/flopy/mf6/modflow/mfgwelke.py b/flopy/mf6/modflow/mfgwelke.py index 691058ba25..1305cdb8d0 100644 --- a/flopy/mf6/modflow/mfgwelke.py +++ b/flopy/mf6/modflow/mfgwelke.py @@ -131,7 +131,7 @@ class ModflowGwelke(MFPackage): package supports a 'withdrawal' flow term. if this withdrawal term is active, then water will be withdrawn from the lake at the calculated temperature of the lake. - * temperature : string + * temperature : double precision real or character value that defines the temperature for the lake. the specified temperature is only applied if the lake is a constant temperature lake. if the options block includes a timeseriesfile entry (see the 'time- @@ -144,23 +144,23 @@ class ModflowGwelke(MFPackage): fluxes into or out of the lake and the inactive value will be written for the lake temperature. if a lake is constant, then the temperature for the lake will be fixed at the user specified value. - * rainfall : string + * rainfall : double precision real or character value that defines the rainfall temperature for the lake. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * evaporation : string + * evaporation : double precision use of the evaporation keyword is allowed in the lke package; however, the specified value is not currently used in lke calculations. instead, the latent heat of evaporation is multiplied by the simulated evaporation rate for determining the thermal energy lost from a stream reach. - * runoff : string + * runoff : double precision real or character value that defines the temperature of runoff for the lake. users are free to use whatever temperature scale they want, which might include negative temperatures. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * ext_inflow : string + * ext_inflow : double precision real or character value that defines the temperature of external inflow for the lake. users are free to use whatever temperature scale they want, which might include negative temperatures. if the options block includes a timeseriesfile @@ -218,6 +218,7 @@ class ModflowGwelke(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -234,6 +235,7 @@ class ModflowGwelke(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -249,6 +251,7 @@ class ModflowGwelke(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -256,6 +259,7 @@ class ModflowGwelke(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -263,6 +267,7 @@ class ModflowGwelke(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -270,6 +275,7 @@ class ModflowGwelke(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -279,6 +285,7 @@ class ModflowGwelke(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal temp_filerec", ], [ "block options", @@ -309,6 +316,7 @@ class ModflowGwelke(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -349,6 +357,7 @@ class ModflowGwelke(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -451,12 +460,13 @@ class ModflowGwelke(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray lakeno strt ktf rbthcnd aux boundname", - "shape (maxbound)", + "shape (nlakes)", "reader urword", ], [ @@ -468,6 +478,7 @@ class ModflowGwelke(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -539,6 +550,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name lakeno", + "mf6internal ifno", "type integer", "shape", "tagged false", @@ -567,7 +579,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name temperature", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -577,7 +589,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -587,7 +599,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -597,7 +609,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -607,7 +619,7 @@ class ModflowGwelke(MFPackage): [ "block period", "name ext-inflow", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -630,6 +642,7 @@ class ModflowGwelke(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwemwe.py b/flopy/mf6/modflow/mfgwemwe.py index ceddcec313..a7693cc6b1 100644 --- a/flopy/mf6/modflow/mfgwemwe.py +++ b/flopy/mf6/modflow/mfgwemwe.py @@ -131,7 +131,7 @@ class ModflowGwemwe(MFPackage): aquifer well package supports a 'withdrawal' flow term. if this withdrawal term is active, then water will be withdrawn from the well at the calculated temperature of the well. - * temperature : string + * temperature : double precision real or character value that defines the temperature for the well. the specified temperature is only applied if the well is a constant temperature well. if the options block includes a timeseriesfile entry (see the 'time- @@ -144,7 +144,7 @@ class ModflowGwemwe(MFPackage): mass fluxes into or out of the well and the inactive value will be written for the well temperature. if a well is constant, then the temperature for the well will be fixed at the user specified value. - * rate : string + * rate : double precision real or character value that defines the injection temperature :math:`(e.g.,:^{circ}c:or:^{circ}f)` for the well. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values @@ -202,6 +202,7 @@ class ModflowGwemwe(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -218,6 +219,7 @@ class ModflowGwemwe(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -233,6 +235,7 @@ class ModflowGwemwe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -240,6 +243,7 @@ class ModflowGwemwe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -247,6 +251,7 @@ class ModflowGwemwe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -254,6 +259,7 @@ class ModflowGwemwe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -263,6 +269,7 @@ class ModflowGwemwe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal temp_filerec", ], [ "block options", @@ -293,6 +300,7 @@ class ModflowGwemwe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -333,6 +341,7 @@ class ModflowGwemwe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -435,12 +444,13 @@ class ModflowGwemwe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray mawno strt ktf fthk aux boundname", - "shape (maxbound)", + "shape (nmawwells)", "reader urword", ], [ @@ -452,6 +462,7 @@ class ModflowGwemwe(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -523,6 +534,7 @@ class ModflowGwemwe(MFPackage): [ "block period", "name mawno", + "mf6internal ifno", "type integer", "shape", "tagged false", @@ -551,7 +563,7 @@ class ModflowGwemwe(MFPackage): [ "block period", "name temperature", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -561,7 +573,7 @@ class ModflowGwemwe(MFPackage): [ "block period", "name rate", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -584,6 +596,7 @@ class ModflowGwemwe(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwesfe.py b/flopy/mf6/modflow/mfgwesfe.py index 8323883877..c41da75877 100644 --- a/flopy/mf6/modflow/mfgwesfe.py +++ b/flopy/mf6/modflow/mfgwesfe.py @@ -131,7 +131,7 @@ class ModflowGwesfe(MFPackage): temperatures cannot be specified for all flow terms. for example, the streamflow package supports a 'diversion' flow term. diversion water will be routed using the calculated temperature of the reach. - * temperature : string + * temperature : double precision real or character value that defines the temperature for the reach. the specified temperature is only applied if the reach is a constant temperature reach. if the options block includes a timeseriesfile entry (see the 'time- @@ -144,25 +144,25 @@ class ModflowGwesfe(MFPackage): fluxes into or out of the reach and the inactive value will be written for the reach temperature. if a reach is constant, then the temperature for the reach will be fixed at the user specified value. - * rainfall : string + * rainfall : double precision real or character value that defines the rainfall temperature :math:`(e.g.,:^{circ}c:or:^{circ}f)` for the reach. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * evaporation : string + * evaporation : double precision use of the evaporation keyword is allowed in the sfe package; however, the specified value is not currently used in sfe calculations. instead, the latent heat of evaporation is multiplied by the simulated evaporation rate for determining the thermal energy lost from a stream reach. - * runoff : string + * runoff : double precision real or character value that defines the temperature of runoff :math:`(e.g.,:^{circ}c:or:^{circ}f)` for the reach. users are free to use whatever temperature scale they want, which might include negative temperatures. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * inflow : string + * inflow : double precision real or character value that defines the temperature of inflow :math:`(e.g.,:^{circ}c:or:^{circ}f)` for the reach. users are free to use whatever temperature scale they want, which might include negative @@ -223,6 +223,7 @@ class ModflowGwesfe(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -239,6 +240,7 @@ class ModflowGwesfe(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -254,6 +256,7 @@ class ModflowGwesfe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -261,6 +264,7 @@ class ModflowGwesfe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -268,6 +272,7 @@ class ModflowGwesfe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -275,6 +280,7 @@ class ModflowGwesfe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -284,6 +290,7 @@ class ModflowGwesfe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal temp_filerec", ], [ "block options", @@ -314,6 +321,7 @@ class ModflowGwesfe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -354,6 +362,7 @@ class ModflowGwesfe(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -456,12 +465,13 @@ class ModflowGwesfe(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray rno strt ktf rbthcnd aux boundname", - "shape (maxbound)", + "shape (nreaches)", "reader urword", ], [ @@ -473,6 +483,7 @@ class ModflowGwesfe(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -544,6 +555,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name rno", + "mf6internal ifno", "type integer", "shape", "tagged false", @@ -572,7 +584,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name temperature", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -582,7 +594,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -592,7 +604,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -602,7 +614,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -612,7 +624,7 @@ class ModflowGwesfe(MFPackage): [ "block period", "name inflow", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -635,6 +647,7 @@ class ModflowGwesfe(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgweuze.py b/flopy/mf6/modflow/mfgweuze.py index bf61bdc0e5..d851afc852 100644 --- a/flopy/mf6/modflow/mfgweuze.py +++ b/flopy/mf6/modflow/mfgweuze.py @@ -124,7 +124,7 @@ class ModflowGweuze(MFPackage): infiltration, uzet, and auxiliary. these settings are used to assign the temperature associated with the corresponding flow terms. temperatures cannot be specified for all flow terms. - * temperature : string + * temperature : double precision real or character value that defines the temperature for the unsaturated zone flow cell. the specified temperature is only applied if the unsaturated zone flow cell is a constant temperature cell. if the options block includes a @@ -138,13 +138,13 @@ class ModflowGweuze(MFPackage): energy fluxes into or out of the uzf cell and the inactive value will be written for the uzf cell temperature. if a uzf cell is constant, then the temperature for the uzf cell will be fixed at the user specified value. - * infiltration : string + * infiltration : double precision real or character value that defines the temperature of the infiltration :math:`(e.g.,:^{circ}c:or:^{circ}f)` for the uzf cell. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * uzet : string + * uzet : double precision real or character value that states what fraction of the simulated unsaturated zone evapotranspiration is associated with evaporation. the evaporative losses from the unsaturated zone moisture content will have an evaporative cooling @@ -205,6 +205,7 @@ class ModflowGweuze(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -221,6 +222,7 @@ class ModflowGweuze(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -236,6 +238,7 @@ class ModflowGweuze(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -243,6 +246,7 @@ class ModflowGweuze(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -250,6 +254,7 @@ class ModflowGweuze(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -257,6 +262,7 @@ class ModflowGweuze(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -266,6 +272,7 @@ class ModflowGweuze(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal temp_filerec", ], [ "block options", @@ -296,6 +303,7 @@ class ModflowGweuze(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -336,6 +344,7 @@ class ModflowGweuze(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -438,12 +447,13 @@ class ModflowGweuze(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray uzfno strt aux boundname", - "shape (maxbound)", + "shape (nuzfcells)", "reader urword", ], [ @@ -455,6 +465,7 @@ class ModflowGweuze(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -508,6 +519,7 @@ class ModflowGweuze(MFPackage): [ "block period", "name uzfno", + "mf6internal ifno", "type integer", "shape", "tagged false", @@ -536,7 +548,7 @@ class ModflowGweuze(MFPackage): [ "block period", "name temperature", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -546,7 +558,7 @@ class ModflowGweuze(MFPackage): [ "block period", "name infiltration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -556,7 +568,7 @@ class ModflowGweuze(MFPackage): [ "block period", "name uzet", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -579,6 +591,7 @@ class ModflowGweuze(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwfdisu.py b/flopy/mf6/modflow/mfgwfdisu.py index 8de18ce067..9801d73f67 100644 --- a/flopy/mf6/modflow/mfgwfdisu.py +++ b/flopy/mf6/modflow/mfgwfdisu.py @@ -143,10 +143,28 @@ class ModflowGwfdisu(MFPackage): only needed if horizontal anisotropy is specified in the npf package, if the xt3d option is used in the npf package, or if the save_specific_discharge option is specified in the npf package. angldegx does not need to be specified - if these conditions are not met. angldegx is of size nja; values specified for - vertical connections and for the diagonal position are not used. note that - angldegx is read in degrees, which is different from modflow-usg, which reads a - similar variable (anglex) in radians. + otherwise. angldegx is of size nja; values specified for vertical connections + and for the diagonal position are not used. note that angldegx is read in + degrees, which is different from modflow-usg, which reads a similar variable + (anglex) in radians. only the value specified for the connection from the + lower to the higher cell number is used, but the values for both directions of + every horizontal connection between active cells are checked. because angldegx + is the outward normal of the shared face, the value for the connection from + cell n to cell m must differ by 180 degrees (within 0.01 degrees) from the + value for the connection from cell m to cell n. if vertices and cell2d are + specified, angldegx must be expressed in the same coordinate system as the + vertex coordinates; the centers of two connected cells must lie on opposite + sides of their shared face, so angldegx must be within 90 degrees of the + direction from the cell center to the center of the connected cell, and, for + cells that share exactly two vertices, angldegx must agree within 1 degree with + the outward normal of the face computed from those vertices. a warning is + issued for each of these three conditions that is not met, and the simulation + terminates with an error if any of them is not met and the npf package requires + angldegx to be specified (xt3d, k22, or save_specific_discharge). a warning is + also issued if angldegx deviates by more than 45 degrees from the direction + between the centers of the connected cells; this is not an error, but such + connections depart from the assumptions of the control-volume finite-difference + method and may warrant a closer look at the grid. vertices : [(iv, xv, yv)] * iv : integer is the vertex number. Records in the VERTICES block must be listed in diff --git a/flopy/mf6/modflow/mfgwflak.py b/flopy/mf6/modflow/mfgwflak.py index fab618f1e8..a6cd4e6599 100644 --- a/flopy/mf6/modflow/mfgwflak.py +++ b/flopy/mf6/modflow/mfgwflak.py @@ -84,11 +84,17 @@ class ModflowGwflak(MFPackage): water mover (mvr) package. when the mover option is specified, additional memory is allocated within the package to store the available, provided, and received water. - dev_force_fallback : keyword - keyword that solves every active lake with the substitution fallback under the - implicit formulation, instead of only those lakes that are too poorly - conditioned to solve as a matrix unknown. this development option is not - supported. + dev_force_legacy : keyword + keyword that solves every active lake with the default (legacy) substitution + solver under the implicit formulation, instead of only those lakes that are too + poorly conditioned to solve as a matrix unknown. cannot be specified with + dev_force_legacy_lake. this development option exists to test that path and is + not supported. + dev_force_legacy_lake : integer + integer lake number that is solved with the default (legacy) substitution + solver under the implicit formulation, while the other lakes are solved as + matrix unknowns. cannot be specified with dev_force_legacy. this development + option exists to test that path and is not supported. dev_groundwater_head_conductance : keyword keyword that calculates the conductance for horizontal connections using the groundwater head instead of the lake stage. this development option is not @@ -114,13 +120,16 @@ class ModflowGwflak(MFPackage): the implicit option is active because the implicit formulation solves the lake stage through its lakebed seepage; the program will terminate with an error if a lake with no groundwater connections is specified with the implicit option. - if a lake stage becomes too poorly conditioned to solve as a matrix unknown, - which can happen for a weakly connected or disconnected lake, that lake is - solved automatically by the default substitution solver instead, so the - implicit formulation remains as robust as the default formulation. the newton - formulation is recommended for lakes that may become perched (disconnected from - the aquifer), because results obtained with the standard rewetting (rewet) - capability can otherwise differ between the two approaches. + the implicit formulation does not yet include the density terms of the buoyancy + (buy) package, and the program will terminate with an error if the implicit + option is specified and the buy package is active. if a lake stage becomes too + poorly conditioned to solve as a matrix unknown, which can happen for a weakly + connected or disconnected lake, that lake is solved automatically by the + default substitution solver instead, so the implicit formulation remains as + robust as the default formulation. the newton formulation is recommended for + lakes that may become perched (disconnected from the aquifer), because results + obtained with the standard rewetting (rewet) capability can otherwise differ + between the two approaches. maximum_iterations : integer integer value that defines the maximum number of newton-raphson iterations allowed for a lake. by default, maximum_iterations is equal to 100. @@ -154,6 +163,9 @@ class ModflowGwflak(MFPackage): and 100.0 when using length units (length_units) of feet, meters, or centimeters in the simulation, respectively. length_conversion does not need to be specified if no lake outlets are specified or length_units are meters. + dev_no_final_check : keyword + flag disabling the final convergence check of the change in lake stages. by + default this check is enabled. nlakes : integer value specifying the number of lakes that will be simulated for all stress periods. @@ -297,26 +309,26 @@ class ModflowGwflak(MFPackage): is defined as a specified flow. MANNING--character keyword to indicate the outlet is defined using Manning's equation. WEIR--character keyword to indicate the outlet is defined using a sharp weir equation. - * invert : string + * invert : double precision real or character value that defines the invert elevation for the lake outlet. A specified INVERT value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is not SPECIFIED. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * width : string + * width : double precision real or character value that defines the width of the lake outlet. A specified WIDTH value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is not SPECIFIED. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * rough : string + * rough : double precision real value that defines the roughness coefficient for the lake outlet. Any value can be specified if COUTTYPE is not MANNING. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * slope : string + * slope : double precision real or character value that defines the bed slope for the lake outlet. A specified SLOPE value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is MANNING. If the Options block includes a TIMESERIESFILE entry (see @@ -336,32 +348,32 @@ class ModflowGwflak(MFPackage): status, stage, rainfall, evaporation, runoff, inflow, withdrawal, and auxiliary. keywords for outlet settings include rate, invert, width, slope, and rough. - * stage : string + * stage : double precision real or character value that defines the stage for the lake. the specified stage is only applied if the lake is a constant stage lake. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * invert : string + * invert : double precision real or character value that defines the invert elevation for the lake outlet. a specified invert value is only used for active lakes if couttype for lake outlet outletno is not specified. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * width : string + * width : double precision real or character value that defines the width of the lake outlet. a specified width value is only used for active lakes if couttype for lake outlet outletno is not specified. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * rough : string + * rough : double precision real value that defines the roughness coefficient for the lake outlet. any value can be specified if couttype is not manning. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * slope : string + * slope : double precision real or character value that defines the bed slope for the lake outlet. a specified slope value is only used for active lakes if couttype for lake outlet outletno is manning. if the options block includes a timeseriesfile entry (see @@ -370,38 +382,38 @@ class ModflowGwflak(MFPackage): * status : string keyword option to define lake status. status can be active, inactive, or constant. by default, status is active. - * rainfall : string + * rainfall : double precision real or character value that defines the rainfall rate :math:`(lt^{-1})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * evaporation : string + * evaporation : double precision real or character value that defines the maximum evaporation rate :math:`(lt^{-1})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * runoff : string + * runoff : double precision real or character value that defines the runoff rate :math:`(l^3 t^{-1})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * inflow : string + * inflow : double precision real or character value that defines the volumetric inflow rate :math:`(l^3 t^{-1})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. by default, inflow rates are zero for each lake. - * withdrawal : string + * withdrawal : double precision real or character value that defines the maximum withdrawal rate :math:`(l^3 t^{-1})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * rate : string + * rate : double precision real or character value that defines the extraction rate for the lake outflow. a positive value indicates inflow and a negative value indicates outflow from the lake. rate only applies to outlets associated with active lakes (status is @@ -485,6 +497,7 @@ class ModflowGwflak(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -499,6 +512,7 @@ class ModflowGwflak(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -506,6 +520,7 @@ class ModflowGwflak(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -519,6 +534,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name stage", + "mf6internal options_stage", "type keyword", "shape", "in_record true", @@ -540,6 +556,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name budget_filerecord", + "mf6internal budfilerec", "type record budget fileout budgetfile", "shape", "reader urword", @@ -580,6 +597,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name budgetcsv_filerecord", + "mf6internal budcsvfilerec", "type record budgetcsv fileout budgetcsvfile", "shape", "reader urword", @@ -610,6 +628,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name package_convergence_filerecord", + "mf6internal pkgconvfilerec", "type record package_convergence fileout package_convergence_filename", "shape", "reader urword", @@ -625,6 +644,7 @@ class ModflowGwflak(MFPackage): "reader urword", "tagged true", "optional false", + "mf6internal pkgconv", ], [ "block options", @@ -635,6 +655,7 @@ class ModflowGwflak(MFPackage): "reader urword", "tagged false", "optional false", + "mf6internal pkgconvfname", ], [ "block options", @@ -661,6 +682,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name filein", + "mf6internal options_filein", "type keyword", "shape", "in_record true", @@ -720,14 +742,24 @@ class ModflowGwflak(MFPackage): ], [ "block options", - "name dev_force_fallback", + "name dev_force_legacy", + "mf6internal iforceleg", "type keyword", "reader urword", "optional true", ], + [ + "block options", + "name dev_force_legacy_lake", + "mf6internal iforceleglak", + "type integer", + "reader urword", + "optional true", + ], [ "block options", "name dev_groundwater_head_conductance", + "mf6internal igwhcopt", "type keyword", "reader urword", "optional true", @@ -735,6 +767,7 @@ class ModflowGwflak(MFPackage): [ "block options", "name dev_maximum_outlet_depth", + "mf6internal outdmax", "type double precision", "reader urword", "optional true", @@ -745,10 +778,12 @@ class ModflowGwflak(MFPackage): "type double precision", "reader urword", "optional true", + "mf6internal surfdep", ], [ "block options", "name implicit", + "mf6internal iimplicit", "type keyword", "reader urword", "optional true", @@ -759,6 +794,7 @@ class ModflowGwflak(MFPackage): "type integer", "reader urword", "optional true", + "mf6internal maxlakit", ], [ "block options", @@ -766,6 +802,7 @@ class ModflowGwflak(MFPackage): "type double precision", "reader urword", "optional true", + "mf6internal dmaxchg", ], [ "block options", @@ -773,6 +810,7 @@ class ModflowGwflak(MFPackage): "type double precision", "reader urword", "optional true", + "mf6internal convtime", ], [ "block options", @@ -780,6 +818,15 @@ class ModflowGwflak(MFPackage): "type double precision", "reader urword", "optional true", + "mf6internal convlength", + ], + [ + "block options", + "name dev_no_final_check", + "type keyword", + "reader urword", + "optional true", + "mf6internal iconvchk", ], [ "block dimensions", @@ -806,7 +853,7 @@ class ModflowGwflak(MFPackage): "block packagedata", "name packagedata", "type recarray ifno strt nlakeconn aux boundname", - "shape (maxbound)", + "shape (nlakes)", "reader urword", ], [ @@ -818,6 +865,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -874,6 +922,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal conndata_ifno", ], [ "block connectiondata", @@ -954,6 +1003,7 @@ class ModflowGwflak(MFPackage): "type recarray ifno tab6 filein tab6_filename", "shape (ntables)", "reader urword", + "optional true", ], [ "block tables", @@ -964,6 +1014,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal tables_ifno", ], [ "block tables", @@ -978,6 +1029,7 @@ class ModflowGwflak(MFPackage): [ "block tables", "name filein", + "mf6internal tables_filein", "type keyword", "shape", "in_record true", @@ -988,6 +1040,7 @@ class ModflowGwflak(MFPackage): [ "block tables", "name tab6_filename", + "mf6internal laktab6_filename", "type string", "preserve_case true", "in_record true", @@ -1001,6 +1054,7 @@ class ModflowGwflak(MFPackage): "type recarray outletno lakein lakeout couttype invert width rough slope", "shape (noutlets)", "reader urword", + "optional true", ], [ "block outlets", @@ -1011,6 +1065,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal outlets_outletno", ], [ "block outlets", @@ -1050,6 +1105,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "time_series true", + "mf6internal outlets_invert", ], [ "block outlets", @@ -1060,6 +1116,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "time_series true", + "mf6internal outlets_width", ], [ "block outlets", @@ -1070,6 +1127,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "time_series true", + "mf6internal outlets_rough", ], [ "block outlets", @@ -1080,6 +1138,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "time_series true", + "mf6internal outlets_slope", ], [ "block period", @@ -1109,6 +1168,7 @@ class ModflowGwflak(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal ifno", ], [ "block period", @@ -1131,8 +1191,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name stage", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "time_series true", @@ -1141,8 +1202,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "reader urword", @@ -1151,8 +1213,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "reader urword", @@ -1161,8 +1224,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "reader urword", @@ -1171,8 +1235,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name inflow", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "reader urword", @@ -1181,8 +1246,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name withdrawal", - "type string", + "type double precision", "shape", + "mf6dimension (nlakes)", "tagged true", "in_record true", "reader urword", @@ -1191,8 +1257,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name rate", - "type string", + "type double precision", "shape", + "mf6dimension (noutlets)", "tagged true", "in_record true", "reader urword", @@ -1201,8 +1268,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name invert", - "type string", + "type double precision", "shape", + "mf6dimension (noutlets)", "tagged true", "in_record true", "reader urword", @@ -1211,8 +1279,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name rough", - "type string", + "type double precision", "shape", + "mf6dimension (noutlets)", "tagged true", "in_record true", "reader urword", @@ -1221,8 +1290,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name width", - "type string", + "type double precision", "shape", + "mf6dimension (noutlets)", "tagged true", "in_record true", "reader urword", @@ -1231,8 +1301,9 @@ class ModflowGwflak(MFPackage): [ "block period", "name slope", - "type string", + "type double precision", "shape", + "mf6dimension (noutlets)", "tagged true", "in_record true", "reader urword", @@ -1241,6 +1312,7 @@ class ModflowGwflak(MFPackage): [ "block period", "name auxiliaryrecord", + "mf6internal auxrecord", "type record auxiliary auxname auxval", "shape", "tagged", @@ -1254,6 +1326,7 @@ class ModflowGwflak(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", @@ -1293,7 +1366,8 @@ def __init__( timeseries=None, observations=None, mover=None, - dev_force_fallback=None, + dev_force_legacy=None, + dev_force_legacy_lake=None, dev_groundwater_head_conductance=None, dev_maximum_outlet_depth=None, surfdep=None, @@ -1302,6 +1376,7 @@ def __init__( maximum_stage_change=None, time_conversion=None, length_conversion=None, + dev_no_final_check=None, nlakes=None, noutlets=None, ntables=None, @@ -1349,8 +1424,9 @@ def __init__( "obs", observations, "continuous", self._obs_filerecord ) self.mover = self.build_mfdata("mover", mover) - self.dev_force_fallback = self.build_mfdata( - "dev_force_fallback", dev_force_fallback + self.dev_force_legacy = self.build_mfdata("dev_force_legacy", dev_force_legacy) + self.dev_force_legacy_lake = self.build_mfdata( + "dev_force_legacy_lake", dev_force_legacy_lake ) self.dev_groundwater_head_conductance = self.build_mfdata( "dev_groundwater_head_conductance", dev_groundwater_head_conductance @@ -1370,6 +1446,9 @@ def __init__( self.length_conversion = self.build_mfdata( "length_conversion", length_conversion ) + self.dev_no_final_check = self.build_mfdata( + "dev_no_final_check", dev_no_final_check + ) self.nlakes = self.build_mfdata("nlakes", nlakes) self.noutlets = self.build_mfdata("noutlets", noutlets) self.ntables = self.build_mfdata("ntables", ntables) diff --git a/flopy/mf6/modflow/mfgwfmaw.py b/flopy/mf6/modflow/mfgwfmaw.py index 7d3aa08915..49e7473c59 100644 --- a/flopy/mf6/modflow/mfgwfmaw.py +++ b/flopy/mf6/modflow/mfgwfmaw.py @@ -473,6 +473,7 @@ class ModflowGwfmaw(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -487,6 +488,7 @@ class ModflowGwfmaw(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -494,6 +496,7 @@ class ModflowGwfmaw(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -528,6 +531,7 @@ class ModflowGwfmaw(MFPackage): [ "block options", "name budget_filerecord", + "mf6internal budfilerec", "type record budget fileout budgetfile", "shape", "reader urword", @@ -568,6 +572,7 @@ class ModflowGwfmaw(MFPackage): [ "block options", "name budgetcsv_filerecord", + "mf6internal budcsvfilerec", "type record budgetcsv fileout budgetcsvfile", "shape", "reader urword", @@ -615,6 +620,7 @@ class ModflowGwfmaw(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal inonvert", ], [ "block options", @@ -640,6 +646,7 @@ class ModflowGwfmaw(MFPackage): [ "block options", "name mfrcsv_filerecord", + "mf6internal mfrcsvfilerec", "type record maw_flow_reduce_csv fileout mfrcsvfile", "shape", "reader urword", @@ -649,6 +656,7 @@ class ModflowGwfmaw(MFPackage): [ "block options", "name maw_flow_reduce_csv", + "mf6internal maw_flw_red_csv", "type keyword", "shape", "in_record true", @@ -755,11 +763,13 @@ class ModflowGwfmaw(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ieffradopt", ], [ "block dimensions", "name nmawwells", "type integer", + "mf6internal maxbound", "reader urword", "optional false", ], @@ -779,6 +789,7 @@ class ModflowGwfmaw(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -805,6 +816,7 @@ class ModflowGwfmaw(MFPackage): "shape", "tagged false", "in_record true", + "time_series true", "reader urword", ], [ @@ -861,6 +873,7 @@ class ModflowGwfmaw(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal conndata_ifno", ], [ "block connectiondata", @@ -932,6 +945,7 @@ class ModflowGwfmaw(MFPackage): "tagged false", "in_record true", "reader urword", + "mf6internal angledata_ifno", "numeric_index true", ], [ @@ -942,6 +956,7 @@ class ModflowGwfmaw(MFPackage): "tagged false", "in_record true", "reader urword", + "mf6internal angledata_icon", "numeric_index true", ], [ @@ -963,6 +978,7 @@ class ModflowGwfmaw(MFPackage): "in_record true", "reader urword", "optional true", + "mf6internal connlen", ], [ "block period", @@ -1014,6 +1030,7 @@ class ModflowGwfmaw(MFPackage): [ "block period", "name flowing_wellrecord", + "mf6internal flwwellrecord", "type record flowing_well fwelev fwcond fwrlen", "shape", "tagged", @@ -1122,6 +1139,7 @@ class ModflowGwfmaw(MFPackage): [ "block period", "name rate_scalingrecord", + "mf6internal ratescalingrec", "type record rate_scaling pump_elevation scaling_length", "shape", "tagged", @@ -1170,6 +1188,7 @@ class ModflowGwfmaw(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwfmvr.py b/flopy/mf6/modflow/mfgwfmvr.py index 9955dfdc15..1ea1af7e2d 100644 --- a/flopy/mf6/modflow/mfgwfmvr.py +++ b/flopy/mf6/modflow/mfgwfmvr.py @@ -60,7 +60,7 @@ class ModflowGwfmvr(MFPackage): provided by the user, then packages are assigned a default value, which is the package acronym followed by a hyphen and the package number. For example, the first Drain Package is named DRN-1. The second Drain Package is named DRN-2, - and so forth. + and so forth. Each package may be listed only once in the PACKAGES block. perioddata : [(mname1, pname1, id1, mname2, pname2, id2, mvrtype, value)] * mname1 : string diff --git a/flopy/mf6/modflow/mfgwfsfr.py b/flopy/mf6/modflow/mfgwfsfr.py index 753582615b..e224c6b062 100644 --- a/flopy/mf6/modflow/mfgwfsfr.py +++ b/flopy/mf6/modflow/mfgwfsfr.py @@ -144,6 +144,9 @@ class ModflowGwfsfr(MFPackage): keyword that deactivates the checks of reach geometry relative to the model grid and of reach parameters for reasonable values. this development option is not supported. + dev_no_final_check : keyword + keyword that deactivates the final convergence check of the change in stream + flow routing stages and flows. this development option is not supported. nreaches : integer integer value specifying the number of stream reaches. there must be nreaches entries in the packagedata block. @@ -184,7 +187,7 @@ class ModflowGwfsfr(MFPackage): block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * man : string + * man : double precision real or character value that defines the Manning's roughness coefficient for the reach. MAN must be greater than zero. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be @@ -307,7 +310,7 @@ class ModflowGwfsfr(MFPackage): that can be used to start the sfrsetting string include: status, bedk, manning, stage, inflow, rainfall, evaporation, runoff, diversion, upstream_fraction, and auxiliary. - * stage : string + * stage : double precision real or character value that defines the stage for the reach. the specified stage is only applied if the reach uses the simple routing option. if stage is not specified for reaches that use the simple routing option, the specified @@ -328,33 +331,33 @@ class ModflowGwfsfr(MFPackage): and inactive downstream reach should be changed to 1.0 and 0.0, respectively, to ensure that the active reach receives all of the downstream outflow from the upstream reach. by default, status is active. - * bedk : string + * bedk : double precision real or character value that defines the hydraulic conductivity of the reach streambed. bedk can be any positive value if the reach is not connected to an underlying gwf cell. otherwise, bedk must be greater than zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * manning : string + * manning : double precision real or character value that defines the manning's roughness coefficient for the reach. manning must be greater than zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * inflow : string + * inflow : double precision real or character value that defines the volumetric inflow rate for the streamflow routing reach. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. by default, inflow rates are zero for each reach. - * rainfall : string + * rainfall : double precision real or character value that defines the volumetric rate per unit area of water added by precipitation directly on the streamflow routing reach. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. by default, rainfall rates are zero for each reach. - * evaporation : string + * evaporation : double precision real or character value that defines the volumetric rate per unit area of water subtracted by evaporation from the streamflow routing reach. a positive evaporation rate should be provided. if the options block includes a @@ -365,7 +368,7 @@ class ModflowGwfsfr(MFPackage): runoff but excluding groundwater leakage into the reach) the volumetric evaporation rate is limited to the sources of water to the reach. by default, evaporation rates are zero for each reach. - * runoff : string + * runoff : double precision real or character value that defines the volumetric rate of diffuse overland runoff that enters the streamflow routing reach. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be @@ -502,6 +505,7 @@ class ModflowGwfsfr(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -516,6 +520,7 @@ class ModflowGwfsfr(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -523,6 +528,7 @@ class ModflowGwfsfr(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -536,6 +542,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name stage", + "mf6internal options_stage", "type keyword", "shape", "in_record true", @@ -557,6 +564,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name budget_filerecord", + "mf6internal budfilerec", "type record budget fileout budgetfile", "shape", "reader urword", @@ -597,6 +605,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name budgetcsv_filerecord", + "mf6internal budcsvfilerec", "type record budgetcsv fileout budgetcsvfile", "shape", "reader urword", @@ -627,6 +636,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name package_convergence_filerecord", + "mf6internal pkgconvfilerec", "type record package_convergence fileout package_convergence_filename", "shape", "reader urword", @@ -636,6 +646,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name package_convergence", + "mf6internal pkgconv", "type keyword", "shape", "in_record true", @@ -646,6 +657,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name package_convergence_filename", + "mf6internal pkgconvfname", "type string", "shape", "in_record true", @@ -678,6 +690,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name filein", + "mf6internal options_filein", "type keyword", "shape", "in_record true", @@ -738,6 +751,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name maximum_picard_iterations", + "mf6internal maxsfrpicard", "type integer", "reader urword", "optional true", @@ -745,6 +759,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name maximum_iterations", + "mf6internal maxsfrit", "type integer", "reader urword", "optional true", @@ -752,6 +767,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name maximum_depth_change", + "mf6internal dmaxchg", "type double precision", "reader urword", "optional true", @@ -767,6 +783,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name length_conversion", + "mf6internal lengthconv", "type double precision", "reader urword", "optional true", @@ -781,6 +798,7 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name dev_storage_weight", + "mf6internal storage_weight", "type double precision", "reader urword", "optional true", @@ -788,6 +806,15 @@ class ModflowGwfsfr(MFPackage): [ "block options", "name dev_no_check", + "mf6internal icheck", + "type keyword", + "reader urword", + "optional true", + ], + [ + "block options", + "name dev_no_final_check", + "mf6internal iconvchk", "type keyword", "reader urword", "optional true", @@ -796,6 +823,7 @@ class ModflowGwfsfr(MFPackage): "block dimensions", "name nreaches", "type integer", + "mf6internal maxbound", "reader urword", "optional false", ], @@ -803,12 +831,13 @@ class ModflowGwfsfr(MFPackage): "block packagedata", "name packagedata", "type recarray ifno cellid rlen rwid rgrd rtp rbth rhk man ncon ustrf ndv aux boundname", - "shape (maxbound)", + "shape (nreaches)", "reader urword", ], [ "block packagedata", "name ifno", + "mf6internal packagedata_ifno", "type integer", "shape", "tagged false", @@ -883,7 +912,7 @@ class ModflowGwfsfr(MFPackage): [ "block packagedata", "name man", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -945,12 +974,13 @@ class ModflowGwfsfr(MFPackage): "type recarray ifno tab6 filein tab6_filename", "shape", "valid", - "optional false", + "optional true", "reader urword", ], [ "block crosssections", "name ifno", + "mf6internal crosssect_ifno", "type integer", "shape", "tagged false", @@ -961,6 +991,7 @@ class ModflowGwfsfr(MFPackage): [ "block crosssections", "name tab6", + "mf6internal xs_tab6", "type keyword", "shape", "in_record true", @@ -971,6 +1002,7 @@ class ModflowGwfsfr(MFPackage): [ "block crosssections", "name filein", + "mf6internal xs_filein", "type keyword", "shape", "in_record true", @@ -981,6 +1013,7 @@ class ModflowGwfsfr(MFPackage): [ "block crosssections", "name tab6_filename", + "mf6internal sfrtab6_filename", "type string", "preserve_case true", "in_record true", @@ -992,12 +1025,13 @@ class ModflowGwfsfr(MFPackage): "block connectiondata", "name connectiondata", "type recarray ifno ic", - "shape (maxbound)", + "shape (nreaches)", "reader urword", ], [ "block connectiondata", "name ifno", + "mf6internal conndata_ifno", "type integer", "shape", "tagged false", @@ -1021,12 +1055,14 @@ class ModflowGwfsfr(MFPackage): "block diversions", "name diversions", "type recarray ifno idv iconr cprior", - "shape (maxbound)", + "shape", + "optional true", "reader urword", ], [ "block diversions", "name ifno", + "mf6internal diversions_ifno", "type integer", "shape", "tagged false", @@ -1037,6 +1073,7 @@ class ModflowGwfsfr(MFPackage): [ "block diversions", "name idv", + "mf6internal diversions_idv", "type integer", "shape", "tagged false", @@ -1067,14 +1104,15 @@ class ModflowGwfsfr(MFPackage): "block initialstages", "name initialstages", "type recarray ifno initialstage", - "shape (maxbound)", + "shape (nreaches)", "valid", - "optional false", + "optional true", "reader urword", ], [ "block initialstages", "name ifno", + "mf6internal initstage_ifno", "type integer", "shape", "tagged false", @@ -1143,7 +1181,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name bedk", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1153,7 +1191,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name manning", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1163,7 +1201,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name stage", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1173,7 +1211,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name inflow", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1183,7 +1221,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1193,7 +1231,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1203,7 +1241,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -1250,15 +1288,18 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name upstream_fraction", + "mf6internal upstream_frac", "type double precision", "shape", "tagged true", "in_record true", "reader urword", + "time_series true", ], [ "block period", "name cross_sectionrecord", + "mf6internal xsectionrecord", "type record cross_section tab6 filein tab6_filename", "shape", "tagged", @@ -1317,6 +1358,7 @@ class ModflowGwfsfr(MFPackage): [ "block period", "name auxiliary", + "mf6internal period_auxiliary", "type keyword", "shape", "in_record true", @@ -1370,6 +1412,7 @@ def __init__( time_conversion=None, dev_storage_weight=None, dev_no_check=None, + dev_no_final_check=None, nreaches=None, packagedata=None, crosssections=None, @@ -1436,6 +1479,9 @@ def __init__( "dev_storage_weight", dev_storage_weight ) self.dev_no_check = self.build_mfdata("dev_no_check", dev_no_check) + self.dev_no_final_check = self.build_mfdata( + "dev_no_final_check", dev_no_final_check + ) self.nreaches = self.build_mfdata("nreaches", nreaches) self.packagedata = self.build_mfdata("packagedata", packagedata) self.crosssections = self.build_mfdata("crosssections", crosssections) diff --git a/flopy/mf6/modflow/mfgwfuzf.py b/flopy/mf6/modflow/mfgwfuzf.py index 89358247bd..9374ca6d99 100644 --- a/flopy/mf6/modflow/mfgwfuzf.py +++ b/flopy/mf6/modflow/mfgwfuzf.py @@ -178,12 +178,12 @@ class ModflowGwfuzf(MFPackage): * ifno : integer integer value that defines the feature (UZF object) number associated with the specified PERIOD data on the line. - * finf : string + * finf : double precision real or character value that defines the applied infiltration rate of the UZF cell (:math:`LT^{-1}`). If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * pet : string + * pet : double precision real or character value that defines the potential evapotranspiration rate of the UZF cell and specified GWF cell. Evapotranspiration is first removed from the unsaturated zone and any remaining potential evapotranspiration is applied @@ -193,7 +193,7 @@ class ModflowGwfuzf(MFPackage): specified in the OPTIONS block. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * extdp : string + * extdp : double precision real or character value that defines the evapotranspiration extinction depth of the UZF cell. If IVERTCON is greater than zero and EXTDP extends below the GWF cell bottom then remaining potential evapotranspiration is applied to the @@ -202,7 +202,7 @@ class ModflowGwfuzf(MFPackage): TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * extwc : string + * extwc : double precision real or character value that defines the evapotranspiration extinction water content of the UZF cell. EXTWC is always specified, but is only used if SIMULATE_ET and UNSAT_ETWC are specified in the OPTIONS block. The @@ -213,20 +213,20 @@ class ModflowGwfuzf(MFPackage): If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * ha : string + * ha : double precision real or character value that defines the air entry potential (head) of the UZF cell. HA is always specified, but is only used if SIMULATE_ET and UNSAT_ETAE are specified in the OPTIONS block. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * hroot : string + * hroot : double precision real or character value that defines the root potential (head) of the UZF cell. HROOT is always specified, but is only used if SIMULATE_ET and UNSAT_ETAE are specified in the OPTIONS block. If the Options block includes a TIMESERIESFILE entry (see the 'Time-Variable Input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * rootact : string + * rootact : double precision real or character value that defines the root activity function of the UZF cell. ROOTACT is the length of roots in a given volume of soil divided by that volume. Values range from 0 to about 3 :math:`cm^{-2}`, depending on the plant @@ -303,6 +303,7 @@ class ModflowGwfuzf(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -317,6 +318,7 @@ class ModflowGwfuzf(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -356,6 +358,7 @@ class ModflowGwfuzf(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -396,6 +399,7 @@ class ModflowGwfuzf(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -426,6 +430,7 @@ class ModflowGwfuzf(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal pkgconvfilerec", ], [ "block options", @@ -436,6 +441,7 @@ class ModflowGwfuzf(MFPackage): "reader urword", "tagged true", "optional false", + "mf6internal pkgconv", ], [ "block options", @@ -446,6 +452,7 @@ class ModflowGwfuzf(MFPackage): "reader urword", "tagged false", "optional false", + "mf6internal pkgconvfname", ], [ "block options", @@ -532,6 +539,7 @@ class ModflowGwfuzf(MFPackage): [ "block options", "name dev_no_final_check", + "mf6internal iconvchk", "type keyword", "reader urword", "optional true", @@ -591,6 +599,7 @@ class ModflowGwfuzf(MFPackage): "type integer", "reader urword", "optional false", + "mf6internal maxbound", ], [ "block dimensions", @@ -618,6 +627,7 @@ class ModflowGwfuzf(MFPackage): [ "block packagedata", "name ifno", + "mf6internal packagedata_ifno", "type integer", "shape", "tagged false", @@ -749,7 +759,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name finf", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -759,7 +769,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name pet", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -769,7 +779,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name extdp", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -779,7 +789,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name extwc", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -789,7 +799,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name ha", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -799,7 +809,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name hroot", - "type string", + "type double precision", "shape", "tagged false", "in_record true", @@ -809,7 +819,7 @@ class ModflowGwfuzf(MFPackage): [ "block period", "name rootact", - "type string", + "type double precision", "shape", "tagged false", "in_record true", diff --git a/flopy/mf6/modflow/mfgwtfmi.py b/flopy/mf6/modflow/mfgwtfmi.py index 15e4130f53..349eead71c 100644 --- a/flopy/mf6/modflow/mfgwtfmi.py +++ b/flopy/mf6/modflow/mfgwtfmi.py @@ -33,12 +33,12 @@ class ModflowGwtfmi(MFPackage): convergence tolerances. packagedata : [(flowtype, filein, fname)] * flowtype : string - is the word GWFBUDGET, GWFHEAD, GWFMOVER or the name of an advanced GWF stress - package. If GWFBUDGET is specified, then the corresponding file must be a - budget file. If GWFHEAD is specified, the file must be a head file. If GWFGRID - is specified, the file must be a binary grid file. If an advanced GWF stress - package name appears then the corresponding file must be the budget file saved - by a LAK, SFR, MAW or UZF Package. + is the word GWFBUDGET, GWFHEAD, GWFGRID, GWFMOVER or the name of an advanced + GWF stress package. If GWFBUDGET is specified, then the corresponding file + must be a budget file. If GWFHEAD is specified, the file must be a head file. + If GWFGRID is specified, the file must be a binary grid file. If an advanced + GWF stress package name appears then the corresponding file must be the budget + file saved by a LAK, SFR, MAW or UZF Package. * filein : keyword keyword to specify that an input filename is expected next. * fname : string diff --git a/flopy/mf6/modflow/mfgwtlkt.py b/flopy/mf6/modflow/mfgwtlkt.py index 38ee2e07c5..9bfb5faef7 100644 --- a/flopy/mf6/modflow/mfgwtlkt.py +++ b/flopy/mf6/modflow/mfgwtlkt.py @@ -125,7 +125,7 @@ class ModflowGwtlkt(MFPackage): the lake package supports a 'withdrawal' flow term. if this withdrawal term is active, then water will be withdrawn from the lake at the calculated concentration of the lake. - * concentration : string + * concentration : double precision real or character value that defines the concentration for the lake. the specified concentration is only applied if the lake is a constant concentration lake. if the options block includes a timeseriesfile entry (see the 'time- @@ -138,12 +138,12 @@ class ModflowGwtlkt(MFPackage): mass fluxes into or out of the lake and the inactive value will be written for the lake concentration. if a lake is constant, then the concentration for the lake will be fixed at the user specified value. - * rainfall : string + * rainfall : double precision real or character value that defines the rainfall solute concentration :math:`(ml^{-3})` for the lake. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * evaporation : string + * evaporation : double precision real or character value that defines the concentration of evaporated water :math:`(ml^{-3})` for the lake. if this concentration value is larger than the simulated concentration in the lake, then the evaporated water will be removed @@ -151,13 +151,13 @@ class ModflowGwtlkt(MFPackage): timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * runoff : string + * runoff : double precision real or character value that defines the concentration of runoff :math:`(ml^{-3})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * ext_inflow : string + * ext_inflow : double precision real or character value that defines the concentration of external inflow :math:`(ml^{-3})` for the lake. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable @@ -215,6 +215,7 @@ class ModflowGwtlkt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -231,6 +232,7 @@ class ModflowGwtlkt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -246,6 +248,7 @@ class ModflowGwtlkt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -253,6 +256,7 @@ class ModflowGwtlkt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -260,6 +264,7 @@ class ModflowGwtlkt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -267,6 +272,7 @@ class ModflowGwtlkt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -276,6 +282,7 @@ class ModflowGwtlkt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal conc_filerec", ], [ "block options", @@ -306,6 +313,7 @@ class ModflowGwtlkt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -346,6 +354,7 @@ class ModflowGwtlkt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -448,12 +457,13 @@ class ModflowGwtlkt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray ifno strt aux boundname", - "shape (maxbound)", + "shape (nlakes)", "reader urword", ], [ @@ -465,6 +475,7 @@ class ModflowGwtlkt(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -546,7 +557,7 @@ class ModflowGwtlkt(MFPackage): [ "block period", "name concentration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -556,7 +567,7 @@ class ModflowGwtlkt(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -566,7 +577,7 @@ class ModflowGwtlkt(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -576,7 +587,7 @@ class ModflowGwtlkt(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -586,7 +597,7 @@ class ModflowGwtlkt(MFPackage): [ "block period", "name ext-inflow", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -609,6 +620,7 @@ class ModflowGwtlkt(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwtmwt.py b/flopy/mf6/modflow/mfgwtmwt.py index 2cdcd4f487..91dee8056d 100644 --- a/flopy/mf6/modflow/mfgwtmwt.py +++ b/flopy/mf6/modflow/mfgwtmwt.py @@ -125,7 +125,7 @@ class ModflowGwtmwt(MFPackage): supports a 'withdrawal' flow term. if this withdrawal term is active, then water will be withdrawn from the well at the calculated concentration of the well. - * concentration : string + * concentration : double precision real or character value that defines the concentration for the well. the specified concentration is only applied if the well is a constant concentration well. if the options block includes a timeseriesfile entry (see the 'time- @@ -138,7 +138,7 @@ class ModflowGwtmwt(MFPackage): mass fluxes into or out of the well and the inactive value will be written for the well concentration. if a well is constant, then the concentration for the well will be fixed at the user specified value. - * rate : string + * rate : double precision real or character value that defines the injection solute concentration :math:`(ml^{-3})` for the well. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a @@ -195,6 +195,7 @@ class ModflowGwtmwt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -211,6 +212,7 @@ class ModflowGwtmwt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -226,6 +228,7 @@ class ModflowGwtmwt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -233,6 +236,7 @@ class ModflowGwtmwt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -240,6 +244,7 @@ class ModflowGwtmwt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -247,6 +252,7 @@ class ModflowGwtmwt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -256,6 +262,7 @@ class ModflowGwtmwt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal conc_filerec", ], [ "block options", @@ -286,6 +293,7 @@ class ModflowGwtmwt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -326,6 +334,7 @@ class ModflowGwtmwt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -428,12 +437,13 @@ class ModflowGwtmwt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray ifno strt aux boundname", - "shape (maxbound)", + "shape (nmawwells)", "reader urword", ], [ @@ -445,6 +455,7 @@ class ModflowGwtmwt(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -526,7 +537,7 @@ class ModflowGwtmwt(MFPackage): [ "block period", "name concentration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -536,7 +547,7 @@ class ModflowGwtmwt(MFPackage): [ "block period", "name rate", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -559,6 +570,7 @@ class ModflowGwtmwt(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwtsft.py b/flopy/mf6/modflow/mfgwtsft.py index 33b7d56511..4c28adf83c 100644 --- a/flopy/mf6/modflow/mfgwtsft.py +++ b/flopy/mf6/modflow/mfgwtsft.py @@ -125,7 +125,7 @@ class ModflowGwtsft(MFPackage): terms. concentrations cannot be specified for all flow terms. for example, the streamflow package supports a 'diversion' flow term. diversion water will be routed using the calculated concentration of the reach. - * concentration : string + * concentration : double precision real or character value that defines the concentration for the reach. the specified concentration is only applied if the reach is a constant concentration reach. if the options block includes a timeseriesfile entry (see @@ -138,12 +138,12 @@ class ModflowGwtsft(MFPackage): mass fluxes into or out of the reach and the inactive value will be written for the reach concentration. if a reach is constant, then the concentration for the reach will be fixed at the user specified value. - * rainfall : string + * rainfall : double precision real or character value that defines the rainfall solute concentration :math:`(ml^{-3})` for the reach. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * evaporation : string + * evaporation : double precision real or character value that defines the concentration of evaporated water :math:`(ml^{-3})` for the reach. if this concentration value is larger than the simulated concentration in the reach, then the evaporated water will be removed @@ -151,13 +151,13 @@ class ModflowGwtsft(MFPackage): timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * runoff : string + * runoff : double precision real or character value that defines the concentration of runoff :math:`(ml^{-3})` for the reach. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * inflow : string + * inflow : double precision real or character value that defines the concentration of inflow :math:`(ml^{-3})` for the reach. value must be greater than or equal to zero. if the options block includes a timeseriesfile entry (see the 'time-variable @@ -217,6 +217,7 @@ class ModflowGwtsft(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -233,6 +234,7 @@ class ModflowGwtsft(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -248,6 +250,7 @@ class ModflowGwtsft(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -255,6 +258,7 @@ class ModflowGwtsft(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -262,6 +266,7 @@ class ModflowGwtsft(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -269,6 +274,7 @@ class ModflowGwtsft(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -278,6 +284,7 @@ class ModflowGwtsft(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal conc_filerec", ], [ "block options", @@ -308,6 +315,7 @@ class ModflowGwtsft(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -348,6 +356,7 @@ class ModflowGwtsft(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -450,12 +459,13 @@ class ModflowGwtsft(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray ifno strt aux boundname", - "shape (maxbound)", + "shape (nreaches)", "reader urword", ], [ @@ -467,6 +477,7 @@ class ModflowGwtsft(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -548,7 +559,7 @@ class ModflowGwtsft(MFPackage): [ "block period", "name concentration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -558,7 +569,7 @@ class ModflowGwtsft(MFPackage): [ "block period", "name rainfall", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -568,7 +579,7 @@ class ModflowGwtsft(MFPackage): [ "block period", "name evaporation", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -578,7 +589,7 @@ class ModflowGwtsft(MFPackage): [ "block period", "name runoff", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -588,7 +599,7 @@ class ModflowGwtsft(MFPackage): [ "block period", "name inflow", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -611,6 +622,7 @@ class ModflowGwtsft(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfgwtuzt.py b/flopy/mf6/modflow/mfgwtuzt.py index f4be064b9e..dd7e0b15ec 100644 --- a/flopy/mf6/modflow/mfgwtuzt.py +++ b/flopy/mf6/modflow/mfgwtuzt.py @@ -125,7 +125,7 @@ class ModflowGwtuzt(MFPackage): infiltration, uzet, and auxiliary. these settings are used to assign the concentration of associated with the corresponding flow terms. concentrations cannot be specified for all flow terms. - * concentration : string + * concentration : double precision real or character value that defines the concentration for the unsaturated zone flow cell. the specified concentration is only applied if the unsaturated zone flow cell is a constant concentration cell. if the options block includes a @@ -139,13 +139,13 @@ class ModflowGwtuzt(MFPackage): solute mass fluxes into or out of the uzf cell and the inactive value will be written for the uzf cell concentration. if a uzf cell is constant, then the concentration for the uzf cell will be fixed at the user specified value. - * infiltration : string + * infiltration : double precision real or character value that defines the infiltration solute concentration :math:`(ml^{-3})` for the uzf cell. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * uzet : string + * uzet : double precision real or character value that defines the concentration of unsaturated zone evapotranspiration water :math:`(ml^{-3})` for the uzf cell. if this concentration value is larger than the simulated concentration in the uzf cell, @@ -205,6 +205,7 @@ class ModflowGwtuzt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal flow_pkg_name", ], [ "block options", @@ -221,6 +222,7 @@ class ModflowGwtuzt(MFPackage): "shape", "reader urword", "optional true", + "mf6internal fp_aux_name", ], [ "block options", @@ -236,6 +238,7 @@ class ModflowGwtuzt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprpak", ], [ "block options", @@ -243,6 +246,7 @@ class ModflowGwtuzt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprconc", ], [ "block options", @@ -250,6 +254,7 @@ class ModflowGwtuzt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal iprflow", ], [ "block options", @@ -257,6 +262,7 @@ class ModflowGwtuzt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal ipakcb", ], [ "block options", @@ -266,6 +272,7 @@ class ModflowGwtuzt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal conc_filerec", ], [ "block options", @@ -296,6 +303,7 @@ class ModflowGwtuzt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budfilerec", ], [ "block options", @@ -336,6 +344,7 @@ class ModflowGwtuzt(MFPackage): "reader urword", "tagged true", "optional true", + "mf6internal budcsvfilerec", ], [ "block options", @@ -438,12 +447,13 @@ class ModflowGwtuzt(MFPackage): "type keyword", "reader urword", "optional true", + "mf6internal dev_nonexpanding", ], [ "block packagedata", "name packagedata", "type recarray ifno strt aux boundname", - "shape (maxbound)", + "shape (nuzfcells)", "reader urword", ], [ @@ -455,6 +465,7 @@ class ModflowGwtuzt(MFPackage): "in_record true", "reader urword", "numeric_index true", + "mf6internal packagedata_ifno", ], [ "block packagedata", @@ -536,7 +547,7 @@ class ModflowGwtuzt(MFPackage): [ "block period", "name concentration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -546,7 +557,7 @@ class ModflowGwtuzt(MFPackage): [ "block period", "name infiltration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -556,7 +567,7 @@ class ModflowGwtuzt(MFPackage): [ "block period", "name uzet", - "type string", + "type double precision", "shape", "tagged true", "in_record true", @@ -579,6 +590,7 @@ class ModflowGwtuzt(MFPackage): "shape", "in_record true", "reader urword", + "mf6internal period_auxiliary", ], [ "block period", diff --git a/flopy/mf6/modflow/mfmvr.py b/flopy/mf6/modflow/mfmvr.py index adb12a30c9..0e68ed857a 100644 --- a/flopy/mf6/modflow/mfmvr.py +++ b/flopy/mf6/modflow/mfmvr.py @@ -52,7 +52,7 @@ class ModflowMvr(MFPackage): provided by the user, then packages are assigned a default value, which is the package acronym followed by a hyphen and the package number. For example, the first Drain Package is named DRN-1. The second Drain Package is named DRN-2, - and so forth. + and so forth. Each package may be listed only once in the PACKAGES block. perioddata : [(mname1, pname1, id1, mname2, pname2, id2, mvrtype, value)] * mname1 : string diff --git a/flopy/mf6/modflow/mfnam.py b/flopy/mf6/modflow/mfnam.py index debd1e89ff..77e839de24 100644 --- a/flopy/mf6/modflow/mfnam.py +++ b/flopy/mf6/modflow/mfnam.py @@ -201,7 +201,6 @@ class ModflowNam(MFPackage): "preserve_case true", "type string", "reader urword", - "optional", ], [ "block models", diff --git a/flopy/mf6/modflow/mfprtprp.py b/flopy/mf6/modflow/mfprtprp.py index cd8bd64553..cb43fa1f17 100644 --- a/flopy/mf6/modflow/mfprtprp.py +++ b/flopy/mf6/modflow/mfprtprp.py @@ -181,10 +181,14 @@ class ModflowPrtprp(MFPackage): releasetimes block or configured via release_time_frequency, only one particle is released at that time. coincidence is evaluated up to the tolerance specified in release_time_tolerance, or :math:`epsilon times 10^{11}` by - default, where :math:`epsilon` is machine precision. if no release times are - configured via this setting, the releasetimes block, or the - release_time_frequency option, a single release time is configured at the - beginning of the first time step of the simulation's first stress period. + default, where :math:`epsilon` is machine precision. when such a coincidence + occurs, the release is reported under the earlier of the two time steps, i.e. + whichever time step the coincident release time falls in per the usual + convention that a time exactly on a time step boundary belongs to the earlier + time step, not the later one. if no release times are configured via this + setting, the releasetimes block, or the release_time_frequency option, a single + release time is configured at the beginning of the first time step of the + simulation's first stress period. * all : keyword keyword to indicate release at the start of all time steps in the period. * first : keyword diff --git a/flopy/mf6/modflow/mfutlncf.py b/flopy/mf6/modflow/mfutlncf.py index b3689207f6..da2f7fa665 100644 --- a/flopy/mf6/modflow/mfutlncf.py +++ b/flopy/mf6/modflow/mfutlncf.py @@ -21,7 +21,12 @@ class ModflowUtlncf(MFPackage): processing purposes only. wkt : [string] is the coordinate reference system (crs) well-known text string in wkt1 format - (ogc 01-009). written as the 'wkt' attribute on the projection variable. + (ogc 01-009). written as the 'wkt' attribute on the projection variable. for + supported projected crs types, individual cf-standard numeric grid_mapping + parameters are also written to the projection variable, for tools that read + these directly rather than parsing 'wkt'/'crs_wkt' text. for a + netcdf_structured export with a rotated grid (angrot not 0.0), wkt1 cannot + represent the rotation; provide 'crs_wkt' instead for correct grid positioning. ignored when latitude and longitude griddata arrays have been provided. to also embed the crs in the binary grid file, provide the same string as the 'crs' option in the discretization input; the discretization 'crs' option @@ -29,13 +34,23 @@ class ModflowUtlncf(MFPackage): (e.g. 'epsg:26918'). crs_wkt : [string] is the coordinate reference system (crs) well-known text string in wkt2 format - (iso 19162:2019), as required by cf-1.11. written as the 'crs_wkt' attribute on + (iso 19162:2019), as required by cf-1.13. written as the 'crs_wkt' attribute on the projection variable. when omitted, the 'wkt' value is written to 'crs_wkt' as a compatibility alias. ignored when latitude and longitude griddata arrays have been provided. to also embed the crs in the binary grid file, provide the same string as the 'crs' option in the discretization input; the discretization 'crs' option additionally accepts other pyproj-compatible identifiers such as - an epsg code (e.g. 'epsg:26918'). + an epsg code (e.g. 'epsg:26918'). for a netcdf_structured export with a rotated + grid (angrot not 0.0), a valid wkt2 projcrs string is required here to + correctly position the grid: it is wrapped in a derived crs that encodes the + rotation, resolved by any crs-aware consumer through ordinary reprojection; + without it the 'x'/'y' coordinate values remain grid-local rather than real- + world. for grid mappings transverse_mercator, lambert_conformal_conic (two- + standard-parallel variant only), and albers_conical_equal_area, individual cf- + standard numeric grid_mapping parameters (e.g. 'false_easting', + 'scale_factor_at_central_meridian') are also written to the projection + variable, for tools -- such as arcgis's classic netcdf connector -- that read + these directly rather than parsing 'wkt'/'crs_wkt' text. deflate : integer is the variable deflate level (0=min, 9=max) in the netcdf file. defining this parameter activates per-variable compression at the level specified. use with @@ -52,9 +67,8 @@ class ModflowUtlncf(MFPackage): parameter (netcdf_mesh2d) or the chunk_z, chunk_y, and chunk_x parameter set (netcdf_structured) to have an effect. chunk_face : integer - is the keyword used to provide a data chunk size for the face dimension in a - netcdf_mesh2d output file. must be used in combination with the the chunk_time - parameter to have an effect. + is the chunk size for the face dimension in a netcdf_mesh2d output file. must + be used in combination with the chunk_time parameter to have an effect. chunk_z : integer is the chunk size for the layer (vertical) dimension in a netcdf_structured output file. must be used in combination with the chunk_time, chunk_x, and @@ -74,11 +88,14 @@ class ModflowUtlncf(MFPackage): ncpl : integer is the number of cells in a projected plane layer. latitude : [double precision] - cell center latitude array. only supported for netcdf_structured export type. - providing latitude and longitude arrays enables geographic-coordinate output as - a special case for external tooling that requires it; the preferred approach is - to use the 'wkt' and 'crs_wkt' options instead. when latitude and longitude - arrays are provided, any 'wkt' or 'crs_wkt' option values are ignored. + cell center latitude array. only supported for netcdf_structured export type; + ignored, with a warning, for netcdf_mesh2d. providing latitude and longitude + arrays enables geographic-coordinate output for external tooling that requires + it, and is the more reliable option for a grid with non-uniform row or column + spacing, or a rotated grid, since some tools cannot otherwise correctly + position such a grid from 'wkt'/'crs_wkt' alone. for other grids, + 'wkt'/'crs_wkt' is the simpler choice. when latitude and longitude arrays are + provided, any 'wkt' or 'crs_wkt' option values are ignored. longitude : [double precision] cell center longitude array. only supported for netcdf_structured export type. see the latitude description for usage notes. diff --git a/flopy/mf6/modflow/mfutlspc.py b/flopy/mf6/modflow/mfutlspc.py index 1299b52e9a..cf7dfa9393 100644 --- a/flopy/mf6/modflow/mfutlspc.py +++ b/flopy/mf6/modflow/mfutlspc.py @@ -39,12 +39,12 @@ class ModflowUtlspc(MFPackage): line of information that is parsed into a keyword and values. keyword values that can be used to start the spcsetting string include: concentration and temperature. - * concentration : string + * concentration : double precision is the boundary concentration. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. by default, the concentration for each boundary feature is zero. - * temperature : string + * temperature : double precision is the user-supplied boundary temperature. if the options block includes a timeseriesfile entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a @@ -165,24 +165,22 @@ class ModflowUtlspc(MFPackage): [ "block period", "name concentration", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal concentration_in", ], [ "block period", "name temperature", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal temperature_in", ], ] diff --git a/flopy/mf6/modflow/mfutltas.py b/flopy/mf6/modflow/mfutltas.py index 1cd78d7559..22e5f0b761 100644 --- a/flopy/mf6/modflow/mfutltas.py +++ b/flopy/mf6/modflow/mfutltas.py @@ -118,7 +118,7 @@ class ModflowUtltas(MFPackage): "block attributes", "name interpolation_method", "type string", - "valid stepwise linear linearend", + "valid stepwise linear", "shape", "tagged false", "reader urword", diff --git a/flopy/mf6/modflow/mfutlts.py b/flopy/mf6/modflow/mfutlts.py index 6cc63f5ff5..8323df251b 100644 --- a/flopy/mf6/modflow/mfutlts.py +++ b/flopy/mf6/modflow/mfutlts.py @@ -22,7 +22,8 @@ class ModflowUtlts(MFPackage): time_series_namerecord : (names, time_series_names) xxx * names : keyword - xxx + keyword to specify the names of the time series. The 'NAME' keyword can be used + as a substitute for 'NAMES'. * time_series_names : [string] Name by which a package references a particular time-array series. The name must be unique among all time-array series used in a package. @@ -34,12 +35,13 @@ class ModflowUtlts(MFPackage): * interpolation_method : [string] Interpolation method, which is either STEPWISE or LINEAR. - interpolation_methodrecord_single : record + interpolation_methodrecord_single : (method, interpolation_method_single) xxx - method : keyword - xxx - interpolation_method_single : string - interpolation method, which is either stepwise or linear. + * method : keyword + xxx + * interpolation_method_single : string + Interpolation method, which is either STEPWISE or LINEAR. + sfacrecord : (sfacs, sfacval) xxx * sfacs : keyword @@ -48,14 +50,14 @@ class ModflowUtlts(MFPackage): Scale factor, which will multiply all array values in time series. SFAC is an optional attribute; if omitted, SFAC = 1.0. - sfacrecord_single : (sfacval) + sfacrecord_single : (sfac, sfacval) xxx + * sfac : keyword + xxx * sfacval : [double precision] Scale factor, which will multiply all array values in time series. SFAC is an optional attribute; if omitted, SFAC = 1.0. - sfac : keyword - xxx timeseries : [(ts_time, ts_array)] xxx * ts_time : double precision @@ -166,6 +168,7 @@ class ModflowUtlts(MFPackage): "shape", "reader urword", "optional false", + "in_record true", ], [ "block attributes", @@ -176,6 +179,7 @@ class ModflowUtlts(MFPackage): "tagged false", "reader urword", "optional false", + "in_record true", ], [ "block attributes", @@ -219,9 +223,9 @@ class ModflowUtlts(MFPackage): "name sfac", "type keyword", "shape", - "tagged false", "reader urword", "optional false", + "in_record true", ], [ "block timeseries", diff --git a/flopy/mf6/modflow/mfutltvk.py b/flopy/mf6/modflow/mfutltvk.py index 02df213e87..ba7ed0a138 100644 --- a/flopy/mf6/modflow/mfutltvk.py +++ b/flopy/mf6/modflow/mfutltvk.py @@ -38,13 +38,13 @@ class ModflowUtltvk(MFPackage): line of information that is parsed into a property name keyword and values. property name keywords that can be used to start the tvksetting string include: k, k22, and k33. - * k : string + * k : double precision is the new value to be assigned as the cell's hydraulic conductivity from the start of the specified stress period, as per k in the npf package. if the options block includes a ts6 entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * k22 : string + * k22 : double precision is the new value to be assigned as the cell's hydraulic conductivity of the second ellipsoid axis (or the ratio of k22/k if the k22overk npf package option is specified) from the start of the specified stress period, as per k22 in the @@ -52,7 +52,7 @@ class ModflowUtltvk(MFPackage): direction. if the options block includes a ts6 entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * k33 : string + * k33 : double precision is the new value to be assigned as the cell's hydraulic conductivity of the third ellipsoid axis (or the ratio of k33/k if the k33overk npf package option is specified) from the start of the specified stress period, as per k33 in the @@ -168,35 +168,32 @@ class ModflowUtltvk(MFPackage): [ "block period", "name k", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal k_in", ], [ "block period", "name k22", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal k22_in", ], [ "block period", "name k33", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal k33_in", ], ] diff --git a/flopy/mf6/modflow/mfutltvs.py b/flopy/mf6/modflow/mfutltvs.py index 2a167c4249..c852412de3 100644 --- a/flopy/mf6/modflow/mfutltvs.py +++ b/flopy/mf6/modflow/mfutltvs.py @@ -44,7 +44,7 @@ class ModflowUtltvs(MFPackage): line of information that is parsed into a property name keyword and values. property name keywords that can be used to start the tvssetting string include: ss and sy. - * ss : string + * ss : double precision is the new value to be assigned as the cell's specific storage (or storage coefficient if the storagecoefficient sto package option is specified) from the start of the specified stress period, as per ss in the sto package. specific @@ -52,7 +52,7 @@ class ModflowUtltvs(MFPackage): includes a ts6 entry (see the 'time-variable input' section), values can be obtained from a time series by entering the time-series name in place of a numeric value. - * sy : string + * sy : double precision is the new value to be assigned as the cell's specific yield from the start of the specified stress period, as per sy in the sto package. specific yield values must be greater than or equal to 0. if the options block includes a ts6 @@ -174,24 +174,22 @@ class ModflowUtltvs(MFPackage): [ "block period", "name ss", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal ss_in", ], [ "block period", "name sy", - "type string", + "type double precision", "shape", "tagged true", "in_record true", "reader urword", "time_series true", - "mf6internal sy_in", ], ]