def solar_elevation_angle_refracted(h0):
    """    
    Computes the solar elevation angle corrected for refraction

    .. math ::
       h_{0}^{ref} = h_{0} + \\Delta h_{0}^{ref}

    where

    .. math ::
        \\Delta h_{0}^{ref}=0.61359\\left(0.1594+1.123h_{0}+
        0.065656h_{0}^{2}\\right)/\\left(
        1+28.9344h_{0}+277.3971h_{0}^{2}\\right)

    Parameters
    ----------
    h0 : float
        solar elevation angle
        :math:`h_{0}`
        [degrees]
        
    Returns
    -------
    h0ref : float
        solar elevation angle corrected for refrection
        :math:`h_{0}^{ref}`
        [degrees]
    """
    delta_h0ref = 0.061359 * (0.1594 + 1.123 * h0 + 0.065656 * h0 ** 2) / (
        1 + 28.9344 * h0 + 277.3971 * h0 ** 2)
    h0ref = h0 + delta_h0ref
    return h0ref
