Merge branch 'fix_make_reservoir_axis_error' into 'main'

Fix axis error in make reservoir

See merge request !16

digger/make/reservoir.py

@@ -436,8 +436,8 @@ def reservoir(
         # determine the x and y values of the extent of eta for setting the regions.
 
         h_dam = dam_eta - z
-        hdam_x = np.max(h_dam, axis=0)
-        hdam_y = np.max(h_dam, axis=1)
+        hdam_x = np.max(h_dam, axis=1)
+        hdam_y = np.max(h_dam, axis=0)
 
         hdam_xinds = np.nonzero(hdam_x > 0)
         row_max = hdam_xinds[0][-1] + 1

tests/make/test_make_reservoir.py

@@ -50,5 +50,5 @@ def test_make_reservoir(tmp_path, datadir):
     # ensure dam is in the correct location.
     np.testing.assert_almost_equal(
         regions["1"],
-        [443198.9327530123, 443318.9327530123, 3909173.820103951, 3909043.820103951],
+        [443958.932753, 444088.932753, 3909933.820104, 3909813.820104],
     )