mHM_plotQmulti_d: Plot daily river flow of multiple simulations.

In JBrenn/mHMr: pre- and post-prossesing for mHM

Description

mHM_plotQmulti_d is ploting a comparison of daily river flow data - simulation vs. observation in pdf file. Additionally GOFs (KGE, NSE) for cross-validation simulations are calculated, if calibration period is given.

Usage

mHM_plotQmulti_d(ts1, ts2, windows = c(start = as.Date("1989-10-01"), end =
  as.Date("2008-03-30")), calibTime = c(start = as.Date("2000-01-01"), end =
  as.Date("2008-03-30")), validTime = c(start = as.Date("1990-01-01"), end =
  as.Date("1999-12-31")), rollsteps = 10, outfile = "out.pdf",
  ylims = NULL, cols = c(rgb(1, 0, 0, 0.75), grey.colors(1, 0.2, 0.2, alpha
  = 0.75), rgb(0, 0, 1, 0.7)))

Arguments

ts1	zoo object, observed and simulated river flow data, as retrived from mhm_readQ.
ts2	zoo object, observed and simulated river flow data, as retrived from mhm_readQ.
windows	vector of Date objects, defining start and end of series being ploted.
calibTime	vector of Date objects, defining start and end of calibration period.
validTime	vector of Date objects, defining start and end of validation period.
rollsteps	integer width of the rolling window, see zoo::rollmean
outfile	output pdf
ylims	numeric vector for y axes limitations, c(min,max)
cols	color definition, lenght = 3

Author(s)

Johannes Brenner [e-mail]

Examples

data("Qout_mHM")

# devide calibration and validation period
gof <- mHM_plotQ_d(ts = Qout_mHM, windows = c(start=as.Date("1989-10-01"), end=as.Date("2008-03-30")), 
                   calibTime = c(start=as.Date("2000-01-01"), end=as.Date("2008-03-30")),
                   validTime = c(start=as.Date("1990-01-01"), end=as.Date("1999-12-31"))
                   rollsteps = 10)
                   
# do not consider calibration/validatio period                  
gof <- mHM_plotQ_d(ts = Qout_mHM, windows = c(start=as.Date("1989-10-01"), end=as.Date("2008-03-30")), 
                   calibTime = NULL, validTime = NULL,
                   rollsteps = 10)
                   

JBrenn/mHMr documentation built on May 7, 2019, 7:39 a.m.