reachRouting: channel routing computation

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/reachRouting.R

Description

function for flood routing using parameteric Muskingum and muskingum-cunge techniques.

Usage

1
2
reachRouting(inflow,routingMethod,
             routingParams,simulation)

Arguments

inflow

a vector of runoff (cms) presenting a runoff event generated by excess rainfall computed by loss methods or an object inherited from any of the following classes :transform ; reachRouting ; reservoirRouting.

routingMethod

a string: the type of channel routing method: "muskingum" or "muskingumcunge". default to "muskingum"

routingParams

a list : parameters associated to the routingMethod:

  • k and x for "muskingum",

  • bedWith (m), sideSlope (m/m), channelSlope (m/m), manningRoughness, riverLength (Km) for "muskingumcunge"

simulation

a list of simulation time and dates as below:

  • start: the date which simulation starts, must be in 'YYYY-MM-DD' format

  • start: the date which simulation ends, must be in 'YYYY-MM-DD' format

  • by: the interval of each steps in seconds

Value

a data.frame: including inflow time series routing resaults and simulation details

Author(s)

Rezgar Arabzadeh

References

Chow, V. T., Maidment, D. R., & Mays, L. W. (1988). Applied hydrology.

See Also

reservoirRouting

Examples

1
2
3
4
5
6
7
8
9
inflow<-c(100,500,1500,2500,5000,11000,22000,28000,28500,26000,
          22000,17500,14000,10000,7000,4500,2500,1500,1000,500,100)
routingMethod<-c("muskingum","muskingumcunge")
routingParams<-list(k=3,x=0.2,bedWith=50,sideSlope=2,channelSlope=0.0001,
                    manningRoughness=0.01,riverLength=100)
simulation<-list(start='2000-01-01',end='2000-01-04',by=3600)

reachRouting(inflow,routingMethod[1],routingParams,simulation)
reachRouting(inflow,routingMethod[2],routingParams,simulation)

Example output

Loading required package: pso
Loading required package: Hmisc
Loading required package: lattice
Loading required package: survival
Loading required package: Formula
Loading required package: ggplot2

Attaching package: 'Hmisc'

The following objects are masked from 'package:base':

    format.pval, units

Loading required package: network
network: Classes for Relational Data
Version 1.14-377 created on 2019-03-04.
copyright (c) 2005, Carter T. Butts, University of California-Irvine
                    Mark S. Handcock, University of California -- Los Angeles
                    David R. Hunter, Penn State University
                    Martina Morris, University of Washington
                    Skye Bender-deMoll, University of Washington
 For citation information, type citation("network").
 Type help("network-package") to get started.


Attaching package: 'network'

The following object is masked from 'package:Hmisc':

    is.discrete

Loading required package: GGally

Attaching package: 'RHMS'

The following objects are masked from 'package:base':

    transform, transform.default

$operation
                        I            O
