Description Usage Arguments Value Examples
Wrapper function for sinus optimisation
1 2 3 4 5 6 7 8 9 10 11 12 | run_optimisation(
data_sw_selected,
data_gw_selected,
retardation_factor = 2,
sw_monitoring_id = ifelse(!is.null(attr(data_sw_selected, "monitoring_id")),
attr(data_sw_selected, "monitoring_id"), "surface-water monitoring point"),
gw_monitoring_id = ifelse(!is.null(attr(data_gw_selected, "monitoring_id")),
attr(data_gw_selected, "monitoring_id"), "groundwater monitoring point"),
limits = c(100, 500),
tolerance = 0.001,
debug = FALSE
)
|
data_sw_selected |
data.frame with daily data temperature data of surface water monitoring point with columns "date" (format: "YYYY-MM-DD") and "value" (format: double, temperature in degree Celsius) for selected time period |
data_gw_selected |
data.frame with daily data temperature data of groundwater monitoring point with columns "date" (format: "YYYY-MM-DD") and "value" (format: double, temperature in degree Celsius) for selected time period |
retardation_factor |
hydraulic retardation factor (default: 2) |
sw_monitoring_id |
optional label for surface water monitoring id (default: "surface-water monitoring point"
or attr(data_sw_selected, "monitoring_id") if data imported with |
gw_monitoring_id |
optional label for groundwater monitoring id (default: "surface-water monitoring point"
or attr(data_sw_selected, "monitoring_id") if data imported with |
limits |
minimum/maximum period length for sinus optimisation in days (default: c(100, 500)) |
tolerance |
the desired accuracy (default: 0.001 |
debug |
show debug messages (default: FALSE) |
list with sim/observation data ("data") fit parameters ("paras"), goodness-of-fit values ("gof")
traveltimes ("traveltimes") and special (min, max, turning) points ("points") as returned by get_predictions
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 | load_temp <- function(base_name) {
kwb.heatsine::load_temperature_from_csv(
kwb.heatsine::extdata_file(base_name)
)
}
data_sw <- load_temp("temperature_surface-water_Txxsxx-mxxxxsxxx.csv")
data_gw <- load_temp("temperature_groundwater_Txxxx3.csv")
data_sw_selected <- kwb.heatsine::select_timeperiod(
data_sw,
date_start = "2015-10-10",
date_end = "2016-10-14"
)
data_gw_selected <- kwb.heatsine::select_timeperiod(
data_gw,
date_start = "2015-12-28",
date_end = "2016-12-26"
)
kwb.heatsine::run_optimisation(data_sw_selected = data_sw_selected,
data_gw_selected = data_gw_selected,
retardation_factor = 1.8,
sw_monitoring_id = attr(data_sw_selected, "monitoring_id"),
gw_monitoring_id = attr(data_gw_selected, "monitoring_id"),
limits = c(100, 500),
tolerance = 0.001,
debug = FALSE)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.