Calculate convective parameters

A core function for calculating convective parameters commonly used in the operational prediction of severe convective storms. Returns a vector of parameters.

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

Value

Named vector of 200+ convective indices

Details

1. MU_CAPE

2. MU_CAPE_M10

3. MU_CAPE_M10_PT

4. MU_02km_CAPE

5. MU_03km_CAPE

6. MU_HGL_CAPE

7. MU_CIN

8. MU_LCL_HGT

9. MU_LFC_HGT

10. MU_EL_HGT

11. MU_LI

12. MU_LI_M10

13. MU_WMAX

14. MU_EL_TEMP

15. MU_LCL_TEMP

16. MU_LFC_TEMP

17. MU_MIXR

18. MU_CAPE_500

19. MU_CAPE_500_M10

20. MU_CAPE_500_M10_PT

21. MU_CIN_500

22. MU_LI_500

23. MU_LI_500_M10

24. SB_CAPE

25. SB_CAPE_M10

26. SB_CAPE_M10_PT

27. SB_02km_CAPE

28. SB_03km_CAPE

29. SB_HGL_CAPE

30. SB_CIN

31. SB_LCL_HGT

32. SB_LFC_HGT

33. SB_EL_HGT

34. SB_LI

35. SB_LI_M10

36. SB_WMAX

37. SB_EL_TEMP

38. SB_LCL_TEMP

39. SB_LFC_TEMP

40. SB_MIXR

41. ML_CAPE

42. ML_CAPE_M10

43. ML_CAPE_M10_PT

44. ML_02km_CAPE

45. ML_03km_CAPE

46. ML_HGL_CAPE

47. ML_CIN

48. ML_LCL_HGT

49. ML_LFC_HGT

50. ML_EL_HGT

51. ML_LI

52. ML_LI_M10

53. ML_WMAX

54. ML_EL_TEMP

55. ML_LCL_TEMP

56. ML_LFC_TEMP

57. ML_MIXR

58. LR_0500m

59. LR_01km

60. LR_02km

61. LR_03km

62. LR_04km

63. LR_06km

64. LR_16km

65. LR_26km

66. LR_24km

67. LR_36km

68. LR_26km_MAX

69. LR_500700hPa

70. LR_500800hPa

71. LR_600800hPa

72. FRZG_HGT

73. FRZG_wetbulb_HGT

74. HGT_max_thetae_03km

75. HGT_min_thetae_04km

76. Delta_thetae

77. Delta_thetae_min04km

78. Thetae_01km

79. Thetae_02km

80. DCAPE

81. Cold_Pool_Strength

82. Wind_Index

83. PRCP_WATER

84. Moisture_Flux_02km

85. RH_01km

86. RH_02km

87. RH_14km

88. RH_25km

89. RH_36km

90. RH_HGL

91. BS_0500m

92. BS_01km

93. BS_02km

94. BS_03km

95. BS_06km

96. BS_08km

97. BS_36km

98. BS_26km

99. BS_16km

100. BS_18km

101. BS_EFF_MU

102. BS_EFF_SB

103. BS_EFF_ML

104. BS_SFC_to_M10

105. BS_1km_to_M10

106. BS_2km_to_M10

107. BS_MU_LFC_to_M10

108. BS_SB_LFC_to_M10

109. BS_ML_LFC_to_M10

110. BS_MW02_to_SM

111. BS_MW02_to_RM

112. BS_MW02_to_LM

113. BS_HGL_to_SM

114. BS_HGL_to_RM

115. BS_HGL_to_LM

116. MW_0500m

117. MW_01km

118. MW_02km

119. MW_03km

120. MW_06km

121. MW_13km

122. SRH_100m_RM

123. SRH_250m_RM

124. SRH_500m_RM

125. SRH_1km_RM

126. SRH_3km_RM

127. SRH_36km_RM

128. SRH_100m_LM

129. SRH_250m_LM

