The simple abc hydrological Model

Share:

Description

Some description todo

Usage

1
abc(pars, init, prec, R.implementation = FALSE)

Arguments

pars

A vector with the three named (!) parameters a,b,c

init

A vector with the named (!) initial condition for the storage G

prec

The precipitation data as vector. ToDo: change to xts

R.implementation

boolean, determining if R or Fortran code is run.

Details

Give some more details about the model here. Todo These are our model equations (Q: stream flow, G: base flow) Q[t] = (1-a-b) * P[t] + c * G[t-1] G[t] = (1-c) * G[t-1] + a * P[t] Rekursionsvorschrift fur Gt: Gt = G0*(1-c)^t + Summe_von_i=1_bis_t_uber(a*Pi*(1-c)^(n-t)) die Zeitreihe fuer G wird um einen Zeitindex nach "rechts" verschoben (also G[1]=G0) auf diese Weise kann dann Q ohne Schleife berechnet werden

Value

q

The vector of simulated discharges

G

The vector of changes in the store G

Author(s)

Maik Heistermann (?), Dominik Reusser

References

Vogel and Sankarasubramanian, WRR, 2003

Examples

1
2
3
4
data(huagrahuma)
abc(pars=c(a=0.5,b=0.3,c=0.7), init=c(G=1), prec = huagrahuma$inputs$P, R.implement=TRUE)
#Todo: make more sensible example
abc(pars=c(a=0.5,b=0.3,c=0.7), init=c(G=1), prec = huagrahuma$inputs$P)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.