The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
An R package for easy calculation of rates from isotope pool dilution experiments.
To load the package, type library(PoolDilutionR).
This code calculates gross production and consumption from pool dilution time series data. Target pool size (e.g., methane, nitrate) and isotopic composition (e.g., 13C, 15N) are measured from a closed system over time, the model optimizes production rate (P) and the first order rate constant (k), by minimizing error in the model-predicted total pool size (unit volume or mass), as well as the isotopic signature. The optimization weighs the signal:noise ratio of pool and signatures using measurement precision as well as the magnitude of change over time. The calculations used here are based on von Fischer and Hedin 2002, 10.1029/2001GB001448, with some modifications (manuscript link pending).
Original code was written for data collected on a Picarro CRDS equipped with a SSIM2 unit.
Time series data of pool and isotopic composition.
Prediction of P, gross production rate, and k, first order rate constant for consumption.
`pred <- pdr_predict(time = dat$time_days,
m0 = dat$cal12CH4ml[1] + dat$cal13CH4ml[1],
n0 = dat$cal13CH4ml[1],
P = P,
k = result$par["k"],
pool = "CH4")`Optimization of predictions using data quality and deviation in pool size (Nm) and isotopic composition (Nd) over time.
`result <- pdr_optimize(time = dat$time_days,
m = dat$cal12CH4ml + dat$cal13CH4ml,
n = dat$cal13CH4ml,
Nm = dat$Nm,
Nd = dat$Nd,
pool = "CH4",
params = params)`
Optimized rates based on total pool size and atom percent composition.
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.