130. SRH_500m_LM

131. SRH_1km_LM

132. SRH_3km_LM

133. SRH_36km_LM

134. SV_500m_RM

135. SV_01km_RM

136. SV_03km_RM

137. SV_500m_LM

138. SV_01km_LM

139. SV_03km_LM

140. MW_SR_500m_RM

141. MW_SR_01km_RM

142. MW_SR_03km_RM

143. MW_SR_500m_LM

144. MW_SR_01km_LM

145. MW_SR_03km_LM

146. MW_SR_VM_500m_RM

147. MW_SR_VM_01km_RM

148. MW_SR_VM_03km_RM

149. MW_SR_VM_500m_LM

150. MW_SR_VM_01km_LM

151. MW_SR_VM_03km_LM

152. SV_FRA_500m_RM

153. SV_FRA_01km_RM

154. SV_FRA_03km_RM

155. SV_FRA_500m_LM

156. SV_FRA_01km_LM

157. SV_FRA_03km_LM

158. Bunkers_RM_A

159. Bunkers_RM_M

160. Bunkers_LM_A

161. Bunkers_LM_M

162. Bunkers_MW_A

163. Bunkers_MW_M

164. Corfidi_downwind_A

165. Corfidi_downwind_M

166. Corfidi_upwind_A

167. Corfidi_upwind_M

168. K_Index

169. Showalter_Index

170. TotalTotals_Index

171. SWEAT_Index

172. STP_fix

173. STP_new

174. STP_fix_LM

175. STP_new_LM

176. SCP_fix

177. SCP_new

178. SCP_fix_LM

179. SCP_new_LM

180. SHIP

181. HSI

182. DCP

183. MU_WMAXSHEAR

184. SB_WMAXSHEAR

185. ML_WMAXSHEAR

186. MU_EFF_WMAXSHEAR

187. SB_EFF_WMAXSHEAR

188. ML_EFF_WMAXSHEAR

189. EHI_500m

190. EHI_01km

191. EHI_03km

192. EHI_500m_LM

193. EHI_01km_LM

194. EHI_03km_LM

195. SHERBS3

196. SHERBE

197. SHERBS3_v2

198. SHERBE_v2

199. DEI

200. DEI_eff

