Calcula a carga poluidora diaria (ou na unidade desejada) como o produto entre concentracao e vazao: \(L = C \times Q \times f\), onde \(f\) e o fator de conversao de unidades.
Usage
compute_load(
df,
param,
flow_col = "vazao",
unit_out = c("kg_dia", "t_dia", "kg_ano", "g_s"),
unit_factor = NULL,
col_name = NULL
)Arguments
- df
Data frame com ao menos a coluna de
parame deflow_col.- param
Character; nome da coluna de concentracao (mg/L).
- flow_col
Character; nome da coluna de vazao. Default
"vazao".- unit_out
Character; unidade de saida. Uma de
"kg_dia"(default),"t_dia","kg_ano","g_s".- unit_factor
Numeric; fator de conversao personalizado. Se fornecido, sobrepoe
unit_out.- col_name
Character; nome da coluna de saida. Default: composto automaticamente como
"{param}_carga_{unit_out}".
Details
Unidades de saida suportadas:
"kg_dia"mg/L * m3/s -> kg/dia (fator = 86.4).
"t_dia"mg/L * m3/s -> t/dia (fator = 0.0864).
"kg_ano"mg/L * m3/s -> kg/ano (fator = 31536).
"g_s"mg/L * m3/s -> g/s (fator = 1).
As unidades pressupõem vazao em m³/s e concentracao em mg/L.
Ajuste o fator com unit_factor se necessario.
See also
Other load-tools:
exceedance_prob()
Examples
# \donttest{
data("wq_demo", package = "tikatuwq")
# Adiciona vazao ficticia para demonstracao
d <- wq_demo
d$vazao <- runif(nrow(d), 2, 10)
d <- compute_load(d, param = "p_total", flow_col = "vazao", unit_out = "kg_dia")
head(d[, c("ponto","p_total","vazao","p_total_carga_kg_dia")])
#> # A tibble: 6 × 4
#> ponto p_total vazao p_total_carga_kg_dia
#> <chr> <dbl> <dbl> <dbl>
#> 1 FBS-BRH-250 0.02 5.11 8.83
#> 2 FBS-BRH-250 0.04 9.80 33.9
#> 3 FBS-BRH-250 0.02 4.32 7.46
#> 4 FBS-BRH-250 0.08 7.43 51.3
#> 5 FBS-BRH-250 0.02 7.88 13.6
#> 6 FBS-BRH-300 0.03 3.57 9.25
# }