def blackbody(nu, ref_freq=353.0):
    '\n    The ratio of the blackbody function for dust at frequency nu\n    over the value for reference frequency ref_freq\n\n    Arguments\n    ---------\n    nu : float\n        Frequency in GHz.\n    ref_freq : float\n        Reference frequency in GHz.\n\n    Returns\n    -------\n    blackbody_ratio : float\n        B(nu, T_dust) / B(nu_ref, T_dust)\n    '
    k = 1.38064852e-23
    h = 6.62607004e-34
    T = 19.6
    nu_ref = (ref_freq * 1000000000.0)
    nu *= 1000000000.0
    x = (((h * nu) / k) / T)
    x_ref = (((h * nu_ref) / k) / T)
    return ((((x ** 3) / (x_ref ** 3)) * (np.exp(x_ref) - 1)) / (np.exp(x) - 1))
The ratio of the blackbody function for dust at frequency nu
over the value for reference frequency ref_freq

Arguments
---------
nu : float
    Frequency in GHz.
ref_freq : float
    Reference frequency in GHz.

Returns
-------
blackbody_ratio : float
    B(nu, T_dust) / B(nu_ref, T_dust)
def blackbody(nu, ref_freq=353.0):
    '\n    The ratio of the blackbody function for dust at frequency nu\n    over the value for reference frequency ref_freq\n\n    Arguments\n    ---------\n    nu : float\n        Frequency in GHz.\n    ref_freq : float\n        Reference frequency in GHz.\n\n    Returns\n    -------\n    blackbody_ratio : float\n        B(nu, T_dust) / B(nu_ref, T_dust)\n    '
    k = 1.38064852e-23
    h = 6.62607004e-34
    T = 19.6
    nu_ref = (ref_freq * 1000000000.0)
    nu *= 1000000000.0
    x = (((h * nu) / k) / T)
    x_ref = (((h * nu_ref) / k) / T)
    return ((((x ** 3) / (x_ref ** 3)) * (np.exp(x_ref) - 1)) / (np.exp(x) - 1))