| sounding_compute {thunder} | R Documentation |
Calculate convective parameters
Description
A core function for calculating convective parameters commonly used in the operational prediction of severe convective storms. Returns a vector of parameters.
Usage
sounding_compute(
pressure,
altitude,
temp,
dpt,
wd,
ws,
accuracy = 2,
interpolate_step = 5,
meanlayer_bottom_top = c(0, 500),
storm_motion = c(999, 999)
)
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 [metres] |
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 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 |
Details
MU_CAPE
MU_CAPE_M10
MU_CAPE_M10_PT
MU_02km_CAPE
MU_03km_CAPE
MU_HGL_CAPE
MU_CIN
MU_LCL_HGT
MU_LFC_HGT
MU_EL_HGT
MU_LI
MU_LI_M10
MU_WMAX
MU_EL_TEMP
MU_LCL_TEMP
MU_LFC_TEMP
MU_MIXR
MU_CAPE_500
MU_CAPE_500_M10
MU_CAPE_500_M10_PT
MU_CIN_500
MU_LI_500
MU_LI_500_M10
SB_CAPE
SB_CAPE_M10
SB_CAPE_M10_PT
SB_02km_CAPE
SB_03km_CAPE
SB_HGL_CAPE
SB_CIN
SB_LCL_HGT
SB_LFC_HGT
SB_EL_HGT
SB_LI
SB_LI_M10
SB_WMAX
SB_EL_TEMP
SB_LCL_TEMP
SB_LFC_TEMP
SB_MIXR
ML_CAPE
ML_CAPE_M10
ML_CAPE_M10_PT
ML_02km_CAPE
ML_03km_CAPE
ML_HGL_CAPE
ML_CIN
ML_LCL_HGT
ML_LFC_HGT
ML_EL_HGT
ML_LI
ML_LI_M10
ML_WMAX
ML_EL_TEMP
ML_LCL_TEMP
ML_LFC_TEMP
ML_MIXR
LR_0500m
LR_01km
LR_02km
LR_03km
LR_04km
LR_06km
LR_16km
LR_26km
LR_24km
LR_36km
LR_26km_MAX
LR_500700hPa
LR_500800hPa
LR_600800hPa
FRZG_HGT
FRZG_wetbulb_HGT
HGT_max_thetae_03km
HGT_min_thetae_04km
Delta_thetae
Delta_thetae_min04km
Thetae_01km
Thetae_02km
DCAPE
Cold_Pool_Strength
Wind_Index
PRCP_WATER
Moisture_Flux_02km
RH_01km
RH_02km
RH_14km
RH_25km
RH_36km
RH_HGL
BS_0500m
BS_01km
BS_02km
BS_03km
BS_06km
BS_08km
BS_36km
BS_26km
BS_16km
BS_18km
BS_EFF_MU
BS_EFF_SB
BS_EFF_ML
BS_SFC_to_M10
BS_1km_to_M10
BS_2km_to_M10
BS_MU_LFC_to_M10
BS_SB_LFC_to_M10
BS_ML_LFC_to_M10
BS_MW02_to_SM
BS_MW02_to_RM
BS_MW02_to_LM
BS_HGL_to_SM
BS_HGL_to_RM
BS_HGL_to_LM
MW_0500m
MW_01km
MW_02km
MW_03km
MW_06km
MW_13km
SRH_100m_RM
SRH_250m_RM
SRH_500m_RM
SRH_1km_RM
SRH_3km_RM
SRH_36km_RM
SRH_100m_LM
SRH_250m_LM
SRH_500m_LM
SRH_1km_LM
SRH_3km_LM
SRH_36km_LM
SV_500m_RM
SV_01km_RM
SV_03km_RM
SV_500m_LM
SV_01km_LM
SV_03km_LM
MW_SR_500m_RM
MW_SR_01km_RM
MW_SR_03km_RM
MW_SR_500m_LM
MW_SR_01km_LM
MW_SR_03km_LM
MW_SR_VM_500m_RM
MW_SR_VM_01km_RM
MW_SR_VM_03km_RM
MW_SR_VM_500m_LM
MW_SR_VM_01km_LM
MW_SR_VM_03km_LM
SV_FRA_500m_RM
SV_FRA_01km_RM
SV_FRA_03km_RM
SV_FRA_500m_LM
SV_FRA_01km_LM
SV_FRA_03km_LM
Bunkers_RM_A
Bunkers_RM_M
Bunkers_LM_A
Bunkers_LM_M
Bunkers_MW_A
Bunkers_MW_M
Corfidi_downwind_A
Corfidi_downwind_M
Corfidi_upwind_A
Corfidi_upwind_M
K_Index
Showalter_Index
TotalTotals_Index
SWEAT_Index
STP_fix
STP_new
STP_fix_LM
STP_new_LM
SCP_fix
SCP_new
SCP_fix_LM
SCP_new_LM
SHIP
HSI
DCP
MU_WMAXSHEAR
SB_WMAXSHEAR
ML_WMAXSHEAR
MU_EFF_WMAXSHEAR
SB_EFF_WMAXSHEAR
ML_EFF_WMAXSHEAR
EHI_500m
EHI_01km
EHI_03km
EHI_500m_LM
EHI_01km_LM
EHI_03km_LM
SHERBS3
SHERBE
SHERBS3_v2
SHERBE_v2
DEI
DEI_eff
TIP
Value
Named vector of 200+ convective indices
Examples
old_options = options(scipen = 99)
pressure = c(1000, 855, 700, 500, 300, 100, 10)
altitude = c(0, 1500, 2500, 6000, 8500, 12000, 25000)
temp = c(25, 10, 0, -15, -30, -50, -92)
dpt = c(20, 5, -5, -30, -55, -80, -99)
wd = c(0, 90, 135, 180, 270, 350, 0)
ws = c(5, 10, 20, 30, 40, 5, 0)
accuracy = 2
sounding_compute(pressure, altitude, temp, dpt, wd, ws, accuracy)
options(old_options)