pdr_predict {PoolDilutionR}

Predict total pool, heavy isotope pool, and atom percent

Description

Predict total pool, heavy isotope pool, and atom percent

Usage

pdr_predict(
  time,
  m0,
  n0,
  P,
  k,
  pool = "CH4",
  frac_P = frac_P_default(pool),
  frac_k = frac_k_default(pool)
)

Arguments

time

Vector of numeric time values (e.g. days); first should be zero

m0

total pool size at time zero, as a volume

n0

pool size of heavy isotope at time zero, as a volume

P

production rate, unit gas/unit time

k

first-order rate constant for consumption, 1/unit time

pool

Name of pool; see pdr_fractionation

frac_P

Fractionation value for production; see pdr_fractionation

frac_k

Fractionation value for consumption; see pdr_fractionation

Value

Returns a data frame with mt, nt, and AP_pred (atom percent) for each time step

Note

This is Eq. 11 from von Fischer and Hedin 2002 with a few modifications.

Author(s)

K.A. Morris & B. Bond-Lamberty

Examples

pdr_predict(time = 0:5, m0 = 10, n0 = 1, P = 0.5, k = 0.3)