apply maxbound to nodelist

Author: mjreno

autotest/test_gwf_disv_uzf.py

@@ -173,7 +173,9 @@ def get_model(ws, name, array_input=False):
 
     # general-head boundary
     if array_input:
-        ghb = flopy.mf6.ModflowGwfghba(gwf, print_flows=True, bhead=abhead, cond=acond)
+        ghb = flopy.mf6.ModflowGwfghba(
+            gwf, print_flows=True, maxbound=20, bhead=abhead, cond=acond
+        )
     else:
         ghb = flopy.mf6.ModflowGwfghb(gwf, print_flows=True, stress_period_data=ghb_spd)
 

autotest/test_gwf_sfr_inactive02.py

@@ -71,7 +71,7 @@ def get_model(ws, name, array_input=False):
         cond = np.full(nlay * nrow * ncol, DNODATA, dtype=float)
         bhead[0] = 1.0
         cond[0] = 1e6
-        flopy.mf6.ModflowGwfghba(gwf, bhead=bhead, cond=cond)
+        flopy.mf6.ModflowGwfghba(gwf, maxbound=1, bhead=bhead, cond=cond)
     else:
         flopy.mf6.ModflowGwfghb(gwf, stress_period_data=[((0, 0, 0), 1.0, 1e6)])
 

autotest/test_gwf_uzf01.py

@@ -109,6 +109,7 @@ def get_model(ws, name, array_input=False):
             gwf,
             print_input=True,
             print_flows=True,
+            maxbound=1,
             bhead=bhead,
             cond=cond,
             save_flows=False,

autotest/test_gwf_vsc01.py

@@ -145,6 +145,7 @@ def get_model(idx, ws, array_input=False):
             temp[0][0, i, ncol - 1] = initial_temperature
         flopy.mf6.ModflowGwfghba(
             gwf,
+            maxbound=nrow,
             pname="GHB-1",
             auxiliary="temperature",
             bhead=bhead,

autotest/test_gwt_henry_nr.py

@@ -226,6 +226,7 @@ def get_model(ws, name, array_input=False):
             print_input=True,
             print_flows=True,
             save_flows=False,
+            maxbound=20,
             pname="GHB-1",
             auxiliary=["CONCENTRATION", "DENSITY"],
             bhead=bheadspd,

autotest/test_netcdf_gwf_disv_uzf.py

@@ -88,6 +88,9 @@ def check_output(idx, test):
         f.write("  PRINT_INPUT\n")
         f.write("  PRINT_FLOWS\n")
         f.write("END options\n\n")
+        f.write("BEGIN dimensions\n")
+        f.write("  MAXBOUND  20\n")
+        f.write("END dimensions\n\n")
         f.write("BEGIN period 1\n")
         f.write("  bhead NETCDF\n")
         f.write("  cond NETCDF\n")

autotest/test_netcdf_gwf_uzf01.py

@@ -96,6 +96,9 @@ def check_output(idx, test, export):
         f.write("  PRINT_FLOWS\n")
         f.write("  OBS6  FILEIN  gwf_uzf01a.ghb.obs\n")
         f.write("END options\n\n")
+        f.write("BEGIN dimensions\n")
+        f.write("  MAXBOUND  1\n")
+        f.write("END dimensions\n\n")
         f.write("BEGIN period 1\n")
         f.write("  bhead NETCDF\n")
         f.write("  cond NETCDF\n")

autotest/test_netcdf_gwf_vsc01.py

@@ -94,6 +94,9 @@ def check_output(idx, test, export):
         f.write("BEGIN options\n")
         f.write("  auxiliary  TEMPERATURE\n")
         f.write("END options\n\n")
+        f.write("BEGIN dimensions\n")
+        f.write("  MAXBOUND  10\n")
+        f.write("END dimensions\n\n")
         f.write("BEGIN period 1\n")
         f.write("  bhead NETCDF\n")
         f.write("  cond NETCDF\n")

autotest/test_netcdf_gwt_henry_nr.py

@@ -101,6 +101,9 @@ def check_output(idx, test, export):
         f.write("  PRINT_INPUT\n")
         f.write("  PRINT_FLOWS\n")
         f.write("END options\n\n")
+        f.write("BEGIN dimensions\n")
+        f.write("  MAXBOUND  20\n")
+        f.write("END dimensions\n\n")
         for i in range(1001):
             f.write(f"BEGIN period {i + 1}\n")
             f.write("  bhead NETCDF\n")

doc/mf6io/mf6ivar/dfn/gwf-ghba.dfn

@@ -109,6 +109,16 @@ longname export array variables to netcdf output files.
 description keyword that specifies input griddata arrays should be written to the model output netcdf file.
 extended true
 
+# --------------------- gwf ghba dimensions ---------------------
+
+block dimensions
+name maxbound
+type integer
+reader urword
+optional false
+longname maximum number of general-head boundaries in any stress period
+description REPLACE maxbound {'{#1}': 'general-head boundary'}
+
 # --------------------- gwf ghba period ---------------------
 
 block period