onlinerouting05:   Muskingum-Cunge routing   ⇒   D = qo/(SocΔx)


Formulas

V = Qp/Ap         qo = Qp/Tp         c = βV

C = c(Δt/Δx)         D = qo/(SocΔx)

C0 = (-1+C+D)/(1+C+D)         C1 = (1+C-D)/(1+C+D)         C2 = (1-C+D)/(1+C+D)

O2 = C0I2 + C1I1 + C2O1

Reference


SAMPLE INPUT DATA:

Peak discharge Qp (m3/s):   1000        Peak flow area Ap (m2):   400        Peak top width Tp (m):   100

Time interval Δt (hr):   1        Reach length Δx (km):   14.4        Rating exponent β :   1.6       

Bottom slope So (m/m):   0.000868        Number of inflow hydrograph ordinates n:   11

Inflow hydrograph ordinates Qi:      0, 200,400,600,800,1000,800,600,400,200,0


SAMPLE OUTPUT:

Time interval Time (hr) Time
(days)
Inflow
(m3/s)
C0I2 C1I1 C2O1 Outflow
(m3/s)
0 0 0 0 - - - 0
1 1 0.042 200 18.183 0 0 18.183
2 2 0.083 400 36.366 163.634 1.653 201.653
3 3 0.125 600 54.549 327.268 18.333 400.15
4 4 0.167 800 72.732 490.903 36.379 600.014
5 5 0.208 1000 90.914 654.537 54.550 800.001
6 6 0.250 800 72.732 818.171 72.732 963.634
7 7 0.292 600 54.549 654.537 87.608 796.694
8 8 0.333 400 36.366 490.903 72.431 599.699
9 9 0.375 200 18.183 327.268 54.521 399.973
10 10 0.417 0 0 163.634 36.363 199.998
11 11 0.458 0 0 0 18.183 18.183
12 12 0.500 0 0 0 1.653 1.653
13 13 0.542 0 0 0 0.150 0.150
14 14 0.583 0 0 0 0.014 0.014
15 15 0.625 0 0 0 0.001 0.001
16 16 0.667 0 0 0 0 0
17 17 0.708 0 0 0 0 0
18 18 0.75 0 0 0 0 0
19 19 0.792 0 0 0 0 0
20 20 0.833 0 0 0 0 0


            


Thank you for running onlinerouting05.   Please call again.   [070423]