fill_mds: Gap-fill a Variable Using the MDS Algorithm

Description Usage Arguments Details Value References

Description

MDS gap filling algorithm adapted after the PV-Wave code and paper by Markus Reichstein. Original name: sEddyProc_sMDSGapFill.

Usage

1
2
3
fill_mds(data, var, qc_var = NULL, qc_val = 0, v1 = "rg", t1 = 50,
  v2 = "vpd", t2 = 5, v3 = "tair", t3 = 2.5, fill_all = TRUE,
  verbose = TRUE, max_run = NULL, only_fill = FALSE)

Arguments

var

A variable to be filled.

qc_var

The quality flag of the variable to be filled

qc_val

The value of quality flag for _good_ (original) data, other data is set to missing

v1

The 1st condition variable (default: Global radiation "Rg" in W m-2)

t1

The tolerance interval for the 1st condition variable (default: 50 W m-2)

v2

The 2nd condition variable (default: Vapor pressure deficit "VPD" in hPa)

t2

The tolerance interval for the 2nd condition variable (default: 5 hPa)

v3

The 3rd condition variable (default: Air temperature "Tair" in degC)

t3

The tolerance interval for the 3rd condition variable (default: 2.5 degC)

fill_all

A logical value whether to fill all values to estimate uncertainties.

verbose

A logical value whether to print status information to screen.

max_run

A scalar integer indicating how many subsequent numerically equal values to allow until a warning is produced. Set to Inf or NA for no warnings. Defaults for "NEE" to records across 4 hours and no warning for others.

only_fill

A logical value whether to only output a vector of the gap- filled variable.

Details

Runs of numerically equal numbers hint to problems of the data and cause unreasonable estimates of uncertainty. The user is warned if this occurs.

MDS gap filling algorithm calls the subroutines Look Up Table fill_lut and Mean Diurnal Course fill_mdc with window sizes as described in the reference. To run dataset only with MDC algorithm fill_mdc, set condition variable v1 to "none".

If meteo condition variables (v1, v2, v3) are at default values but do not exist as columns, they are set to "none" (= disabled identifier). This allows running MDS with less variables than prescribed in the default setting. If meteo condition variable are same as variable to fill, also set to "none". This prevents filling artificial gaps (for uncertainty estimates) with itself as meteo condition variable.

Value

Data frame (or vector if only_fill) with gap filling results.

References

Reichstein, M. et al. (2005) On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Global Change Biology, 11, 1424-1439.


grahamstewart12/tidyflux documentation built on June 4, 2019, 7:44 a.m.