R/wrf_get.R
Read a variable
wrf_get(file = file.choose(), name = "E_CO", as_raster = FALSE, raster_crs = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
file | name of file interactively (default) or specified |
---|---|
name | name of the variable (any variable) |
as_raster | return a raster instead of an array |
raster_crs | crs to use if as_raster is TRUE |
array or raster object
{ # create the folder and emission file dir.create(file.path(tempdir(), "EMISS")) wrf_create(wrfinput_dir = system.file("extdata", package = "eixport"), wrfchemi_dir = file.path(tempdir(), "EMISS")) # get the name of created file files <- list.files(path = file.path(tempdir(), "EMISS"), pattern = "wrfchemi", full.names = TRUE) # open, put some numbers and write CO <- wrf_get(file = files[1], name = "E_CO") CO[] = rnorm(length(CO)) wrf_put(file = files[1], name = "E_CO", POL = CO) COr <- wrf_get(file = files[1], name = "E_CO", as_raster = TRUE) }