def METHOD_NAME(self):
    with vega_spectrum.set('Bohlin2014'):
        vega = Vega.from_default()
        f = vega(0.55 * u.um)
        assert np.isclose(f.value, 3.546923511485616e-08)  # W/(m2 μm)