# NOT RUN {
var.values<-ext.output(dycd.output=system.file("extdata", "dysim.nc", package = "dycdtools"),
var.extract=c("TEMP"))
for(i in 1:length(var.values)){
expres<-paste0(names(var.values)[i],"<-data.frame(var.values[[",i,"]])")
eval(parse(text=expres))
}
# interpolate temperature for depths from 0 to 13 m at increment of 0.5 m
temp.interpolated<-interpol(layerHeights = dyresmLAYER_HTS_Var,
var = dyresmTEMPTURE_Var,
min.dept = 0,max.dept = 13,by.value = 0.5)
data(obs_temp)
# time series plot of temperature sim and obs
plot_ts(sim = temp.interpolated,
obs = obs_temp,
target.depth=c(1,6),
sim.start="2017-06-06",
sim.end="2017-06-15",
plot.start="2017-06-06",
plot.end="2017-06-15",
ylabel="Temperature \u00B0C",
min.depth=0,
max.depth=13,
by.value=0.5)
# }
Run the code above in your browser using DataLab