aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPaul Garlick <pgarlick@tourbillion-technology.com>2019-10-14 13:27:50 +0100
committerPaul Garlick <pgarlick@tourbillion-technology.com>2019-10-14 13:27:50 +0100
commitb1ac742d3e1b12f2839de566052c621f3641fe97 (patch)
treef5f52e6b4b0c96bed552ff273da6dda08c639e06
parente1a2c82bb7c4e2faa3e391b3dd6d9a051a0e8dd9 (diff)
downloadfullSWOF-utils-b1ac742d3e1b12f2839de566052c621f3641fe97.tar.gz
read number of height intervals from boundary definition file.
-rw-r--r--boundaryDefinition.txt3
-rwxr-xr-xmakeBoundary.py3
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