extract_serie: Extract time series of wrf file list of lat/lon

In eva3dm: Evaluation of 3D Meteorological and Air Quality Models

Extract time series of wrf file list of lat/lon

Description

Read and extract data from a list of wrf output files and a table of lat/lon points based on the distance of the points and the center of model grid points, points outside the domain (and points on domain boundary) are not extracted.

Usage

extract_serie(
  filelist,
  point,
  variable = "o3",
  field = "4d",
  prefix = "serie",
  new = "check",
  return.nearest = FALSE,
  fast = FALSE,
  use_ij = FALSE,
  latitude = "XLAT",
  longitude = "XLONG",
  use_TFLAG = FALSE,
  use_datesec = FALSE,
  id = "id",
  verbose = TRUE
)

Arguments

filelist	list of files to be read
point	data.frame with lat/lon
variable	variable name
field	'4d' (defoult), '3d', '2d' or '2dz' see notes
prefix	to output file, default is serie
new	TRUE, FALSE of 'check' see notes
return.nearest	return the data.frame of nearest points instead of extract the serie
fast	faster calculation of grid distances but less precise
use_ij	logical, use i and j from input instead of calculate
latitude	name of latitude coordinate variable in the netcdf
longitude	name of longitude coordinate variable in the netcdf
use_TFLAG	use the variable TFLAG (CMAQ / smoke) instead of Times (WRF)
use_datesec	use the variable date and datesec (WACCM / CAM-Chem) instead of Times (WRF)
id	name of the column with station names, if point is a SpatVector (points) from terra package
verbose	display additional information

Value

No return value

Note

The field argument '4d' or '2dz' is used to read a 4d/3d variable droping the 3rd dimention (z).

new = TRUE create a new file, new = FALSE append the data in a old file, and new = 'check' check if the file exist and append if the file exist and create if the file doesnt exist

FOR CAMx time-series, use the options: use_TFLAG=T, latitude='latitude', longitude='longitude', new=T

FOR WACCM time-series, use the options: use_datesec=T, latitude='lat', longitude='lon', new=T

The site-list of two global data-sets (METAR and AERONET) are provided on examples and site-list for stations on Brazil (INMET and Air Quality stations).

Examples

cat('Example 1: METAR site list\n')
sites <- readRDS(paste0(system.file("extdata",package="eva3dm"),"/sites_METAR.Rds"))

cat('Example 2: AERONET site list\n')
sites <- readRDS(paste0(system.file("extdata",package="eva3dm"),"/sites_AERONET.Rds"))

cat('Example 3: list of INMET stations on Brazil\n')
sites <- readRDS(paste0(system.file("extdata",package="eva3dm"),"/sites_INMET.Rds"))

cat('Example 4: list of Air Quality stations on Brazil\n')
sites <- readRDS(paste0(system.file("extdata",package="eva3dm"),"/sites_AQ_BR.Rds"))

files    <- dir(path = system.file("extdata",package="eva3dm"),
                pattern = 'wrf.day',
                full.names = TRUE)
dir.create(file.path(tempdir(),"SERIE"))
folder <- file.path(tempdir(),"SERIE")

# extract data for 3 locations
extract_serie(filelist = files, point = sites[1:3,],prefix = paste0(folder,'/serie'))

eva3dm documentation built on April 4, 2025, 5:11 a.m.