Side_VS_HV01
How to use this plot
Make sure you have the required datafields (air_pressure, air_temperature, eastward_wind, northward_wind)
You can use it as is by appending this code into your mswms_settings.py:
from mslib.mswms.mpl_vsec_styles import VS_HorizontalVelocityStyle_01
register_vertical_layers = [] if not register_vertical_layers else register_vertical_layers
register_vertical_layers.append((VS_HorizontalVelocityStyle_01, [next(iter(data))]))
If you want to modify the plot
Download this
file
Put this file into your mswms_settings.py directory, e.g. ~/mss
Append this code into your mswms_settings.py:
from Side_VS_HV01 import VS_HorizontalVelocityStyle_01
register_vertical_layers = [] if not register_vertical_layers else register_vertical_layers
register_vertical_layers.append((VS_HorizontalVelocityStyle_01, [next(iter(data))]))
Plot Code
"""
This file is part of MSS.
:copyright: Copyright 2021-2024 by the MSS team, see AUTHORS.
:license: APACHE-2.0, see LICENSE for details.
Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at
http://www.apache.org/licenses/LICENSE-2.0
Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
"""
import warnings
import matplotlib
import matplotlib.pyplot as plt
import mpl_toolkits.axes_grid1.inset_locator
from matplotlib import patheffects
import numpy as np
from mslib.mswms.mpl_vsec import AbstractVerticalSectionStyle
from mslib.mswms.utils import make_cbar_labels_readable
import mslib.mswms.generics as generics
from mslib.utils import thermolib
from mslib.utils.units import convert_to
class VS_HorizontalVelocityStyle_01(AbstractVerticalSectionStyle):
"""
Vertical sections of horizontal velocity.
"""
name = "VS_HV01"
title = "Horizontal Wind (m/s) Vertical Section"
abstract = "Horizontal wind speed (m/s) with temperature (K) and potential temperature (K)"
# NOTE: This style is used for the flight performance computations. Make sure
# that it always requests air_pressure, air_temperature, eastward_wind,
# northward_wind! (mr, 2012Nov09)
# Variables with the highest number of dimensions first (otherwise
# MFDatasetCommonDims will throw an exception)!
required_datafields = [
("ml", "air_pressure", "Pa"),
("ml", "air_temperature", "K"),
("ml", "eastward_wind", "m/s"),
("ml", "northward_wind", "m/s")]
def _prepare_datafields(self):
"""
Computes potential temperature from pressure and temperature and
total horizontal wind speed.
"""
self.data['air_potential_temperature'] = thermolib.pot_temp(
self.data['air_pressure'], self.data['air_temperature'])
self.data["horizontal_wind"] = np.hypot(
self.data["eastward_wind"], self.data["northward_wind"])
def _plot_style(self):
"""
Make a horizontal velocity vertical section with temperature/potential
temperature overlay.
"""
ax = self.ax
curtain_p = self.data["air_pressure"]
curtain_t = self.data["air_temperature"]
curtain_pt = self.data["air_potential_temperature"]
curtain_v = self.data["horizontal_wind"]
# Contour spacing for temperature lines.
delta_t = 2 if (np.log(self.p_bot) - np.log(self.p_top)) < 2.2 else 4
delta_pt = 5 if (np.log(self.p_bot) - np.log(self.p_top)) < 2.2 else 10
wind_contours = np.arange(10, 90, 5)
# Filled contour plot of relative humidity.
# INFO on COLORMAPS:
# http://matplotlib.sourceforge.net/examples/pylab_examples/show_colormaps.html
cs = ax.contourf(self.horizontal_coordinate, curtain_p, curtain_v,
wind_contours, cmap=plt.cm.hot_r)
# Contour line plot of relative humidity.
cs_v1 = ax.contour(self.horizontal_coordinate,
curtain_p, curtain_v, np.arange(5, 90, 5),
colors="orange", linestyles="solid", linewidths=1)
ax.clabel(cs_v1, fontsize=8, fmt='%.0f')
cs_v2 = ax.contour(self.horizontal_coordinate,
curtain_p, curtain_v, np.arange(95, 150, 5),
colors="black", linestyles="solid", linewidths=0.5)
ax.clabel(cs_v2, fontsize=8, fmt='%.0f')
# Contour line plot of temperature.
cs_t = ax.contour(self.horizontal_coordinate,
curtain_p, curtain_t, np.arange(236, 330, delta_t),
colors='green', linestyles='solid', linewidths=1)
ax.clabel(cs_t, fontsize=8, fmt='%.0f')
cs_t = ax.contour(self.horizontal_coordinate, curtain_p, curtain_t, [234],
colors='green', linestyles='solid', linewidths=2)
ax.clabel(cs_t, fontsize=8, fmt='%.0f')
cs_t = ax.contour(self.horizontal_coordinate,
curtain_p, curtain_t, np.arange(160, 232, delta_t),
colors='green', linestyles='dashed', linewidths=1)
ax.clabel(cs_t, fontsize=8, fmt='%.0f')
# Contour line plot of potential temperature.
cs_pt = ax.contour(self.horizontal_coordinate, curtain_p, curtain_pt,
np.arange(200, 700, delta_pt), colors='grey',
linestyles='solid', linewidths=1)
ax.clabel(cs_pt, fontsize=8, fmt='%.0f')
# Pressure decreases with index, i.e. orography is stored at the
# zero-p-index (data field is flipped in mss_plot_driver.py if
# pressure increases with index).
self._latlon_logp_setup(orography=curtain_p[0, :])
self.add_colorbar(cs, "Horizontal wind speed (m/s)")