diff options
| -rwxr-xr-x | fullswof-utils/makeBoundary | 14 |
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, |