knitr::opts_chunk$set( collapse = TRUE, comment = "#>", fig.path = "README-" ) library(glimmr) library(dplyr)
Gas Fluxes and Dynamic Chamber Measurements
The glimmr
package provides functions to convert high frequency concentration data
obtained by (dynamic-) chamber measurements into gasfluxes. Data recorded with different devices can be processed automatically. The data is separated into chunks covering the single chamber applications and containing additional meta information (e.g. IDs, temperature...) provided by a metadata file.
Furthermore glimmr
contains functions to calculate gasfluxes with the 'Boundary Layer Equation' (BLE).
To calculate gasfluxes from high frequency chamber measurements, glimmr offers a cusomizable system of gas analyzer definitions (analyzer()
), as well as two preconfigured devices with read_
, process_
and inspect_
functions:
process_chamber()
& inspect_chamber()
: custom devicesread_gasmet()
, process_gasmet()
& inspect_gasmet()
: GASMETread_losgatos()
, process_losgatos()
& inspect_losgatos()
: LosGatos# The package is not released on CRAN yet. # Development version from GitHub: # install.packages("devtools") devtools::install_github("tekknosol/glimmr")
Calculate fluxes from device records and a meta file containing information about single measurements.
Fluxes are reported in [mmol/m²/d] using linear model fitting and robust linear model fitting with robust::lmRob()
# read data gasmet <- read_gasmet("path/to/gasmet_file.txt") meta_gasmet <- read_csv("path/to/meta_file.csv") # create diagnostic plot(s) inspect_gasmet(gasmet, meta_gasmet)
# compute fluxes in mmol/m²/d # F_LM: Liner Model # F_RLM: Robust Liner Model process_gasmet(gasmet, meta_gasmet)
Calculate gas transfer coefficient used with the Boundary Layer Equation (BLE).
# Windpseed in m/s windspeed <- c(2, 4.3, 1.8) # Water temperature in °C watertemperature <- c(18, 21, 19.3) # Calculate gas transfer coefficient calc_kW(windspeed, watertemperature)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.