thunder (version 1.1.4)

sounding_export: Sounding export

Description

Internal package function for exporting interpolated profile with 5 m (or user-defined) steps

Usage

sounding_export(
  pressure,
  altitude,
  temp,
  dpt,
  wd,
  ws,
  accuracy = 3,
  interpolate_step = 5,
  meanlayer_bottom_top = c(0, 500),
  storm_motion = c(999, 999)
)

Value

Data frame of computed values for visualizing parcel trajectories

1. pressure pressure [hPa]

2. altitude altitude [m]

3. temp temperature [degree Celsius]

4. tempV virtual temperature [degree Celsius]

5. dpt dew point temperature [degree Celsius]

6. wd wind direction [azimuth in degrees]

7. ws wind speed [knots]

8. MU temperature for most unstable CAPE trajectory [degree Celsius]

9. SB temperature for surface based CAPE trajectory [degree Celsius]

10. ML temperature for mixed layer CAPE trajectory [degree Celsius]

Arguments

pressure

pressure [hPa]

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 [meters]

temp

temperature [degree Celsius]

dpt

dew point temperature [degree Celsius]

wd

wind direction [azimuth in degrees]

ws

wind speed [knots]

accuracy

accuracy of computations where 3 = high (slow), 2 = medium (recommended), 1 = low (fast)

interpolate_step

interpolation step to be used for vertical interpolation. Valid only if `accuracy` is set to 3 (default is set to 5 m)

meanlayer_bottom_top

(optional) vector of length 2 for bottom and top heights used for computing parcel starting parameters; default: 0, 500

storm_motion

(optional) for moving storms only - one can define vector of length two with wind speed [m/s] and wind directions [degrees] that will be used to compute adjusted SRH parameters

Examples

data("sounding_vienna")
attach(sounding_vienna)
skewt_plot(close_par = FALSE)
output = sounding_export(pressure, altitude, temp, dpt, wd, ws)
skewt_lines(output$dpt, output$pressure, col = "forestgreen", lwd = 2.5)
skewt_lines(output$temp, output$pressure, col = "red", lwd = 2.5)
skewt_lines(output$MU, output$pressure, col = "orange", lty = 1, lwd = 2)
skewt_lines(output$tempV, output$pressure, col = "red3", lty = 3, lwd = 1.5)