【MDS】

path="D:/Fluxnet/OUTCOME/MDS_TRY"
outpath="C:/Users/Lenovo/Desktop/R/outcome_try"
getwd()
setwd("D:/Fluxnet/OUTCOME/MDS_TRY")
#+++ load libraries used in this vignette
library(REddyProc)
library(dplyr)



file <- dir(path)
for(i in 1:length(file)){
  
  EddyData <- fLoadTXTIntoDataframe(file[i]) 
  
  #+++ Replace long runs of equal NEE values by NA
  EddyData <- filterLongRuns(EddyData, "LE")
  
  #+++ Add time stamp in POSIX time format and filter long runs of equal values 
  EddyDataWithPosix <- fConvertTimeToPosix(
    EddyData, 'YDH',Year = 'Year',Day = 'DoY', Hour = 'Hour') 
  
  EProcMDS <- sEddyProc$new(
    file[i], EddyDataWithPosix, c('LE', 'Rg', 'Tair', 'VPD'))
    
  #Initialize 'NEE' as variable to fill
  EProcMDS$sFillInit('LE')
  # Set variables and tolerance