Source code for bsrn.qc.erl

"""
BSRN Level 2 (extremely rare) quality tests.
"""

import numpy as np
import pandas as pd




def ghi_erl_test(ghi, zenith, bni_extra):
    """
    Check global horizontal irradiance (GHI, $G_h$) against extremely rare limits [1]_.

    Parameters
    ----------
    ghi : numeric or Series
        Global horizontal irradiance ($G_h$). [W/m^2]
    zenith : numeric or Series
        Solar zenith angle ($Z$). [degrees]
    bni_extra : numeric or Series
        Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]

    Returns
    -------
    flags : Series or ndarray
        Boolean flags (True = Pass). [bool]

    References
    ----------
    .. [1] Long, C. N., & Shi, Y. (2008). An automated quality assessment
       and control algorithm for surface radiation measurements. The Open
       Atmospheric Science Journal, 2(1), 23-37.
    """
    mu0 = np.cos(np.radians(zenith))
    mu0 = np.maximum(mu0, 0)

    # Upper limit: 1.2 * E_0n * mu0^1.2 + 50
    upper_limit = 1.2 * bni_extra * (mu0 ** 1.2) + 50
    # Lower limit: -2 W/m^2
    lower_limit = -2

    return (ghi >= lower_limit) & (ghi <= upper_limit)





def bni_erl_test(bni, zenith, bni_extra):
    """
    Check beam normal irradiance (BNI, $B_n$) against extremely rare limits [1]_.

    Parameters
    ----------
    bni : numeric or Series
        Beam normal irradiance ($B_n$). [W/m^2]
    zenith : numeric or Series
        Solar zenith angle ($Z$). [degrees]
    bni_extra : numeric or Series
        Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]

    Returns
    -------
    flags : Series or ndarray
        Boolean flags (True = Pass). [bool]

    References
    ----------
    .. [1] Long, C. N., & Shi, Y. (2008). An automated quality assessment
       and control algorithm for surface radiation measurements. The Open
       Atmospheric Science Journal, 2(1), 23-37.
    """
    mu0 = np.cos(np.radians(zenith))
    mu0 = np.maximum(mu0, 0)

    # Upper limit: 0.95 * E_0n * mu0^0.2 + 10
    upper_limit = 0.95 * bni_extra * (mu0 ** 0.2) + 10
    # Lower limit: -2 W/m^2
    lower_limit = -2

    return (bni >= lower_limit) & (bni <= upper_limit)





def dhi_erl_test(dhi, zenith, bni_extra):
    """
    Check diffuse horizontal irradiance (DHI, $D_h$) against extremely rare limits [1]_.

    Parameters
    ----------
    dhi : numeric or Series
        Diffuse horizontal irradiance ($D_h$). [W/m^2]
    zenith : numeric or Series
        Solar zenith angle ($Z$). [degrees]
    bni_extra : numeric or Series
        Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]

    Returns
    -------
    flags : Series or ndarray
        Boolean flags (True = Pass). [bool]

    References
    ----------
    .. [1] Long, C. N., & Shi, Y. (2008). An automated quality assessment
       and control algorithm for surface radiation measurements. The Open
       Atmospheric Science Journal, 2(1), 23-37.
    """
    mu0 = np.cos(np.radians(zenith))
    mu0 = np.maximum(mu0, 0)

    # Upper limit: 0.75 * E_0n * mu0^1.2 + 30
    upper_limit = 0.75 * bni_extra * (mu0 ** 1.2) + 30
    # Lower limit: -2 W/m^2
    lower_limit = -2

    return (dhi >= lower_limit) & (dhi <= upper_limit)





def lwd_erl_test(lwd):
    """
    Check downward longwave radiation (LWD, $L_d$) against extremely rare limits [1]_.

    Parameters
    ----------
    lwd : numeric or Series
        Downward longwave radiation ($L_d$). [W/m^2]

    Returns
    -------
    flags : Series or ndarray
        Boolean flags (True = Pass). [bool]

    References
    ----------
    .. [1] Long, C. N., & Shi, Y. (2008). An automated quality assessment
       and control algorithm for surface radiation measurements. The Open
       Atmospheric Science Journal, 2(1), 23-37.
    """
    # Range: [60, 500] W/m^2
    return (lwd >= 60) & (lwd <= 500)