Function for plotting wind direction and wind speed profile with the use of wind barbs. Can be launched as standalone function or coupled with pre-drawn Skew-T diagram.

sounding_barbs(
  pressure,
  ws,
  wd,
  altitude,
  ptop = 100,
  interpolate = TRUE,
  showaxis = FALSE,
  barb_cex = 0.3,
  ...
)

Arguments

pressure

pressure [hPa]

ws

wind speed [knots]

wd

wind direction [azimuth in degrees]

altitude

altitude [m] (can be above sea level or above ground level as function always consider first level as surface, i.e h = 0 m) - altitude [m]

ptop

Pressure top level [hPa] to be used for plotting wind speed. Valid options should be < 200 hPa (100 by default)

interpolate

logical, draw wind barbs only at interpolated altitudes with 500 m interval (default = TRUE) instead of all wind barbs for a given input dataset

showaxis

logical, drawing bounding box with left axis for pressure heighs (default FALSE)

barb_cex

size of wind barbs (default = 0.3)

...

extra graphic arguments

Value

wind barbs plot for a given vertical profile of atmosphere

Examples

# load examplary dataset:
data("sounding_vienna")
attach(sounding_vienna)
#> The following objects are masked from sounding_vienna (pos = 3):
#> 
#>     altitude, dpt, pressure, temp, wd, ws
#> The following objects are masked from sounding_vienna (pos = 4):
#> 
#>     altitude, dpt, pressure, temp, wd, ws
#> The following object is masked from package:datasets:
#> 
#>     pressure
sounding_barbs(pressure = pressure, ws = ws, wd = wd, altitude = altitude,
              interpolate = TRUE, showaxis = TRUE)