# 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)
# scatter plot of sim and obs temperature
plot_scatter(sim=temp.interpolated,
obs=obs_temp,
sim.start="2017-06-06",
sim.end="2017-06-15",
plot.start="2017-06-06",
plot.end="2017-06-15",
plot.save=FALSE,
min.depth = 0,max.depth = 13,by.value = 0.5)
# }
Run the code above in your browser using DataLab