From fe0db16806da76a18b4f9a567f3e044c6a166b91 Mon Sep 17 00:00:00 2001
From: Paul Garlick <pgarlick@tourbillion-technology.com>
Date: Mon, 14 Oct 2019 12:07:52 +0100
Subject: calculate maximum height in channel.

---
 makeBoundary.py | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

(limited to 'makeBoundary.py')

diff --git a/makeBoundary.py b/makeBoundary.py
index 358bd29..d8c4c98 100755
--- a/makeBoundary.py
+++ b/makeBoundary.py
@@ -110,8 +110,11 @@ for p in range(num_panels):
 # print(zregion)
 # print(xin[12:20])
 
-zmin = zregion[channel].min()      # minimum height
-zmax = zregion[channel].max()-ztol # overtopping height
+# minimum height in channel:
+zmin = zregion[channel].min()
+# overtopping height (minimum of left bank and right bank height):
+zmax = min(zregion[channel][0], zregion[channel][-1]) - ztol
+
 zmax_west = zmax
 zmax_east = zmax
 
-- 
cgit