#' @title test
#'
#' @description test
#'
#' @param test test
#'
#' @return test
#'
#' @details missing
#' @references Marvin Reich (2017), mreich@@posteo.de
#' @examples missing
plot_waterDistributionOverTime <- function(
soilmoisture_mod,
output_dir = dir_output,
...
){
## load data
sm_mod_data = load(file = paste0(output_dir, soilmoisture_mod))
sm_mod_data = get(sm_mod_data)
# format data
data_plot = dplyr::group_by(sm_mod_data, datetime, Depth) %>%
dplyr::summarize(value = max(value, na.rm = T))
# get value ranges
data_qt = round(stats::quantile(data_plot$value, na.rm = T, names = F), 2)
# maximal vertical extent
# min has to be used as maximum extents are given in negative distances (below surface)
maxy = min(data_plot$Depth, na.rm=T)
# miny = max(data_plot$Depth, na.rm=T)
# actual plotting
data.gg = ggplot(data_plot, aes(x = datetime, y = Depth)) +
geom_tile(aes(fill = value)) +
ylim(maxy,0) +
# theme(legend.position=leg) +
ylab("Depth [m]") + xlab("Time since start of experiment") +
# scale_fill_gradientn(colours = rev(viridis(7))) +
scale_fill_gradientn(breaks = data_qt, colours=rev(viridis(7)), na.value="red") +
labs(fill = expression(Delta * "Soil moisture [%VWC]")) +
theme(legend.position ="bottom",
# legend.text=element_text(size=17),
# legend.title=element_text(size=19),
panel.grid.major = element_line(colour = "black", linetype = "dotted"),
panel.grid.minor = element_line(colour = "black", linetype = "dotted"))
return(data.gg)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.