-
Notifications
You must be signed in to change notification settings - Fork 0
/
example_UH_Mykonos.py
44 lines (35 loc) · 8.54 KB
/
example_UH_Mykonos.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
from hydrograph import SnyderUH
'''
This is a example by zhanliang Roger
'''
# Actually this name should be the catchment's name XD
name= "Roger"
# descriptions for variables
'''
A: Area
L: River Length
Lc: Distance from catchment center to basin exit point
C1: Unit correction factor
Ct: Unit Hydrograph regional parameter
Cp: Unit Hydrograph regional parameter
tR: Effective rainfall duration
'''
A = 2.7
L = 2.10
Lc = 1.37
C1 = 1
Ct = 1.9
Cp = 0.65
tR = 0.25
#Make ShyderUH
mykonosUH = SnyderUH(name, A, L, Lc, C1, Ct, Cp, tR)
# Calculate ShyderUh
print(mykonosUH.calc())
# View plot
mykonosUH.plot()
# Calculate UH & View plot for various rainfall durations
for item in [0.25,0.5,1,2,3,4]:
UH = SnyderUH(name, A, L, Lc, C1, Ct, Cp, item)
UH.calc()
UH.plot()
print("tR: {} hr, Q: {}, Tb:{} hr".format(UH.tR, UH.QPR, UH.Tb))