datC<-
cbind(
P = c(P_E$P),
Q=as.numeric(Qtrue),
E = c(P_E$E))
datCz1 <- zoo(datC, as.Date("2010-12-31") + 1:nrow(datC))
# split the data to use for calibration
Data_Cal1<- window(datCz1, start = "2011-01-01",end = "2017-12-31")
#defining the parameter ranges in the object
modCz1 <- hydromad(Data_Cal1,
sma="sacramento", uztwm=c(1.0,150.0), uzfwm=c(1.0,150.0), uzk=c(0.1,0.5),
pctim=c(0.0,0.1),adimp=c(0.0,0.40), zperc=c(1.0,250.0), rexp=c(0.0,5.0),
lztwm=c(1.0,500.0), lzfsm=c(1.0,1000), lzfpm=c(1.0,1000), lzsk=c(0.01,0.25),
lzpk=c(0.0001,0.025),
pfree=c(0.0,0.6))