2000-01-01 12:00:00   100 0.000000e+00
2000-01-01 13:00:00   500 2.068966e+01
2000-01-01 14:00:00  1500 1.514863e+02
2000-01-01 15:00:00  2500 5.820083e+02
2000-01-01 16:00:00  5000 1.157178e+03
2000-01-01 17:00:00 11000 2.275392e+03
2000-01-01 18:00:00 22000 4.904567e+03
2000-01-01 19:00:00 28000 1.059265e+04
2000-01-01 20:00:00 28500 1.657794e+04
2000-01-01 21:00:00 26000 2.077520e+04
2000-01-01 22:00:00 22000 2.271479e+04
2000-01-01 23:00:00 17500 2.262348e+04
2000-01-02 00:00:00 14000 2.097745e+04
2000-01-02 01:00:00 10000 1.870937e+04
2000-01-02 02:00:00  7000 1.580958e+04
2000-01-02 03:00:00  4500 1.285800e+04
2000-01-02 04:00:00  2500 1.004490e+04
2000-01-02 05:00:00  1500 7.477692e+03
2000-01-02 06:00:00  1000 5.433661e+03
2000-01-02 07:00:00   500 3.922053e+03
2000-01-02 08:00:00   100 2.755828e+03
2000-01-02 09:00:00     0 1.843474e+03
2000-01-02 10:00:00     0 1.207793e+03
2000-01-02 11:00:00     0 7.913127e+02
2000-01-02 12:00:00     0 5.184462e+02
2000-01-02 13:00:00     0 3.396717e+02
2000-01-02 14:00:00     0 2.225435e+02
2000-01-02 15:00:00     0 1.458044e+02
2000-01-02 16:00:00     0 9.552700e+01
2000-01-02 17:00:00     0 6.258666e+01
2000-01-02 18:00:00     0 4.100505e+01
2000-01-02 19:00:00     0 2.686538e+01
2000-01-02 20:00:00     0 1.760145e+01
2000-01-02 21:00:00     0 1.153199e+01
2000-01-02 22:00:00     0 7.555440e+00
2000-01-02 23:00:00     0 4.950116e+00
2000-01-03 00:00:00     0 3.243179e+00
2000-01-03 01:00:00     0 2.124842e+00
2000-01-03 02:00:00     0 1.392138e+00
2000-01-03 03:00:00     0 9.120902e-01
2000-01-03 04:00:00     0 5.975763e-01
2000-01-03 05:00:00     0 3.915155e-01
2000-01-03 06:00:00     0 2.565102e-01
2000-01-03 07:00:00     0 1.680584e-01
2000-01-03 08:00:00     0 1.101072e-01
2000-01-03 09:00:00     0 7.213921e-02
2000-01-03 10:00:00     0 4.726362e-02
2000-01-03 11:00:00     0 3.096582e-02
2000-01-03 12:00:00     0 2.028795e-02
2000-01-03 13:00:00     0 1.329211e-02
2000-01-03 14:00:00     0 8.708621e-03
2000-01-03 15:00:00     0 5.705649e-03
2000-01-03 16:00:00     0 3.738184e-03
2000-01-03 17:00:00     0 2.449155e-03
2000-01-03 18:00:00     0 1.604619e-03
2000-01-03 19:00:00     0 1.051302e-03
2000-01-03 20:00:00     0 6.887840e-04
2000-01-03 21:00:00     0 4.512723e-04
2000-01-03 22:00:00     0 2.956611e-04
2000-01-03 23:00:00     0 1.937090e-04
2000-01-04 00:00:00     0 1.269128e-04
2000-01-04 01:00:00     0 8.314977e-05
2000-01-04 02:00:00     0 5.447743e-05
2000-01-04 03:00:00     0 3.569211e-05
2000-01-04 04:00:00     0 2.338449e-05
2000-01-04 05:00:00     0 1.532087e-05
2000-01-04 06:00:00     0 1.003781e-05
2000-01-04 07:00:00     0 6.576498e-06
2000-01-04 08:00:00     0 4.308740e-06
2000-01-04 09:00:00     0 2.822967e-06
2000-01-04 10:00:00     0 1.849530e-06
2000-01-04 11:00:00     0 1.211761e-06
2000-01-04 12:00:00     0 7.939126e-07

$call
reachRouting.default(inflow = inflow, routingMethod = routingMethod[1], 
    routingParams = routingParams, simulation = simulation)

attr(,"class")
[1] "reachRouting"
$operation
                    inflow        32 Km        64 Km        95 Km
