Improve defaulted aquifer start points

sintefmath · Sep 19, 2024 · 5ac9d4b · 5ac9d4b
1 parent 3cdfcac
commit 5ac9d4b
Showing 1 changed file with 7 additions and 20 deletions.
diff --git a/src/CornerPointGrid/nnc_and_aquifers.jl b/src/CornerPointGrid/nnc_and_aquifers.jl
@@ -115,26 +115,13 @@ function mesh_add_numerical_aquifers!(mesh, AQUNUM, AQUCON)
     new_faces_neighbors = Tuple{Int, Int}[]
     for (i, aqucon) in enumerate(AQUCON)
         id, I_start, I_stop, J_start, J_stop, K_start, K_stop, dir, tranmult, opt, = aqucon
-        # Defaulted means start at 1
-        if I_start < 1
-            I_start = 1
-        end
-        if J_start < 1
-            J_start = 1
-        end
-        if K_start < 1
-            K_start = 1
-        end
-        # Defaulted means no upper bound
-        if I_stop < 1
-            I_stop = dims[1]
-        end
-        if J_stop < 1
-            J_stop = dims[2]
-        end
-        if K_stop < 1
-            K_stop = dims[3]
-        end
+        # Defaulted means start and stop at 1?
+        I_start = max(I_start, 1)
+        J_start = max(J_start, 1)
+        K_start = max(K_start, 1)
+        I_stop = max(I_stop, 1)
+        J_stop = max(J_stop, 1)
+        K_stop = max(K_stop, 1)
         bfaces = find_faces_for_aquifer(mesh, I_start:I_stop, J_start:J_stop, K_start:K_stop, dir, IJK)
         prm = aquifer_parameters[id]