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 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)
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)