From b1ac742d3e1b12f2839de566052c621f3641fe97 Mon Sep 17 00:00:00 2001 From: Paul Garlick Date: Mon, 14 Oct 2019 13:27:50 +0100 Subject: read number of height intervals from boundary definition file. --- boundaryDefinition.txt | 3 +++ makeBoundary.py | 3 +-- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/boundaryDefinition.txt b/boundaryDefinition.txt index fd6f3c0..31d2ea5 100644 --- a/boundaryDefinition.txt +++ b/boundaryDefinition.txt @@ -27,3 +27,6 @@ ztol: 0.01 n_co_chan: 0.035 n_co_west: 0.040 n_co_east: 0.040 + +# number of height intervals: +numH: 50 diff --git a/makeBoundary.py b/makeBoundary.py index 646af68..f874f05 100755 --- a/makeBoundary.py +++ b/makeBoundary.py @@ -32,6 +32,7 @@ height_data = definition_dict["height_data"] # topography markers = definition_dict["markers"] # distances from corner point channel = definition_dict["channel"] # identifier of channel panel ztol = definition_dict["ztol"] # tolerance in overtopping height +numH = definition_dict["numH"] # number of height intervals # print(len(markers)) @@ -120,8 +121,6 @@ print(zmin[1]) # overtopping height (minimum of left bank and right bank heights): zmax = min(zregion[channel][0], zregion[channel][-1]) - ztol -numH = 50 # number of height intervals - def conveyance(numH, n_co, xregion, zregion, zmin, zmax): p_i = [] # wetted perimeter A_i = [] # area -- cgit