# 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)