diff options
| -rwxr-xr-x | makeBoundary | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/makeBoundary b/makeBoundary index cb27fd7..fc8390c 100755 --- a/makeBoundary +++ b/makeBoundary @@ -301,7 +301,7 @@ h_extra = h_i[panel[ind_p]][ind_q-1] \ print('heights:', h_i[panel[ind_p]][ind_q-1], h_extra, h_i[panel[ind_p]][ind_q]) # find area at target flow by linear interpolation A_extra = A_i[panel[ind_p]][ind_q-1] \ - + (h_extra-h_i[panel[ind_p]][ind_q-1])\ + + (h_extra - h_i[panel[ind_p]][ind_q-1])\ *(A_i[panel[ind_p]][ind_q] - A_i[panel[ind_p]][ind_q-1])\ /(h_i[panel[ind_p]][ind_q] - h_i[panel[ind_p]][ind_q-1]) print('hydraulic radii:', r_h[panel[ind_p]][ind_q-1], r_h[panel[ind_p]][ind_q]) |