selepe: QC for Atmospheric Pressure (PP)
In INQC: Quality Control of Climatological Daily Time Series

Description Usage Arguments Value Examples

This function will centralize temperature-like QC routines. It will create a file in the folder QC with an additional 0/1 column where "1" means test failed.

selepe(
  element = "PP",
  large = 15000,
  small = 8000,
  maxjump = 2000,
  maxseq = 3,
  margina = 0.999,
  level = 5,
  window = 30,
  roundmax = 10,
  blockmanymonth = 15,
  blockmanyyear = 180,
  blocksizeround = 20,
  qjump = 0.999,
  tjump = 1.5,
  inisia = FALSE
)

`element`	two-letters ECA&D code for the element (PP for sea level pressure)
`large`	value above which the observation is considered physically impossible for the region
`small`	value below which the observation is considered physically impossible for the region
`maxjump`	forcing for jumps2() in absolute mode (in the same units of the variable). Passed on to jumps2(). See ?jumps2 for further details.
`maxseq`	maximum number of consecutive repeated values, for flat function (11.1,11.1,11.1 would be 3 consecutive values)
`margina`	tolerance margin, expressed as quantile of the differences, FUNCTION: newfriki(). Passed on to newfriki(). See ?newfriki for details
`level`	number of IQRs for IQR outliers
`window`	window, in days, for IQR outliers
`roundmax`	maximum number of consecutive decimal part value, for flat function (10.0, 11.0, 12.0 would be 3 consecutive values)
`blockmanymonth`	maximum number of equal values in a month, FUNCTION: toomany()
`blockmanyyear`	maximum number of equal values in a year, FUNCTION: toomany()
`blocksizeround`	maximum number of values in a month with the same decimal, for rounding function
`qjump`	quantile for jumps2() in quantile mode. Passed on to jumps2(). See ?jumps2 for further details
`tjump`	factor to multiply the quantile value for jumps2(). Passed on to jumps2(). See ?jumps2 for further details
`inisia`	a logical flag. If it is TRUE inithome() will be called

results of QC for PP

#Set a temporal working directory:
wd <- tempdir()
wd0 <- setwd(wd)
#Create subdirectory where raw data files have to be located
dir.create(file.path(wd, 'raw'))
options("homefolder"='./'); options("blend"=FALSE)
#Extract the ECA&D data and station files from the example data folder
path2pplist<-system.file("extdata", "ECA_blend_source_pp.txt", package = "INQC")
pplist<-readr::read_lines_raw(path2pplist)
readr::write_lines(pplist,'ECA_blend_source_pp.txt')
path2ppdata<-system.file("extdata", "PP_SOUID132729.txt", package = "INQC")
ppdata<-readr::read_lines_raw(path2ppdata)
readr::write_lines(ppdata, file=paste(wd,'/raw/PP_SOUID132729.txt',sep=''))
#Perform QC of Atmospheric Pressure data
selepe(inisia=TRUE)
#Remove some temporary files
list = list.files(pattern = "Rfwf")
file.remove(list)
#Return to user's working directory:
setwd(wd0)
#The QC results can be found in the directory:
print(wd)