use lists of lists to store conveyance data.

1 files changed, 14 insertions, 2 deletions

diff --git a/makeBoundary.py b/makeBoundary.py
index f874f05..5cc4e2d 100755
--- a/makeBoundary.py
+++ b/makeBoundary.py
@@ -191,8 +191,20 @@ def plot_region(xdata, labelx,
     plt.plot(xdata, ydata3)
     plt.show()
 
-p_i, A_i, r_h, h_i, K_i, Q_i = conveyance(
-    numH, n_co_chan, xregion, zregion,zmin, zmax)
+p_i = [[] for _ in range(num_panels)]
+A_i = [[] for _ in range(num_panels)]
+r_h = [[] for _ in range(num_panels)]
+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):
+    p_i[p], A_i[p], r_h[p], h_i[p], K_i[p], Q_i[p] = conveyance(
+        numH,
+        n_co_chan,
+        xregion[p],
+        zregion[p],
+        zmin[p],
+        zmax)
 
 plot_region(h_i-zmin, 'maximum depth / m',
             r_h, 'hydraulic radius / m',