def quickPlot(self, output_file=None, **kwargs):
""" Do a quick plot of the sounding """
my_skew_t_figure = figure()
# Add an Skew-T axes to the Figure
my_skew_t_axes = my_skew_t_figure.add_subplot(111,
projection='skewx',
**kwargs)
pressure, temperature, dew_point_temperature = self.getCleanSounding()
# Add a profile to the Skew-T diagram
my_skew_t_axes.addProfile(pressure,
temperature,
dew_point_temperature,
hPa=True,
celsius=True,
method=0,
diagnostics=True)
if output_file is not None:
# Save the figure
my_skew_t_figure.save_fig(output_file)
else:
# Show the figure
my_skew_t_figure.show_plot()
"""
Simple SkewT Figure example from Sounding
"""
from __future__ import print_function, division
from SkewTplus.skewT import figure
from SkewTplus.sounding import sounding
# Load the sounding data
mySounding = sounding("./armSoundingExample.cdf")
# Create a Figure Manager
mySkewT_Figure = figure()
# Add an Skew-T axes to the Figure
mySkewT_Axes = mySkewT_Figure.add_subplot(111, projection="skewx")
# Extract the data from the Sounding
# The getCleanSounding method remove levels where invalid temperatura or pressure
# values are present
pressure, temperature, dewPointTemperature = mySounding.getCleanSounding()
# Add a profile to the Skew-T diagram
mySkewT_Axes.addProfile(
pressure,
temperature,
dewPointTemperature,
hPa=True,
'''
Plot two soundings in the same Skew-T diagram with out any parcel analysis
'''
from __future__ import print_function, division
from SkewTplus.skewT import figure
from SkewTplus.sounding import sounding
#Load the sounding data
mySounding1 = sounding("./bna_day1.txt")
mySounding2 = sounding("./bna_day2.txt")
# Create a Figure Manager with a suitable size for both plots
mySkewT_Figure = figure(figsize=(5, 6))
# Add the Skew-T axes to the Figure
mySkewT_Axes1 = mySkewT_Figure.add_subplot(111, projection='skewx', tmin=-40)
# Add one profile to the Skew-T diagram
# The line style is set to be a solid line and a label is added
# to the plot. Since the label is not None, a legend will be added
# automatically to the plot
mySkewT_Axes1.addProfile(*mySounding1.getCleanSounding(),
hPa=True,
celsius=True,
parcel=False,
label='Day 1',
linestyle='-')
# Add a second profile to the Skew-T diagram
mySounding.quickPlot()
'''
Simple SkewT Figure example from Sounding
'''
#Load the sounding data
mySounding = sounding("./exampleSounding.txt")
# Create a Figure Manager
mySkewT_Figure = figure()
# Add an Skew-T axes to the Figure
mySkewT_Axes = mySkewT_Figure.add_subplot(111, projection='skewx')
# Extract the data from the Sounding
# The getCleanSounding method remove levels where invalid temperatura or pressure
# values are present
pressure, temperature, dewPointTemperature = mySounding.getCleanSounding()
# Add a profile to the Skew-T diagram
mySkewT_Axes.addProfile(pressure,temperature, dewPointTemperature ,
hPa=True, celsius=True, method=0, diagnostics=True)
# Show the figure
mySkewT_Figure.show_figure()
parcel analysis, and plot them side to side
"""
from __future__ import print_function, division
from SkewTplus.skewT import figure
from SkewTplus.sounding import sounding
# Load the sounding data
mySounding1 = sounding("./bna_day1.txt")
mySounding2 = sounding("./bna_day2.txt")
# Create a Figure Manager with a suitable size for both plots
mySkewT_Figure = figure(figsize=(9, 5))
# Now we want to create two axes side to side
# Add the first Skew-T axes to the Figure
mySkewT_Axes1 = mySkewT_Figure.add_subplot(121, projection="skewx", tmin=-40)
# Extract the data from the Sounding
pressure, temperature, dewPointTemperature = mySounding1.getCleanSounding()
# Add a profile to the Skew-T diagram
mySkewT_Axes1.addProfile(
pressure,
temperature,