# This wrapper script calls the 'create.thresholds.from.file' function to create thresholds
# based on the details specified below.

library(climdex.pcic.ncdf,lib.loc="/data/local/hadhd/climpact2-master/R/x86_64-redhat-linux-gnu-library/3.3")

# list of one to three input files. e.g. c("a.nc","b.nc","c.nc")
input.files = c("./4.apdaba.pag1j0.05216.nc","./4.apdaba.pag1j0.3236X.nc","./4.apdaba.pag1j0.3236N.nc")

# list of variable names according to above file(s)
vars=c(prec="precipitation_flux", tmax="Tmax", tmin="Tmin")

# output file name
output.file = "precis.thresholds.1961-1990.nc"

# author data
author.data=list(institution="Met Office", institution_id="UKMO")

# reference period
base.range=c(1961,1990)

# number of cores to use (or FALSE)
parallel = 4  #FALSE

# print messages?
verbose=TRUE



######################################
# Do not modify without a good reason.

fclimdex.compatible=FALSE

create.thresholds.from.file(input.files,output.file,author.data,variable.name.map=vars,base.range=base.range,parallel=parallel,verbose=verbose,fclimdex.compatible=fclimdex.compatible)