import numpy as np
from shared import *
import crl, crl.tables

plot_setup(width=0.7)

for T in reversed([5000, 5500, 6500, 7500]):
	wvl = np.linspace(380, 780, 100)
	n = crl.CIEDaylightSpectrum(T).Yi(wvl)

	plt.plot(wvl, n,
	         label="$D_{%d}$" % (T // 100))

plt.xlabel("$\\lambda$ [nm]")
plt.ylabel("$n_\\lambda(\\lambda)$")
plt.xlim(380, 780)
plt.ylim(0, 200)
plt.legend()
plot_grid()
plot_export()