Calculate the NLC50 and PLC50 for a set of PAHs — plc50_multi • pahwq

Given the results of a TUV model run (via tuv() or run_tuv()), get the NLC50, Pabs, and PLC50 for a set of PAHs

Usage

plc50_multi(tuv_results, pahs, time_multiplier = 2, ...)

Arguments

tuv_results: data.frame of TUV results
pahs: names of PAHs for which to calculate NLC50, Pabs, and PLC50
time_multiplier: multiplier to get the total exposure time. I.e., if the tuv_results contains 24 hours of data, and you need a 48h exposure, the multiplier would be 2. (this is the default)
...: arguments passed on to nlc50()

Value

a data.frame of NLC50, Pabs, and PLC50 for the given PAHs and TUV results

Examples

tuv_res <- tuv(
  depth_m = 0.25,
  lat = 49.601632,
  lon = -119.605862,
  DOC = 5,
  date = "2023-06-21"
)

plc50_multi(tuv_res, c("Anthracene", "fluorene", "pyrene"))
#>          pah     nlc50         pabs      plc50
#> 1 anthracene  58.40685 1141.8593777  1.9545602
#> 2   fluorene 111.26696    0.3084483 71.1781565
#> 3     pyrene  18.88737  513.4631473  0.9127577