Note
Go to the end to download the full example code.
ASTM Spectra#
Standardized Sun spectra at different airmasses

1347.9343199999998
import os
import matplotlib.pyplot as plt
import numpy as np
import mirage as mr
import mirage.vis as mrv
sun_spectra = np.loadtxt(
os.path.join(os.environ['DATADIR'], 'ASTMG173.csv'), skiprows=2, delimiter=','
)
lambdas_nm = sun_spectra[:, 0]
sun_vacuum_w_per_m2_nm = sun_spectra[:, 1]
sun_global_w_per_m2_nm = sun_spectra[:, 2]
sun_direct_w_per_m2_nm = sun_spectra[:, 3]
ss_proof = mr.sun_spectrum(lambdas_nm)
i0 = np.trapz(sun_vacuum_w_per_m2_nm, lambdas_nm)
print(i0)
plt.plot(lambdas_nm, sun_vacuum_w_per_m2_nm)
plt.plot(lambdas_nm, sun_global_w_per_m2_nm)
plt.plot(lambdas_nm, sun_direct_w_per_m2_nm)
plt.plot(lambdas_nm, ss_proof)
mrv.texit(
'Sun Irradiance Spectra',
'Wavelength [nm]',
r'$\left[ \frac{W}{m^2 \cdot m} \right]$',
['AM0', 'AM1.5 Global', 'AM1.5 Direct', 'PROOF zero airmass'],
)
plt.show()
Total running time of the script: (0 minutes 0.159 seconds)