2000-01-01 12:00:00    100 1.000000e+02 1.000000e+02 1.000000e+02
2000-01-01 13:00:00    500 2.517241e+02 1.575505e+02 1.218295e+02
2000-01-01 14:00:00   1500 7.510381e+02 3.924641e+02 2.282003e+02
2000-01-01 15:00:00   2500 1.492358e+03 8.469703e+02 4.799958e+02
2000-01-01 16:00:00   5000 2.927675e+03 1.703348e+03 9.822051e+02
2000-01-01 17:00:00  11000 6.205191e+03 3.538320e+03 2.026792e+03
2000-01-01 18:00:00  22000 1.269517e+04 7.289061e+03 4.180081e+03
2000-01-01 19:00:00  28000 1.946850e+04 1.247129e+04 7.648472e+03
2000-01-01 20:00:00  28500 2.378184e+04 1.748947e+04 1.188302e+04
2000-01-01 21:00:00  26000 2.511408e+04 2.103621e+04 1.593820e+04
2000-01-01 22:00:00  22000 2.402504e+04 2.259416e+04 1.899326e+04
2000-01-01 23:00:00  17500 2.133935e+04 2.226706e+04 2.060967e+04
2000-01-02 00:00:00  14000 1.815602e+04 2.061118e+04 2.078268e+04
2000-01-02 01:00:00  10000 1.462997e+04 1.808702e+04 1.974234e+04
2000-01-02 02:00:00   7000 1.125415e+04 1.513558e+04 1.782273e+04
2000-01-02 03:00:00   4500 8.249642e+03 1.211986e+04 1.538001e+04
2000-01-02 04:00:00   2500 5.678639e+03 9.273989e+03 1.272475e+04
2000-01-02 05:00:00   1500 3.762940e+03 6.809539e+03 1.012204e+04
2000-01-02 06:00:00   1000 2.479497e+03 4.850149e+03 7.777734e+03
2000-01-02 07:00:00    500 1.574730e+03 3.361112e+03 5.797884e+03
2000-01-02 08:00:00    100 9.035374e+02 2.243077e+03 4.195993e+03
2000-01-02 09:00:00      0 4.772182e+02 1.433906e+03 2.945128e+03
2000-01-02 10:00:00      0 2.465557e+02 8.840003e+02 2.006097e+03
2000-01-02 11:00:00      0 1.273835e+02 5.306895e+02 1.329719e+03
2000-01-02 12:00:00      0 6.581295e+01 3.123980e+02 8.607097e+02
2000-01-02 13:00:00      0 3.400239e+01 1.811455e+02 5.458989e+02
2000-01-02 14:00:00      0 1.756740e+01 1.037902e+02 3.402544e+02
2000-01-02 15:00:00      0 9.076235e+00 5.889382e+01 2.089303e+02
2000-01-02 16:00:00      0 4.689255e+00 3.315056e+01 1.266458e+02
2000-01-02 17:00:00      0 2.422713e+00 1.853412e+01 7.591088e+01
2000-01-02 18:00:00      0 1.251700e+00 1.030253e+01 4.505561e+01
2000-01-02 19:00:00      0 6.466931e-01 5.698343e+00 2.651136e+01
2000-01-02 20:00:00      0 3.341153e-01 3.138074e+00 1.548030e+01
2000-01-02 21:00:00      0 1.726214e-01 1.721530e+00 8.977398e+00
2000-01-02 22:00:00      0 8.918517e-02 9.412200e-01 5.174309e+00
2000-01-02 23:00:00      0 4.607770e-02 5.130396e-01 2.965841e+00
2000-01-03 00:00:00      0 2.380614e-02 2.788867e-01 1.691468e+00
2000-01-03 01:00:00      0 1.229949e-02 1.512294e-01 9.602778e-01
2000-01-03 02:00:00      0 6.354555e-03 8.182294e-02 5.428993e-01
2000-01-03 03:00:00      0 3.283094e-03 4.418041e-02 3.057608e-01
2000-01-03 04:00:00      0 1.696217e-03 2.381086e-02 1.716000e-01
2000-01-03 05:00:00      0 8.763539e-04 1.281081e-02 9.599401e-02
2000-01-03 06:00:00      0 4.527700e-04 6.881644e-03 5.353858e-02
2000-01-03 07:00:00      0 2.339245e-04 3.691250e-03 2.977690e-02
2000-01-03 08:00:00      0 1.208576e-04 1.977271e-03 1.651832e-02
2000-01-03 09:00:00      0 6.244132e-05 1.057820e-03 9.141177e-03
2000-01-03 10:00:00      0 3.226043e-05 5.652576e-04 5.047268e-03
2000-01-03 11:00:00      0 1.666742e-05 3.017199e-04 2.780935e-03
2000-01-03 12:00:00      0 8.611255e-06 1.608846e-04 1.529191e-03
2000-01-03 13:00:00      0 4.449022e-06 8.570479e-05 8.393064e-04
2000-01-03 14:00:00      0 2.298596e-06 4.561430e-05 4.598479e-04
2000-01-03 15:00:00      0 1.187574e-06 2.425632e-05 2.515276e-04
2000-01-03 16:00:00      0 6.135627e-07 1.288836e-05 1.373645e-04
2000-01-03 17:00:00      0 3.169984e-07 6.842872e-06 7.490609e-05
2000-01-03 18:00:00      0 1.637779e-07 3.630486e-06 4.078938e-05
2000-01-03 19:00:00      0 8.461617e-08 1.924833e-06 2.218173e-05
2000-01-03 20:00:00      0 4.371711e-08 1.019854e-06 1.204733e-05
2000-01-03 21:00:00      0 2.258653e-08 5.400253e-07 6.535219e-06
2000-01-03 22:00:00      0 1.166938e-08 2.857813e-07 3.541017e-06
2000-01-03 23:00:00      0 6.029008e-09 1.511504e-07 1.916538e-06
2000-01-04 00:00:00      0 3.114899e-09 7.990093e-08 1.036216e-06
2000-01-04 01:00:00      0 1.609319e-09 4.221547e-08 5.596886e-07
2000-01-04 02:00:00      0 8.314579e-10 2.229352e-08 3.020124e-07
2000-01-04 03:00:00      0 4.295744e-10 1.176744e-08 1.628182e-07
2000-01-04 04:00:00      0 2.219405e-10 6.208546e-09 8.769957e-08
2000-01-04 05:00:00      0 1.146660e-10 3.274242e-09 4.719804e-08
2000-01-04 06:00:00      0 5.924242e-11 1.726045e-09 2.538031e-08
2000-01-04 07:00:00      0 3.060771e-11 9.095377e-10 1.363736e-08
2000-01-04 08:00:00      0 1.581354e-11 4.790971e-10 7.322122e-09
2000-01-04 09:00:00      0 8.170096e-12 2.522707e-10 3.928522e-09
2000-01-04 10:00:00      0 4.221097e-12 1.327873e-10 2.106292e-09
2000-01-04 11:00:00      0 2.180838e-12 6.987117e-11 1.128538e-09
2000-01-04 12:00:00      0 1.126734e-12 3.675335e-11 6.042715e-10

$call
reachRouting.default(inflow = inflow, routingMethod = routingMethod[2], 
    routingParams = routingParams, simulation = simulation)

attr(,"class")
[1] "reachRouting"

RHMS documentation built on Sept. 27, 2021, 5:06 p.m.