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

library(climdex.pcic.ncdf)
# list of one to three input files. e.g. c("a.nc","b.nc","c.nc")
input.files = c("./4.rabtfa.pag0i9.05216.nc","./4.rabtfa.pag0i9.3236X.nc","./4.rabtfa.pag0i9.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.1960-1989.nc"

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

# reference period
base.range=c(1960,1989)

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

# print messages?
verbose=TRUE



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

fclimdex.compatible=TRUE

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)