read.functions: Read data as provided by Arpa-ER

Description Usage Arguments Note

Description

Functions to read some data as provided by Arpa Emilia-Romagna

Usage

1
2
3
4
5
6
7
8
9
read.qaria (file) 
read.sql   (file) 
qaria2long (datafiles, anafile, codes=NULL)

read.ncdf.arpaer(con=NULL, pollutant="pm10", lev=1, 
                 tz.in="UTC", tz.out="Africa/Algiers")
read.grads(filectl) 

read.field(file, coords.col=1:2, data.col=3, coords.fact=1, ...)

Arguments

file

name of the input file

datafiles

vector of names of files (format 'estra_qaria', an internal standard for some people in Arpa-ER)

anafile

file with stations metadata

codes

vector with stations codes, if NULL it is argued from the elements of datafiles, if they are in the form /path/to/files/pollutant_stationcode.asc

con

name of a NetCDF file, or an already open connection to NetCDF

pollutant

pollutant code used in the NetCDF

lev

level

tz.in

timezone used by the CTM

tz.out

timezone used for the observed data

filectl

name of the GrADS ctl file

coords.col

a vector with the numbers of the columns containing the coordinates

data.col

a scalar with the number of the column containing the data

coords.fact

multiplication factor for coordinates (e.g. 1000 if you need to convert kilometers to meters)

...

optional arguments; will be passed to read.table, called by read.field

Note

If your input data are not provided following internal 'standard' formats of Arpa-ER, you don't need these functions, except maybe read.field, which is quite general. Then, you'll need to write your own reading functions, suitable for your data format.

The function read.grads is deprecated, since I wrote it when I was young, and it needs external executables.


jobonaf/pesco documentation built on May 19, 2019, 2:58 p.m.