mHM_plotQ_d: Plot daily river flow.

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

Description

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

Usage

mHM_plotQ_d(ts, 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 = NA, outfile = "out.pdf", basinid,
  ylims = NULL)

Arguments

ts	zoo object, simulated and observed 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, default = NA no aggregation, see zoo::rollmean
outfile	output pdf
basinid	character, basin id
ylims	numeric vector for y axes limitations, c(min,max)

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/validation 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.