Title: | R Package for Generic Calculation of Attenuation of PMT & Other Substances During Riverbank Filtration |
---|---|
Description: | R package for generic calculation of attenuation of PMT & other substances during riverbank filtration. |
Authors: | Dwight Baldwin [aut, cre] , Michael Rustler [ctb] , PROMISCES [fnd], Kompetenzzentrum Wasser Berlin gGmbH (KWB) [cph] |
Maintainer: | Dwight Baldwin <[email protected]> |
License: | MIT + file LICENSE |
Version: | 0.1.0 |
Built: | 2024-11-20 10:37:43 UTC |
Source: | https://github.com/KWB-R/kwb.1dbear |
Calculate Bear 1D
calculate_bear1d( n_values, rs_values, foc_values, log_koc_values, hl_values, C0 = 1, D_values, v_values, t_values = max(hl_values * max(log_koc_values)), x_values )
calculate_bear1d( n_values, rs_values, foc_values, log_koc_values, hl_values, C0 = 1, D_values, v_values, t_values = max(hl_values * max(log_koc_values)), x_values )
n_values |
Porosity (-) |
rs_values |
density solids (g/g) |
foc_values |
fraction of organic matter |
log_koc_values |
min, max Koc for relevant compounds Nguyen et al 2020 |
hl_values |
half life (T) |
C0 |
initial concentration of the solute (M/L3) |
D_values |
coefficient of longitudinal dispersion (L2/T) |
v_values |
average linear ground water velocity (L/T) |
t_values |
time (T), default: max(hl_values * max(log_koc_values)) |
x_values |
flow path distance (L) |
Bear 1D results
bear1d_vienna <- calculate_bear1d(n_values = 0.15, rs_values = 2.7, foc_values = 0.002, log_koc_values = c(1,2.9), hl_values = c(500,2000), C0 = 1, D_values = 10, v_values = 21.3, t_values = seq(0, 350, by=1), x_values = 141 )
bear1d_vienna <- calculate_bear1d(n_values = 0.15, rs_values = 2.7, foc_values = 0.002, log_koc_values = c(1,2.9), hl_values = c(500,2000), C0 = 1, D_values = 10, v_values = 21.3, t_values = seq(0, 350, by=1), x_values = 141 )
Calculate Cx
calculate_Cx(x, t, v, D, retardation, C0, k)
calculate_Cx(x, t, v, D, retardation, C0, k)
x |
x |
t |
t |
v |
v |
D |
D |
retardation |
retardation |
C0 |
initial concentration |
k |
k |
????
Plot Combined
plot_combined(bear1d_list, text_size = 20, line_width = 1.25, xbreaks = 200)
plot_combined(bear1d_list, text_size = 20, line_width = 1.25, xbreaks = 200)
bear1d_list |
Bear1D list as retrieved by |
text_size |
size of plot title and axis labels, default: 20 |
line_width |
line width, default: 1.25 |
xbreaks |
breaks on x axis, default: 200 |
combined plot
Plot Heatmap
plot_heatmap(bear1d_list, rect_data = NULL, measured_C_C0 = NULL)
plot_heatmap(bear1d_list, rect_data = NULL, measured_C_C0 = NULL)
bear1d_list |
Bear1D list as retrieved by |
rect_data |
data.frame with xmin/xmax/ymin/ymax parameters (default: NULL) |
measured_C_C0 |
in percent (default: NULL) |
heatmap plot