lsc: Linear storage cascade, unit hydrograph
In maxdschneider/golmtutorial:

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

Optimize the parameters for unit hydrograph as in the framework of the linear storage cascade. Plot observed & simulated data

lsc(P, Q, area = 50, Qbase = Q[1], n = 2, k = 3, x = 1:length(P),
  fit = 1:length(Q), plot = TRUE, main = "Precipitation and discharge",
  plotsim = TRUE, returnsim = FALSE, type = c("o", "l"),
  legx = "center", legy = NULL, ...)

`P`	Vector with precipitation values in mm in hourly spacing
`Q`	Vector with observed discharge (runoff) in m^3/s with the same length as precipitation.
`area`	Single numeric. Catchment area in km^2
`Qbase`	Baseflow that is added to UH-induced simulated Q DEFAULT: Q[1] thus cutting off baseflow in a very simple manner.
`n`	Numeric. Initial number of storages in cascade. Not necessarily integer. DEFAULT: 2
`k`	Numeric. Initial storage coefficient (resistance to let water run out). High damping, slowly reacting landscape, high k. DEFAULT: 3
`x`	Vector for the x-axis of the plot. DEFAULT: sequence along P
`fit`	Integer vector. Indices for a subset of Q that Qsim is fitted to. DEFAULT: all of Q
`plot`	Logical. plot input data? DEFAULT: TRUE
`main`	Character string. DEFAULT: "Precipitation and discharge"
`plotsim`	Logical. add best fit to plot? DEFAULT: TRUE
`returnsim`	Logical. Return simulated Q instead of parameters of UH? DEFAULT: FALSE
`type`	Vector with two characters: type as in `plot`, repeated if only one is given. 1st for obs, 2nd for sim. DEFAULT: c("o","l")
`legx`	legend position. DEFAULT: "center"
`legy`	legend position. DEFAULT: NULL
`...`	Arguments passed to optim

Either vector with optimized n and k and the Nash-Sutcliffe Index, or simulated discharge, depending on the value of returnsim

Berry Boessenkool, berry-b@gmx.de, July 2013

http://ponce.sdsu.edu/onlineuhcascade.php
Skript 'Abflusskonzentration' zur Vorlesungsreihe Abwasserentsorgung I von Prof. Krebs an der TU Dresden
http://tu-dresden.de/die_tu_dresden/fakultaeten/fakultaet_forst_geo_und_hydrowissenschaften/fachrichtung_wasserwesen/isiw/sww/lehre/dateien/abwasserbehandlung/uebung_ws09_10/uebung_awe_1_abflusskonzentration.pdf
http://www.uni-potsdam.de/fs-g3/file.php?fileserver=klausuren&file=%2FMaster_of_Science%2FHydroII_Lernzettel.pdf

unitHydrograph, superPos, nse, rmse. deconvolution.uh in the package hydromad, http://hydromad.catchment.org

# Simplified for this course. True documentation is in berryFunctions::lsc
 # this will not run directly, as the debugging exercise requires to fix the bugs
data(calib)
lsc(calib$P, calib$Q, area=1.6, quietNA=T)