def phase_LogLinear(phase, slope=0.04):
    """A logLinear phase function, roughly appropriate for cometary nuclei.
    An H-G phase function is likely a better approximation.

    Parameters
    ----------
    phase : float or array
        Phase angle (degrees)
    slope : float, optional
        The slope for the phase function. Default 0.04.

    Returns
    -------
    phi : float or array
        Phase function evaluated at phase
    """
    return 10**(-0.4 * slope * phase)