201. TIP

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)
#>              MU_CAPE          MU_CAPE_M10       MU_CAPE_M10_PT 
#>       2269.925733301        998.144308627        313.074714966 
#>         MU_02km_CAPE         MU_03km_CAPE          MU_HGL_CAPE 
#>        247.979380820        575.629305978       1616.538411895 
#>               MU_CIN           MU_LCL_HGT           MU_LFC_HGT 
#>          0.000000000        730.000000000        730.000000000 
#>            MU_EL_HGT                MU_LI            MU_LI_M10 
#>       8300.000000000        -10.111930240        -10.853947310 
#>              MU_WMAX           MU_EL_TEMP          MU_LCL_TEMP 
#>         67.378419888        -28.800000000         17.700000000 
#>          MU_LFC_TEMP              MU_MIXR          MU_CAPE_500 
#>         17.700000000         14.875917910       1076.032178138 
#>      MU_CAPE_500_M10   MU_CAPE_500_M10_PT           MU_CIN_500 
#>        389.365140754        137.081449189          0.000000000 
#>            MU_LI_500        MU_LI_500_M10              SB_CAPE 
#>         -5.041702392         -6.234631722       2269.925733301 
#>          SB_CAPE_M10       SB_CAPE_M10_PT         SB_02km_CAPE 
#>        998.144308627        313.074714966        247.979380820 
#>         SB_03km_CAPE          SB_HGL_CAPE               SB_CIN 
#>        575.629305978       1616.538411895          0.000000000 
#>           SB_LCL_HGT           SB_LFC_HGT            SB_EL_HGT 
#>        730.000000000        730.000000000       8300.000000000 
#>                SB_LI            SB_LI_M10              SB_WMAX 
#>        -10.111930240        -10.853947310         67.378419888 
#>           SB_EL_TEMP          SB_LCL_TEMP          SB_LFC_TEMP 
#>        -28.800000000         17.700000000         17.700000000 
#>              SB_MIXR              ML_CAPE          ML_CAPE_M10 
#>         14.875917910       1646.063927273        670.100095437 
#>       ML_CAPE_M10_PT         ML_02km_CAPE         ML_03km_CAPE 
#>        225.281641180        164.079828657        422.429000173 
#>          ML_HGL_CAPE               ML_CIN           ML_LCL_HGT 
#>       1250.022127396          0.000000000        975.000000000 
#>           ML_LFC_HGT            ML_EL_HGT                ML_LI 
#>        975.000000000       7900.000000000         -7.620329500 
#>            ML_LI_M10              ML_WMAX           ML_EL_TEMP 
#>         -8.584526572         57.377067323        -26.400000000 
#>          ML_LCL_TEMP          ML_LFC_TEMP              ML_MIXR 
#>         15.250000000         15.250000000         13.048747763 
#>             LR_0500m              LR_01km              LR_02km 
#>        -10.000000000        -10.000000000        -10.000000000 
#>              LR_03km              LR_04km              LR_06km 
#>         -9.047619048         -7.857142857         -6.666666667 
#>              LR_16km              LR_26km              LR_24km 
#>         -6.000000000         -5.000000000         -5.767195767 
#>              LR_36km          LR_26km_MAX         LR_500700hPa 
#>         -4.285714286         -5.714285714         -4.285714286 
#>         LR_500800hPa         LR_600800hPa             FRZG_HGT 
#>         -5.180722892         -5.833333333       2500.000000000 
#>     FRZG_wetbulb_HGT  HGT_max_thetae_03km  HGT_min_thetae_04km 
#>       2275.000000000          0.000000000       3750.000000000 
#>         Delta_thetae Delta_thetae_min04km          Thetae_01km 
#>         28.069785504         28.834554624        330.532272184 
#>          Thetae_02km                DCAPE   Cold_Pool_Strength 
#>        323.619051941        598.310049405         12.632184321 
#>           Wind_Index           PRCP_WATER   Moisture_Flux_02km 
#>         33.906407176         27.104631918         30.425480722 
#>              RH_01km              RH_02km              RH_14km 
#>          0.729103473          0.719706866          0.645204350 
#>              RH_25km              RH_36km               RH_HGL 
#>          0.555017886          0.443603638          0.460317694 
#>             BS_0500m              BS_01km              BS_02km 
#>          1.917219591          3.834439182          8.782108409 
#>              BS_03km              BS_06km              BS_08km 
#>         12.655968060         18.005540000         17.407678201 
#>              BS_36km              BS_26km              BS_16km 
#>          9.369289559         13.330420492         16.647826116 
#>              BS_18km            BS_EFF_MU            BS_EFF_SB 
#>         20.279072460         14.223187917         14.223187917 
#>            BS_EFF_ML        BS_SFC_to_M10        BS_1km_to_M10 
#>         13.896750821         15.510376997         13.649859352 
#>        BS_2km_to_M10     BS_MU_LFC_to_M10     BS_SB_LFC_to_M10 
#>          9.883032095         14.073701016         14.073701016 
#>     BS_ML_LFC_to_M10        BS_MW02_to_SM        BS_MW02_to_RM 
#>         13.686365658          7.304027597         10.140964728 
#>        BS_MW02_to_LM         BS_HGL_to_SM         BS_HGL_to_RM 
#>         10.786958453          4.893411960          7.785950290 
#>         BS_HGL_to_LM             MW_0500m              MW_01km 
#>          9.988490342          2.308638103          2.425112273 
#>              MW_02km              MW_03km              MW_06km 
#>          3.347646594          4.800331697          7.810700539 
#>              MW_13km          SRH_100m_RM          SRH_250m_RM 
#>          6.838886565          4.253473515         10.053664673 
#>          SRH_500m_RM           SRH_1km_RM           SRH_3km_RM 
#>         19.720649935         39.634639574        152.521946573 
#>          SRH_36km_RM          SRH_100m_LM          SRH_250m_LM 
#>        236.590130791          1.502688386          3.551808913 
#>          SRH_500m_LM           SRH_1km_LM           SRH_3km_LM 
#>          6.967009791         14.002323601        -13.130823085 
#>          SRH_36km_LM           SV_500m_RM           SV_01km_RM 
#>        -24.378976558          0.003904164          0.003879744 
#>           SV_03km_RM           SV_500m_LM           SV_01km_LM 
#>          0.004796086          0.001013916          0.001100521 
#>           SV_03km_LM        MW_SR_500m_RM        MW_SR_01km_RM 
#>         -0.001435868         10.086252358         10.150075094 
#>        MW_SR_03km_RM        MW_SR_500m_LM        MW_SR_01km_LM 
#>          9.535939066         13.782126361         12.858486067 
#>        MW_SR_03km_LM     MW_SR_VM_500m_RM     MW_SR_VM_01km_RM 
#>          8.457941382         10.107849533         10.225339777 
#>     MW_SR_VM_03km_RM     MW_SR_VM_500m_LM     MW_SR_VM_01km_LM 
#>         10.535820052         13.764654240         12.837090856 
#>     MW_SR_VM_03km_LM       SV_FRA_500m_RM       SV_FRA_01km_RM 
#>          8.834244774          0.998219354          0.987136746 
#>       SV_FRA_03km_RM       SV_FRA_500m_LM       SV_FRA_01km_LM 
#>          0.956008752          0.259238846          0.280009490 
#>       SV_FRA_03km_LM         Bunkers_RM_A         Bunkers_RM_M 
#>         -0.286213018        209.404566875          7.793284819 
#>         Bunkers_LM_A         Bunkers_LM_M         Bunkers_MW_A 
#>        122.058537262         13.182518634        151.949436108 
#>         Bunkers_MW_M   Corfidi_downwind_A   Corfidi_downwind_M 
#>          7.810700539        218.695496108         14.698235841 
#>     Corfidi_upwind_A     Corfidi_upwind_M              K_Index 
#>        231.328321825          9.179441145         24.354838710 
#>      Showalter_Index    TotalTotals_Index          SWEAT_Index 
#>          3.750147453         44.354838710        106.416793661 
#>              STP_fix              STP_new           STP_fix_LM 
#>          0.359981513          0.200492750          0.127176069 
#>           STP_new_LM              SCP_fix              SCP_new 
#>          0.070831081          6.233760868          4.924259548 
#>           SCP_fix_LM           SCP_new_LM                 SHIP 
#>         -0.536673003         -0.423936243          0.628661858 
#>                  HSI                  DCP         MU_WMAXSHEAR 
#>          1.715903102          1.150670745       1213.184834430 
#>         SB_WMAXSHEAR         ML_WMAXSHEAR     MU_EFF_WMAXSHEAR 
#>       1213.184834430       1033.105080761        958.335927637 
#>     SB_EFF_WMAXSHEAR     ML_EFF_WMAXSHEAR             EHI_500m 
#>        958.335927637        797.354807442          0.279777567 
#>             EHI_01km             EHI_03km          EHI_500m_LM 
#>          0.562298052          2.163834321          0.098841218 
#>          EHI_01km_LM          EHI_03km_LM              SHERBS3 
#>          0.198651467         -0.186287458          0.648168875 
#>               SHERBE           SHERBS3_v2            SHERBE_v2 
#>          0.701454223          0.864225167          0.935272297 
#>                  DEI              DEI_eff                  TIP 
#>          1.519761039          1.188457460          2.435630415 
options(old_options)