skewt.axis {RadioSonde} | R Documentation |
Draws a SKEW-T, log p axis. This is the standard axis for interpreting atmospheric
sounding profiles like those collected by the instrument carried aloft by a
weather balloon (radiosonde). Use skewt.lines
or skewt.points
to layer information on top of the skewt axis.
skewt.axis(BROWN="brown", GREEN="green", redo=FALSE, ... )
BROWN |
Color of lines for both Temperature (solid skewed) and Pressure (dashed horizontal) |
GREEN |
Color of lines for dry adiabats (solid) and constant mixing ratio (dashed) |
redo |
flag to generate the adiabats, should skewt.data become
corrupt. The adiabats are the result of an iterative process
and to make smooth curves, you need a lot of points. Hence this
is time-consuming, so [FALSE] is the preferred value. |
... |
The usual plot parameters. |
Radiosondes record temperature, humidity and winds.
They can be lifted by weather balloons, dropped from aircraft,
there is even something called a glidersonde. The data collected by
radiosondes are plotted versus pressure/height to give details on the
vertical structure of the atmosphere. The type of plot is called a
SKEW-T, log p diagram.
‘skewt.axis’ creates the traditional axis for a SKEW-T, log p
diagram, including moist and dry adiabats, etc ...
Generating the necessary plot data is time-consuming, so there is an
option to read from a specific dataset that Doug will change.
Returns the par()$plt
values, which are used by
plotsonde
if you choose to plot the wind profile
adjacent to the skewt axis.
Tim Hoar
getsonde
,
plotsonde
,
skewt.lines
,
skewt.points
# draw a background, then # draw the temperature (with a solid line) in color 6 # draw the dewpoint in color 7 # overlay the temperature observations in a different color # you get the point ... # filename <- system.file("exampleData", "ExampleSonde.txt", package = "RadioSonde") sample.sonde <- getsonde(filename) skewt.axis( mar=c(5.1, 1.1, 2.1, 5.1) ) skewt.lines( sample.sonde$temp, sample.sonde$press, col = 6) skewt.lines( sample.sonde$dewpt, sample.sonde$press, col = 7) skewt.points(sample.sonde$temp, sample.sonde$press, col = 3) skewt.points(sample.sonde$dewpt, sample.sonde$press, col = 4) # # Changing the moist adiabats: you must edit the this-is-escaped-code{ function # directly and then capture the output in this-is-escaped-codenormal-bracket36bracket-normal to be used in # subsequent calls. skewt.data <- skewt.axis(redo=TRUE) skewt.axis() skewt.axis()