routing_muskingum.Rdrouting_muskingum
routing_muskingum(I1, I2, Q1, deltaT, par)
routing_muskingum_nL(I, Q0, deltaT, par)
calib_routing_muskingum_nL(I, Q0, yobs, maxn = 1000)in hours
X : 0-1/2
K : in the same unit as deltaT
nt : divide river into nt chunks (default 1)