aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rwxr-xr-xfullswof-utils/makeBoundary14
1 files changed, 3 insertions, 11 deletions
diff --git a/fullswof-utils/makeBoundary b/fullswof-utils/makeBoundary
index 512a968..87abf44 100755
--- a/fullswof-utils/makeBoundary
+++ b/fullswof-utils/makeBoundary
@@ -221,9 +221,9 @@ xregion = []
zregion = []
zmin = []
for p in range(num_panels):
- # identify regions:
- xregion.append(xin[marker_ind[p]:marker_ind[p+1]])
- zregion.append(zin[marker_ind[p]:marker_ind[p+1]])
+ # identify regions (note regions overlap by one cell):
+ xregion.append(xin[marker_ind[p]:marker_ind[p+1]+1])
+ zregion.append(zin[marker_ind[p]:marker_ind[p+1]+1])
# identify minimum heights within each panel:
zmin.append(zregion[p].min())
@@ -255,14 +255,6 @@ h_i = [[] for _ in range(num_panels)]
K_i = [[] for _ in range(num_panels)]
Q_i = [[] for _ in range(num_panels)]
for p in range(num_panels):
- if p == panel[0]-1 and zregion[p][-1] < zmax:
- # ensure end node in region to the left of channel is dry:
- xregion[p] = np.append(xregion[p], xin[marker_ind[p]])
- zregion[p] = np.append(zregion[p], zin[marker_ind[p]])
- if p == panel[0]+1 and zregion[p][0] < zmax:
- # ensure start node in region to the right of channel is dry:
- xregion[p] = np.insert(xregion[p], 0, xin[marker_ind[p]-1])
- zregion[p] = np.insert(zregion[p], 0, zin[marker_ind[p]-1])
if zmax > zmin[p]:
p_i[p], A_i[p], r_h[p], h_i[p], K_i[p], Q_i[p] = conveyance(
numH,