Function for plotting a customized version of the Skew-T diagram. Please note that drawing Skew-T may require increasing size or modifying aspect ratio of plotting window in order to provide readable results.
skewt_plot(
ptop = 100,
isoterms_col = "#d8be9b",
temp_stripes = FALSE,
mixing_ratio_col = "#8470FF90",
dry_adiabats_col = "#d6878750",
moist_adiabats_col = "#00FF0095",
deg45 = FALSE,
isotherm0 = TRUE,
close_par = TRUE,
...
)Arguments
- ptop
Pressure top level to be used for plotting diagram. Valid options: 200, 150, 100 (default) and 50 hPa
- isoterms_col
color to be used for drawing dry isoterms
- temp_stripes
logical, whether to draw color stripes for isotherms
- mixing_ratio_col
color to be used for drawing mixing ratio isolines and labels. If set to NA or empty string isolines are not drawn
- dry_adiabats_col
color to be used for drawing dry adiabats. If set to NA or not provided drawing lines skipped
- moist_adiabats_col
color to be used for drawing moist adiabats. If set to NA or not provided drawing lines skipped
- deg45
whether to preserve 45 degrees for diagonal isolines on Skew-T diagram regardless ploting window aspect ratio. [logical, default: FALSE]
- isotherm0
whether to deliminate 0 degree Celsius isother [logical, default: TRUE]
- close_par
if plot will be modified in next steps storing par settings is needed. This logical argument is turned on by default. If you want to modify Skew-T plot in next step set it to FALSE
- ...
additional (mostly graphical) parameters to be passed
Value
Draws a Skew-T log-p diagram
Examples
skewt_plot(ptop = 100)
skewt_plot(ptop = 150, temp_stripes = TRUE) # add color stripes for temperature
skewt_plot(ptop = 100, close_par = FALSE)
title("Your title")
mtext('WMO ID: 11035, 2011-08-23 1200 UTC', padj = -0.5, col = "white")
data("sounding_vienna")
attach(sounding_vienna)
#> The following objects are masked from sounding_vienna (pos = 3):
#>
#> altitude, dpt, pressure, temp, wd, ws
#> The following object is masked from package:datasets:
#>
#> pressure
output = sounding_export(pressure, altitude, temp, dpt, wd, ws)
skewt_lines(output$dpt, output$pressure,type='l',col='forestgreen',lwd = 2.5)
skewt_lines(output$temp,output$pressure,type='l',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)
