TIMEINTERVAL(min) 60
MUSKINGUMK 5
MUSKINGUMX 0.3
NUMBERREACHES 1
NUMBERFLOWS 287
INITIALSUBREACHOU FINALSUBREACH INFLOW(m3/s) OUTFLOW(m 3/s)
TFLOW(m 3/s) OUTFLOW(m 3/s)
0 1.18317E-08 0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0 0
0.023269128 -0.005817282
0.080549348 -0.012865735
0.176426773 -0.01348132
0.260915868 0.012873429
0.315395175 0.061264212
0.33413975 0.120110809
0.316534008 0.17801948
0.267748866 0.224844398
0.211392023 0.249659725
0.158639433 0.253280947
0.116782705 0.240084751
0.093101135 0.215179632
0.083851964 0.1869723
0.081697415 0.161730853
0.082011364 0.141644007
0.082746988 0.12655194
0.082873183 0.115569153
0.086314306 0.10653488
0.097153225 0.098770007
0.118438629 0.09304446
0.145621341 0.092597325
0.171015777 0.09950472
0.188704624 0.112960273
0.195066699 0.130305842
0.192682844 0.147092019
0.183527551 0.160778549
0.167068336 0.170580603
0.141276847 0.176150409
0.111185494 0.174954856
0.081975968 0.166314898
0.057424883 0.151367936
0.039426705 0.132381717
0.027479732 0.112129708
0.019226271 0.093030579
0.013335119 0.07605229
0.009287645 0.061384866
0.006442101 0.049071946
0.004455299 0.038911185
0.003040992 0.030650791
0.002020659 0.024003424
0.001311197 0.018685098
0.00083104 0.014461662
0.000510936 0.011134033
0.000280091 0.00853597
0.000123117 0.006511244
3.07793E-05 0.004937296
0 0.003718362
0 0.002788772
0 0.002091579
0 0.001568684
0 0.001176513
0 0.000882385
0 0.000661789
0 0.000496341
0 0.000372256
0 0.000279192
0 0.000209394
0 0.000157046
0 0.000117784
0 8.83381E-05
0 6.62536E-05
0 4.96902E-05
0 3.72676E-05
0 2.79507E-05
0 2.0963E-05
0 1.57223E-05
0 1.17917E-05
0 8.84379E-06
0 6.63284E-06
0 4.97463E-06
0 3.73097E-06
0 2.79823E-06
0 2.09867E-06
0 1.574E-06
0 1.1805E-06
0 8.85377E-07
0 6.64033E-07
0 4.98025E-07
0 3.73518E-07
0 2.80139E-07
0 2.10104E-07
0 1.57578E-07
0 1.18184E-07
0 8.86377E-08
0 6.64783E-08
0 4.98587E-08
0 3.7394E-08
0 2.80455E-08
0 2.10341E-08
0 1.57756E-08
0 1.18317E-08
Instruction
1. Place "LibHydro.dll" in C:\temp\
2. Input the values for yellow boxes
3. Push the button "Muskingum" to calculate outflows
4. Push the button "Graph" to draw a graph
Muskingum
Graph
6.1 ROUTEMUSKINGUM Muskingum Routing
Purpose:
Performs hydrograph routing using Muskingum method.
Calling Sequence:
CALL RouteMuskingum (timeInterval, muskingumK, muskingumX, numberReaches,
initialSubreachOutflow, finalSubreachOutflow, numberFlows, inflow, outflow, errorCode,
errorMessage)
Declarations:
INTEGER TIMEINTERVAL
DOUBLE PRECISION MUSKINGUMK. MUSKINGUMX
INTEGER NUMBERREACHES
DOUBLE PRECISION INITIALSUBREACHOUTFLOW(NUMBERREACHES)
DOUBLE PRECISION FINALSUBREACHOUTFLOW(NUMBERREACHES)
INTEGER NUMBERFLOWS
DOUBLE PRECISION INFLOW(NUMBERFLOWS), OUTFLOW(NUMBERFLOWS)
INTEGER*4 ERRORCODE
CHARACTER ERRORMESSAGE*80
Augument Description:
TIMEINTERVAL (Input) Time interval, minutes
MUSKINGUMK (Input) Muskingum K, hours. muskingumK is divided by
numberReaches to obtain travel time for each subreach.
MUSKINGUMX (Input) Muskingum x
NUMBERREACHES (Input) Number of routing reaches. muskingumK is divided by
this value and routing is repeated for each reach.
INITIALSUBREACHOUTFLOW (Input) Initial outflow for each subreach, m3/s
FINALSUBREACHOUTFLOW (Output) Final outflow for each subreach, m3/s.
NUMBERFLOWS (Input) Number ofhydrograph ordinates.
INFLOW (Input) Inflow hydrograph ordinates, m3/s.
OUTFLOW (Output) Routed hydrograph ordinates, m3/s.
ERRORCODE (Output) Coded integer indicating error number, severity and
ERRORMESSAGE (Output) Errormessage. See Remarks.
Remarks:
Flows are routed using the formula Oi=(ca-cb)Ii-1+cbIi+(1-ca)Oi-1 where
ca=2Δt/2K(1-x)+Δt, and Δt-2Kx/(2K(1-x)+Δt
Initial outflow is required for each subreach. Initial inflow for the first subreach in the fIrst hydrograph
ordinate. Initial inflow for each subsequent subreach is the initial outflow from the subreach
immediately upstream.
For error codes 1 and 2, outflow is set to undefmedDbl.
Muskingum routing is stable when 1/2(1-K)≤K/Δt≤1/2x
Error severity errorMessage
0 0
11 4 Muskingum routing is unstable with the given .parameters.
12 4 Invalid Muskingum parameters: K = xxx, x = xxx, Number of reaches = nnn
Subroutines Called: