Vertical interpolation (Wind Components)

Interpolates 3D wind components from model levels to pressure levels.

Usage

Vintp_wind(
  wind_ml,
  ps,
  ap_full,
  b_full,
  ap_half = NULL,
  b_half = NULL,
  pl,
  model_levels = "pl",
  orog = NULL,
  ta_ml = NULL,
  nthreads = 1
)

Arguments

wind_ml

ua or va wind [lon, lat, lev, time] on model levels.

ps

surface pressure [lon, lat, time] (hPa).

ap_full

vertical coordinate formula term ap(k) or a(k) [lev] at layer midpoints (hPa). In case model uses hybrid height coordinates, `ap_full` is hybrid height coordinate (meter).

b_full

vertical coordinate formula term b(k) [lev] at layer midpoints. b(k) is dimensionless fraction. In case model uses hybrid height coordinates, `b_full` is dimensionless fraction.

ap_half

vertical coordinate formula term ap(k) or a(k) [lev+1] at layer interfaces (hPa). It must be provided in case, model uses hybrid height levels.

b_half

vertical coordinate formula term b(k) [lev+1] at layer interfaces. b(k) is dimensionless fraction. It must be provided in case, model uses hybrid height levels.

pl

vector of the required pressure levels [plev] e.g., c(1000,950,850) in hPa.

model_levels

Character. Either "pl" for models on hybrid sigma pressure levels or "ml" for models on hybrid height levels.

orog

surface altitude in meter [lon, lat]. It must be provided in case, model uses hybrid height levels.

ta_ml

temperature on model levels [lon, lat, lev, time]. It must be provided in case, model uses hybrid height levels.

nthreads

an integer specifying the number of threads to use for computation (OpenMP). Default is 1.

Value

wind component [lon, lat, plev, time] on pressure levels.

Examples

data("NorESM2_LM")
ap_full <- NorESM2_LM$a * NorESM2_LM$p0
res <- Vintp_wind(
  wind_ml = NorESM2_LM$ua,
  ps = NorESM2_LM$ps / 100,
  ap_full = ap_full / 100,
  b_full = NorESM2_LM$b,
  pl = seq(1000, 50, -25),
  model_levels = "pl",
  nthreads = 1
)