From c9c002df8c28a37d47bb665f18e34060701e5fd3 Mon Sep 17 00:00:00 2001 From: Paul Garlick Date: Mon, 14 Oct 2019 22:29:19 +0100 Subject: use list to store Manning's n coefficients. --- makeBoundary.py | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) (limited to 'makeBoundary.py') diff --git a/makeBoundary.py b/makeBoundary.py index 5cc4e2d..f20ed2e 100755 --- a/makeBoundary.py +++ b/makeBoundary.py @@ -23,10 +23,9 @@ definition_dict = read_definition('boundaryDefinition.txt') slope = abs(definition_dict["slope"]) # slope at top boundary target_flow = definition_dict["target_flow"] # imposed discharge location = definition_dict["location"] # boundary location -n_co_chan = definition_dict["n_co_chan"] # coefficient for inland water -n_co_west = definition_dict["n_co_west"] # coefficient for general surface -n_co_east = definition_dict["n_co_east"] # coefficient for general surface -# TODO: use weighted mean 'n' value. See http://help.floodmodeller.com/isis/ISIS/River_Section.htm (Eq. 4) +n_co = definition_dict["n_co"] # Manning's 'n' coefficients +# TODO: use weighted mean 'n' values. See +# http://help.floodmodeller.com/isis/ISIS/River_Section.htm (Eq. 4) # Note: weighted mean calculation requires roughness map. height_data = definition_dict["height_data"] # topography markers = definition_dict["markers"] # distances from corner point @@ -200,7 +199,7 @@ 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, + n_co[p], xregion[p], zregion[p], zmin[p], -- cgit