fix(get_disu_kwargs): incorrect indexing of delr and delc

modflowpy · Nov 20, 2023 · 7a73645 · 7a73645
1 parent 7cc5e6f
commit 7a73645
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 3 deletions.
diff --git a/autotest/test_gridutil.py b/autotest/test_gridutil.py
@@ -82,6 +82,15 @@ def test_get_lni_infers_layer_count_when_int_ncpl(ncpl, nodes, expected):
             np.array([-10]),
             np.array([-30.0, -50.0]),
         ),
+        (
+            1,  # nlay
+            3,  # nrow
+            4,  # ncol
+            np.array(4 * [4.]),  # delr
+            np.array(3 * [3.]),  # delc
+            np.array([-10]),  # top
+            np.array([-30.0]),  # botm
+        ),
     ],
 )
 def test_get_disu_kwargs(nlay, nrow, ncol, delr, delc, tp, botm):

diff --git a/flopy/utils/gridutil.py b/flopy/utils/gridutil.py
@@ -110,7 +110,7 @@ def get_nn(k, i, j):
                 n = get_nn(k, i, j)
                 ja.append(n)
                 iac[n] += 1
-                area[n] = delr[i] * delc[j]
+                area[n] = delr[j] * delc[i]
                 ihc.append(n + 1)
                 cl12.append(n + 1)
                 hwva.append(n + 1)
@@ -126,7 +126,7 @@ def get_nn(k, i, j):
                     ihc.append(0)
                     dz = botm[k - 1] - botm[k]
                     cl12.append(0.5 * dz)
-                    hwva.append(delr[i] * delc[j])
+                    hwva.append(delr[j] * delc[i])
                 # back
                 if i > 0:
                     ja.append(get_nn(k, i - 1, j))
@@ -165,7 +165,7 @@ def get_nn(k, i, j):
                     else:
                         dz = botm[k - 1] - botm[k]
                     cl12.append(0.5 * dz)
-                    hwva.append(delr[i] * delc[j])
+                    hwva.append(delr[j] * delc[i])
     ja = np.array(ja, dtype=int)
     nja = ja.shape[0]
     hwva = np.array(hwva, dtype=float)