rean2climatol: Append a reanalysis series in CSV to 'climatol' input files.

In climatol: Climate Tools (Series Homogenization and Derived Products)

Append a reanalysis series in CSV to climatol input files.

Description

This function reads a reanalysis series in CSV format and appends their corresponging data to the climatol input file *.dat. A line is also appended to the stations file *.est.

Usage

  rean2climatol(varcli, anyi, anyf, csvfile, X, Y, skip=0, datacol=1:2,
  dateformat='%Y-%m-%d', a=0, b=1, ndec=1, rname='rean', sep=',', dec='.')

Arguments

varcli	Short name of the climatic variable under study.
anyi, anyf	First and last years to study.
csvfile	Name of the CSV (text) file containing the reanalysis series.
X, Y	Longitude and Latitude of the series in degrees.
skip	Number of header lines to be skipped. [0]
datacol	Columns holding dates and data. [1:2]
dateformat	Format of dates if not in separate columns. ['%Y-%m-%d']
a, b	Parameters of the optional transformation a+b*dat to be applied to the reanalyis series to adapt them to the units used in our data files. [0, 1]
ndec	No. of decimals to round. [1]
rname	Name of the reanalysis. ['rean']
sep	Data separator (',' by default: Comma Separated Values).
dec	Decimal point ('.' by default).

Details

If 3 (4) values are provided in the datacol parameter, dates are expected to appear as year, month (and day) in separate columns. Otherwise, dates will be provided as character strings with dateformat format.

The reanalysis series will be assigned an elevation of 99 m. (The user can edit the *.est file to change it.)

See Also

homogen

Examples

## Set a temporal working directory and write input files:
wd <- tempdir()
wd0 <- setwd(wd)

## Save example climatol input files:
load('~/cc/climatol/data/climatol_data.rda')
dat <- as.matrix(RR3st[,2:4])
write(dat,'Prec_1981-1995.dat')
write.table(SIstations,'Prec_1981-1995.est',row.names=FALSE,col.names=FALSE)

## Copy the provided simulated reanalysis series:
file.copy(exampleFiles('bogusERA5.txt'),'.')

## Now run the function. (Reanalysis precipitation rates in kg/m2/s
##  are converted to mm with parameter b):
rean2climatol('Prec', 1981, 1995, 'bogusERA5.txt', 15, 46.2, datacol=c(2:4,1),
  rname='bogusERA5', b=86400)

## Return to user's working directory:
setwd(wd0)

## New output files can be found in directory:
print(wd)

climatol documentation built on May 20, 2026, 9:08 a.m.