sEddyProc_sMDSGapFill: sEddyProc sMDSGapFill

In REddyProc: Post Processing of (Half-)Hourly Eddy-Covariance Measurements

sEddyProc sMDSGapFill

Description

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

Usage

sEddyProc_sMDSGapFill(Var = Var.s, QFVar = if (!missing(QFVar.s)) QFVar.s else "none", 
    QFValue = if (!missing(QFValue.n)) QFValue.n else NA_real_, 
    V1 = if (!missing(V1.s)) V1.s else "Rg", 
    T1 = if (!missing(T1.n)) T1.n else 50, 
    V2 = if (!missing(V2.s)) V2.s else "VPD", 
    T2 = if (!missing(T2.n)) T2.n else 5, 
    V3 = if (!missing(V3.s)) V3.s else "Tair", 
    T3 = if (!missing(T3.n)) T3.n else 2.5, 
    FillAll = if (!missing(FillAll.b)) FillAll.b else TRUE, 
    isVerbose = if (!missing(Verbose.b)) Verbose.b else TRUE, 
    suffix = if (!missing(Suffix.s)) Suffix.s else "", 
    minNWarnRunLength = if (Var == "NEE") 4 * 
        .self$sINFO$DTS/24 else NA_integer_, 
    Var.s, QFVar.s, QFValue.n, V1.s, T1.n, 
    V2.s, T2.n, V3.s, T3.n, FillAll.b, Verbose.b, 
    Suffix.s)

Arguments

Var	Variable to be filled
QFVar
QFValue
V1	Condition variable 1 (default: Global radiation 'Rg' in W m-2)
T1	Tolerance interval 1 (default: 50 W m-2)
V2	Condition variable 2 (default: Vapour pressure deficit 'VPD' in hPa)
T2	Tolerance interval 2 (default: 5 hPa)
V3	Condition variable 3 (default: Air temperature 'Tair' in degC)
T3	Tolerance interval 3 (default: 2.5 degC)
FillAll	Fill all values to estimate uncertainties
isVerbose	Print status information to screen
suffix	String suffix needed for different processing setups on the same dataset (for explanations see below)
minNWarnRunLength	scalar integer: warn if number of subsequent numerically equal values exceeds this number. Set to Inf or NA for no warnings. defaults for "NEE" to records across 4 hours and no warning for others.
Var.s	deprecated
QFVar.s	deprecated
QFValue.n	deprecated
V1.s	deprecated
T1.n	deprecated
V2.s	deprecated
T2.n	deprecated
V3.s	deprecated
T3.n	deprecated
FillAll.b	deprecated
Verbose.b	deprecated
Suffix.s	deprecated

Details

Initialize temporal data frame sTEMP for newly generated gap filled data and qualifiers, see sEddyProc_sFillInit for explanations on suffixes.

Runs of numerically equal numbers hint to problems of the data and cause unreasonable estimates of uncertainty. This routine warns the user.

MDS gap filling algorithm calls the subroutines Look Up Table sEddyProc_sFillLUT and Mean Diurnal Course sEddyProc_sFillMDC with different window sizes as described in the reference.

To run dataset only with MDC algorithm sEddyProc_sFillMDC, set condition variable V1 to 'none'.

Different processing setups on the same dataset

Attention: When processing the same site data set with different setups for the gap filling or flux partitioning (e.g. due to different ustar filters), a string suffix is needed! This suffix is added to the result column names to distinguish the results of the different setups.

Value

Gap filling results in sTEMP data frame (with renamed columns).

Author(s)

AMM, TW Department for Biogeochemical Integration at MPI-BGC, Jena, Germany <[e-mail]> [cph], Thomas Wutzler <[e-mail]> [aut, cre], Markus Reichstein <[e-mail]> [aut], Antje Maria Moffat <[e-mail]> [aut, trl], Olaf Menzer <[e-mail]> [ctb], Mirco Migliavacca <[e-mail]> [aut], Kerstin Sickel <[e-mail]> [ctb, trl], Ladislav <U+0160>igut <[e-mail]> [ctb]

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.

REddyProc documentation built on March 18, 2022, 5:41 p.m.