Unidata / MetPy

    y = y[keep_idx]

    return x, y

@exporter.export

def cape(pressure, temperature, dewpt, parcel_temperature=None, use_virtual_temperature=False):

    r"""Calculate convective available potential energy (CAPE).

        Description

        Parameters

        ----------

        pressure : `pint.Quantity`

            The atmospheric pressure level(s) of interest. The first entry should be the starting

            point pressure.

        temperature : `pint.Quantity`

            The starting temperature

        dewpt : `pint.Quantity`

            The starting dew point

        parcel_temperature : `pint.Quantity`

            The temperature of the parcel

        Returns

        -------

        `pint.Quantity`

            Convective available potential energy (CAPE).

        See Also

        --------

        virtual_temperature, find_intersections, lfc, el

        """

    if parcel_temperature is None:

        parcel_temperature = parcel_profile(pressure, temperature[0], dewpt[0])

    if use_virtual_temperature:

        temperature = virtual_temperature(pressure, temperature, dewpt)

        parcel_temperature = virtual_temperature(pressure, parcel_temperature, dewpt)

    y = (parcel_temperature - temperature).to(units.degK)

    x,y = _find_append_zero_crossings(np.copy(pressure),y)

    # Clip out negative area (temperature parcel < temperature environment)

    y[y <= 0 * units.degK] = 0 * units.degK

    # Only use data between the LFC and EL for calculation

    lfc_pressure = lfc(pressure, temperature, dewpt)[0].magnitude

    el_pressure = el(pressure, temperature, dewpt)[0].magnitude

    p_mask = (x <= lfc_pressure) & (x >= el_pressure)

    x = x[p_mask]

    y = y[p_mask]

    return (Rd * (np.trapz(y, np.log(x)) * units.degK)).to(units('J/kg'))

@exporter.export

    return (Rd * (np.trapz(y, np.log(x)) * units.degK)).to(units('J/kg'))

@exporter.export

def cin(pressure, temperature, dewpt, parcel_temperature=None, use_virtual_temperature=False):

    r"""Calculate convective inhibition (CIN).

        Description

        Parameters

        ----------

        pressure : `pint.Quantity`

            The atmospheric pressure level(s) of interest. The first entry should be the starting

            point pressure.

        temperature : `pint.Quantity`

            The starting temperature

        dewpt : `pint.Quantity`

            The starting dew point

        parcel_temperature : `pint.Quantity`

            The temperature of the parcel

        Returns

        -------

        `pint.Quantity`

            Convective inhibition (CIN).

        See Also

        --------

        virtual_temperature, find_intersections, lfc

        """

    if parcel_temperature is None:

        parcel_temperature = parcel_profile(pressure, temperature[0], dewpt[0])

    if use_virtual_temperature:

        temperature = virtual_temperature(pressure, temperature, dewpt)

        parcel_temperature = virtual_temperature(pressure, parcel_temperature, dewpt)

    y = (parcel_temperature - temperature).to(units.degK)

    x,y = _find_append_zero_crossings(np.copy(pressure),y)

    # Clip out positive area (temperature parcel > temperature environment)

    y[y >= 0 * units.degK] = 0 * units.degK

    # Only use data between the surface and the LFC for calculation

    lfc_pressure = lfc(pressure, temperature, dewpt)[0].magnitude

    p_mask = (x >= lfc_pressure)

    x = x[p_mask]

    y = y[p_mask]

    return (Rd * (np.trapz(y, np.log(x)) * units.degK)).to(units('J/kg'))

@exporter.export