Linear_LS_W01

How to use this plot

Make sure you have the required datafields (air_pressure, air_temperature, lagrangian_tendency_of_air_pressure)

You can use it as is by appending this code into your mswms_settings.py:

from mslib.mswms.mpl_lsec_styles import LS_VerticalVelocityStyle_01
register_linear_layers = [] if not register_linear_layers else register_linear_layers
register_linear_layers.append((LS_VerticalVelocityStyle_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 Linear_LS_W01 import LS_VerticalVelocityStyle_01
register_linear_layers = [] if not register_linear_layers else register_linear_layers
register_linear_layers.append((LS_VerticalVelocityStyle_01, [next(iter(data))]))

Plot Code

"""
    This file is part of MSS.

    :copyright: Copyright 2021-2023 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 numpy as np
from mslib.mswms.mpl_lsec import AbstractLinearSectionStyle
import mslib.utils.thermolib as thermolib
from mslib.utils.units import convert_to

class LS_VerticalVelocityStyle_01(AbstractLinearSectionStyle):
    """
    Linear plot of vertical velocity.
    """

    name = "LS_W01"
    title = "Vertical Velocity (cm/s) Linear Plot"
    abstract = "Vertical velocity (cm/s)"

    # 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", "lagrangian_tendency_of_air_pressure", "Pa/s")]

    def _prepare_datafields(self):
        """
        Computes vertical velocity in cm/s.
        """
        self.data["upward_wind"] = convert_to(
            thermolib.omega_to_w(self.data["lagrangian_tendency_of_air_pressure"],
                                 self.data['air_pressure'], self.data["air_temperature"]),
            "m/s", "cm/s")
        self.variable = "upward_wind"
        self.y_values = self.data[self.variable]
        self.unit = "cm/s"