ThundeR
thundeR - Rapid computation and visualisation of convective parameters from rawinsonde and NWP data
Install / Use
/learn @bczernecki/ThundeRREADME
thundeR <img src="https://github.com/bczernecki/thundeR/raw/master/man/figures/logo.png" align="right" width="150" />
<!-- badges: start -->
Rapid computation and visualisation of convective parameters from rawinsonde and numerical weather prediction data
thundeR is a freeware R package for
rapid computation and visualisation of convective parameters commonly
used in the operational forecasting of severe convective storms. Core
algorithm is based on C++ code implemented into R language via Rcpp.
This solution allows to compute over 200 thermodynamic and kinematic
parameters in less than 0.02s per profile and process large datasets
such as reanalyses or operational NWP models in a reasonable amount of
time. Package has been developed since 2017 by research meteorologists
specializing in severe convective storms and is constantly updated with
new features.
Online browser
Online rawinsonde browser of thundeR package is available at http://rawinsonde.com
Installation
The stable version can be installed from the CRAN repository:
install.packages("thunder")
The development version can be installed directly from the github repository:
remotes::install_github("bczernecki/thunder")
Usage
Draw Skew-T, hodograph and convective parameters on a single layout and export to png file
data("sounding_vienna") # load example dataset (Vienna rawinsonde profile for 23 Aug 2011 12UTC):
pressure = sounding_vienna$pressure # vector of pressure [hPa]
altitude = sounding_vienna$altitude # vector of altitude [meters]
temp = sounding_vienna$temp # vector of temperature [degree Celsius]
dpt = sounding_vienna$dpt # vector of dew point temperature [degree Celsius]
wd = sounding_vienna$wd # vector of wind direction [azimuth in degrees]
ws = sounding_vienna$ws # vector of wind speed [knots]
sounding_save(filename = "Vienna.png", title = "Vienna - 23 August 2011 1200 UTC", pressure, altitude, temp, dpt, wd, ws)
Download LBF North Platte rawinsonde profile for 03 Jul 1999 00UTC and export to png file
profile = get_sounding(wmo_id = 72562, yy = 1999, mm = 7, dd = 3, hh = 0)
sounding_save(filename = "NorthPlatte.png", title = "North Platte - 03 July 1999 0000 UTC", profile$pressure, profile$altitude, profile$temp, profile$dpt, profile$wd, profile$ws)
Compute convective parameters based on a sample vertical profile data:
library("thunder")
pressure = c(1000, 855, 700, 500, 300, 100, 10) # pressure [hPa]
altitude = c(0, 1500, 2500, 6000, 8500, 12000, 25000) # altitude [meters]
temp = c(25, 10, 0, -15, -30, -50, -92) # temperature [degree Celsius]
dpt = c(20, 5, -5, -30, -55, -80, -99) # dew point temperature [degree Celsius]
wd = c(0, 90, 135, 180, 270, 350, 0) # wind direction [azimuth in degress]
ws = c(5, 10, 20, 30, 40, 5, 0) # wind speed [knots]
accuracy = 2 # accuracy of computations where 3 = high (slow), 2 = medium (recommended), 1 = low (fast)
sounding_compute(pressure, altitude, temp, dpt, wd, ws, accuracy)
# MU_CAPE MU_CAPE_M10 MU_CAPE_M10_PT MU_02km_CAPE
# 2269.9257 998.1443 313.0747 247.9794
# MU_03km_CAPE MU_HGL_CAPE MU_CIN MU_LCL_HGT
# 575.6293 1616.5384 0.0000 730.0000
# MU_LFC_HGT MU_EL_HGT MU_LI MU_LI_M10
# 730.0000 8300.0000 -10.1119 -10.8539
# MU_WMAX MU_EL_TEMP MU_LCL_TEMP MU_LFC_TEMP
# 67.3784 -28.8000 17.7000 17.7000
# MU_MIXR MU_CAPE_500 MU_CAPE_500_M10 MU_CAPE_500_M10_PT
# 14.8759 1076.0322 389.3651 137.0814
# MU_CIN_500 MU_LI_500 MU_LI_500_M10 SB_CAPE
# 0.0000 -5.0417 -6.2346 2269.9257
# SB_CAPE_M10 SB_CAPE_M10_PT SB_02km_CAPE SB_03km_CAPE
# 998.1443 313.0747 247.9794 575.6293
# SB_HGL_CAPE SB_CIN SB_LCL_HGT SB_LFC_HGT
# 1616.5384 0.0000 730.0000 730.0000
# SB_EL_HGT SB_LI SB_LI_M10 SB_WMAX
# 8300.0000 -10.1119 -10.8539 67.3784
# SB_EL_TEMP SB_LCL_TEMP SB_LFC_TEMP SB_MIXR
# -28.8000 17.7000 17.7000 14.8759
# ML_CAPE ML_CAPE_M10 ML_CAPE_M10_PT ML_02km_CAPE
# 1646.0639 670.1001 225.2816 164.0798
# ML_03km_CAPE ML_HGL_CAPE ML_CIN ML_LCL_HGT
# 422.4290 1250.0221 0.0000 975.0000
# ML_LFC_HGT ML_EL_HGT ML_LI ML_LI_M10
# 975.0000 7900.0000 -7.6203 -8.5845
# ML_WMAX ML_EL_TEMP ML_LCL_TEMP ML_LFC_TEMP
# 57.3771 -26.4000 15.2500 15.2500
# ML_MIXR LR_0500m LR_01km LR_02km
# 13.0487 -10.0000 -10.0000 -10.0000
# LR_03km LR_04km LR_06km LR_16km
# -9.0476 -7.8571 -6.6667 -6.0000
# LR_26km LR_24km LR_36km LR_26km_MAX
# -5.0000 -5.7672 -4.2857 -5.7143
# LR_500700hPa LR_500800hPa LR_600800hPa FRZG_HGT
# -4.2857 -5.1807 -5.8333 2500.0000
# FRZG_wetbulb_HGT HGT_max_thetae_03km HGT_min_thetae_04km Delta_thetae
# 2275.0000 0.0000 3750.0000 28.0698
# Delta_thetae_min04km Thetae_01km Thetae_02km DCAPE
# 28.8346 330.5323 323.6191 598.3100
# Cold_Pool_Strength Wind_Index PRCP_WATER Moisture_Flux_02km
# 12.6322 33.9064 27.1046 30.4255
# RH_01km RH_02km RH_14km RH_25km
# 0.7291 0.7197 0.6452 0.5550
# RH_36km RH_HGL BS_0500m BS_01km
# 0.4436 0.4603 1.9172 3.8344
# BS_02km BS_03km BS_06km BS_08km
# 8.7821 12.6560 18.0055 17.4077
# BS_36km BS_26km BS_16km BS_18km
# 9.3693 13.3304 16.6478 20.2791
# BS_EFF_MU BS_EFF_SB BS_EFF_ML BS_SFC_to_M10
# 14.2232 14.2232 13.8968 15.5104
# BS_1km_to_M10 BS_2km_to_M10 BS_MU_LFC_to_M10 BS_SB_LFC_to_M10
# 13.6499 9.8830 14.0737 14.0737
# BS_ML_LFC_to_M10 BS_MW02_to_SM BS_MW02_to_RM BS_MW02_to_LM
# 13.6864 7.3040 10.1410 10.7870
# BS_HGL_to_SM BS_HGL_to_RM BS_HGL_to_LM MW_0500m
# 4.8934 7.7860 9.9885 2.3086
# MW_01km MW_02km MW_03km MW_06km
# 2.4251 3.3476 4.8003 7.8107
# MW_13km SRH_100m_RM SRH_250m_RM SRH_500m_RM
# 6.8389 4.2535 10.0537 19.7206
# SRH_1km_RM SRH_3km_RM SRH_36km_RM SRH_100m_LM
# 39.6346 152.5219 236.5901 1.5027
# SRH_250m_LM SRH_500m_LM SRH_1km_LM SRH_3km_LM
# 3.5518 6.9670 14.0023 -13.1308
# SRH_36km_LM SV_500m_RM SV_01km_RM SV_03km_RM
# -24.3790 0.0039 0.0039 0.0048
# SV_500m_LM SV_01km_LM SV_03km_LM MW_SR_500m_RM
# 0.0010 0.0011 -0.0014 10.0863
# MW_S
