read.met {oce} | R Documentation |
Read a met File
Description
Reads some meteorological file formats used by the Environment Canada (see
reference 1). Since the agency does not publish the data formats, this
function has to be adjusted every few years, when a user finds that the
format has changed. Caution: as of March 2022, this function fails on some
on Windows machines, for reasons that seem to be related to the handling of
both file encoding and system encoding. Adjusting the encoding
parameter
of this function might help. If not, try reading the data with
read.csv()
and then using as.met()
to create a met object.
Usage
read.met(
file,
type = NULL,
skip = NULL,
encoding = "latin1",
tz = getOption("oceTz"),
debug = getOption("oceDebug")
)
Arguments
file |
a character string naming a file that holds met data. |
type |
if |
skip |
integer giving the number of header lines that precede the
data. This is ignored unless |
encoding |
a character value that indicates the encoding to be used for
this data file, if it is textual. The default value for most functions is
|
tz |
timezone assumed for time data. This defaults to
|
debug |
a flag that turns on debugging. Set to 1 to get a moderate amount of debugging information, or to 2 to get more. |
Value
A met object.
Sample of Usage
# Example 1: "csv1" Environment Canada format (found to be obsolete as of Oct 2019) csv1 <- read.met(system.file("extdata", "test_met_vsn1.csv", package="oce")) csv1 <- oceSetData(csv1, "time", csv1[["time"]]+4*3600, note="add 4h to local time to get UTC time") # Example 2: "csv2" Environment Canada format (found to be obsolete as of Jan 2022) csv2 <- read.met(system.file("extdata", "test_met_vsn2.csv", package="oce")) csv2 <- oceSetData(csv2, "time", csv2[["time"]]+4*3600, note="add 4h to local time to get UTC time") # Example 3: "csv3" Environment Canada format. Note timezone correction csv3 <- read.met(system.file("extdata", "test_met_vsn3.csv", package="oce")) csv3 <- oceSetData(csv3, "time", csv3[["time"]]+4*3600, note="add 4h to local time to get UTC time") # Example 4: "xml2" format. (Uncertain timezone, so not corrected.) if (requireNamespace("XML", quietly=TRUE)) xml2 <- read.met(system.file("extdata", "test_met_xml2.xml", package="oce")) # Example 5: download and plot library(oce) # Recreate data(met) and plot u(t) and v(t) metFile <- download.met(id=6358, year=2003, month=9, destdir=".") met <- read.met(metFile) met <- oceSetData(met, "time", met[["time"]]+4*3600, note="add 4h to local time to get UTC time") plot(met)
Author(s)
Dan Kelley
References
Environment Canada website for Historical Climate Data
https://climate.weather.gc.ca/index_e.html
See Also
Other things related to met data:
[[,met-method
,
[[<-,met-method
,
as.met()
,
download.met()
,
met-class
,
met
,
plot,met-method
,
subset,met-method
,
summary,met-method