flopy3 adapters (#127)

Author: mjreno

.github/workflows/ci.yml

@@ -63,7 +63,7 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        os: [ ubuntu-22.04, macos-13, windows-2022 ]
+        os: [ ubuntu-latest, macos-latest, windows-latest ]
         python: [ "311", "312", "313" ]
     env:
       GCC_V: 11

flopy4/mf6/gwf/dis.py

@@ -53,19 +53,19 @@ class Dis(Package):
     top: NDArray[np.floating] = array(
         block="griddata",
         default=1.0,
-        dims=("ncol", "nrow"),
+        dims=("nrow", "ncol"),
         converter=Converter(convert_array, takes_self=True, takes_field=True),
     )
     botm: NDArray[np.floating] = array(
         block="griddata",
         default=0.0,
-        dims=("ncol", "nrow", "nlay"),
+        dims=("nlay", "nrow", "ncol"),
         converter=Converter(convert_array, takes_self=True, takes_field=True),
     )
     idomain: NDArray[np.integer] = array(
         block="griddata",
         default=1,
-        dims=("ncol", "nrow", "nlay"),
+        dims=("nlay", "nrow", "ncol"),
         converter=Converter(convert_array, takes_self=True, takes_field=True),
     )
     nnodes: int = dim(

flopy4/mf6/ims.py

@@ -29,6 +29,7 @@ class Ims(Solution):
     inner_maximum: Optional[int] = field(block="linear", default=None)
     inner_dvclose: Optional[float] = field(block="linear", default=None)
     inner_rclose: Optional[float] = field(block="linear", default=None)
+    inner_hclose: Optional[float] = field(block="linear", default=None)
     rclose_option: Optional[str] = field(block="linear", default=None)
     linear_acceleration: Optional[str] = field(block="linear", default=None)
     relaxation_factor: Optional[float] = field(block="linear", default=None)