def test_angstrom_aod(): """Test Angstrom turbidity model functions.""" aod550 = 0.15 aod1240 = 0.05 alpha = atmosphere.angstrom_alpha(aod550, 550, aod1240, 1240) assert np.isclose(alpha, 1.3513924317859232) aod700 = atmosphere.angstrom_aod_at_lambda(aod550, 550, alpha) assert np.isclose(aod700, 0.10828110997681031)
def to_wvl(wvl, aod550, alpha=1.14): """ Convert AOD from tau550 to specified wavelength. """ # aod = aod550 * ((wvl / 550.0) ** (-alpha)) return angstrom_aod_at_lambda(aod550, 550.0, alpha, wvl)