def net_radiation(r0, ra_24, l_net, int_wm2):
    """
    Computes the net radiation

    .. math ::
        Q^{*} = \\left[\\left(1-\\alpha_{0}\\right)S^{\\downarrow}-L^{*}-I\\right]

    Parameters
    ----------
    r0 : float
        albedo
        :math:`\\alpha_{0}`
        [-]
    ra_24 : float
        daily solar radiation
        :math:`S^{\\downarrow}`
        [Wm-2]
    l_net : float
        daily net longwave radiation
        :math:`L^{*}`
        [wm-2]
    int_wm2 : float
        interception
        :math:`I`
        [Wm-2]

    Returns
    -------
    rn_24 : float
        daily net radiation
        :math:`Q^{*}`
        [Wm-2]

    Examples
    --------
    >>> import ETLook.radiation as rad
    >>> rad.net_radiation(r0=0.10, ra_24=123., l_net=24., int_wm2=0)
    86.7
    """
    return (1 - r0) * ra_24 - l_net - int_wm2
