xtractor: Extract several variables from WRF outputs

In atmoschem/wrftools: Provide tools to run WRF

Description

Return data.framee or list of raster and df from points. When is data.frame contains cell value, interpolation for lat long based on neighour cells,

Return data.framee or list of raster and df from points. When is data.frame contains cell value, interpolation for lat long based on neighour cells,

Usage

xtractor(
  atmos,
  vars,
  level = 1,
  points,
  fac_res = 1.2,
  stations,
  crs_points = 4326,
  model = "WRF",
  tz = "America/Sao_Paulo",
  return_list = FALSE,
  verbose = TRUE
)

xtractor(
  atmos,
  vars,
  level = 1,
  points,
  fac_res = 1.2,
  stations,
  crs_points = 4326,
  model = "WRF",
  tz = "America/Sao_Paulo",
  return_list = FALSE,
  verbose = TRUE
)

Arguments

atmos	Character; path to output of model (wrfout)
vars	Character; Variable with has 3 or 4 dimensions.
level	Integer; Which level
points	data.frame, matrix, SpatialPointsDataFrame or sf 'POINTS', with coordinates East-weast (lat), north-south (long) and 'Station'.
fac_res	Numeric. Expansion factor for buffer distance. The buffer distane is calculated as the average of the resolution times fac_res
stations	Character; names of stations for each point. (DEPRECATED)
crs_points	Integer, crs points.
model	Character; Currently, only "WRF"
tz	Character; Time zone
return_list	Logical; If TRUE, return a list with raster and data.frames if FALSE, only data.frame.
verbose	logical

Value

Return data.framee or list of raster and df

Return data.framee or list of raster and df

Note

Based on NCL scrip xtractor from DCA/IAG/USP

Based on NCL scrip xtractor from DCA/IAG/USP

Examples

## Not run: 
# do not run
data(cetesb)
cetesb <- cetesb[!is.na(cetesb$Station), "Station"]
#use your wrfout
#wrf <- "~/Documentos/wrfo/wrfoA.nc"
t2 = c("T2", "o3", "no")
df <- xtractor(atmos = wrf, vars = t2, points = cetesb,
stations = cetesb$Station[1:3])
r <- xtractor(atmos = wrf, vars = t2, points = cetesb return_list = T)

## End(Not run)
## Not run: 
# do not run
data(cetesb)
cetesb <- cetesb[!is.na(cetesb$Station), "Station"]
#use your wrfout
#wrf <- "~/Documentos/wrfo/wrfoA.nc"
t2 = c("T2", "o3", "no")
df <- xtractor(atmos = wrf, vars = t2, points = cetesb,
stations = cetesb$Station[1:3])
r <- xtractor(atmos = wrf, vars = t2, points = cetesb return_list = T)

## End(Not run)

atmoschem/wrftools documentation built on Oct. 28, 2021, 4:17 p.m.