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 update_data(cls, data):
data['pwat'] = data[TCWV] / 10.0 # convert kg / m^2 to cm
data['alpha'] = angstrom_alpha(data[AOD1240], 1240.0,
data[AOD550], 550.0)
# -np.log(data[AOD1240] / data[AOD550]) / np.log(1240.0 / 550.0)
data['tau380'] = cls.to_wvl(380.0, data[AOD550], data['alpha'])
data['tau500'] = cls.to_wvl(500.0, data[AOD550], data['alpha'])
data['tau700'] = cls.to_wvl(700.0, data[AOD550], data['alpha'])
return data
def update_data(cls, data):
data['pwat'] = data[TCWV] / 10.0 # convert kg / m^2 to cm
data['alpha'] = angstrom_alpha(data[AOD1240], 1240.0, data[AOD550],
550.0)
# -np.log(data[AOD1240] / data[AOD550]) / np.log(1240.0 / 550.0)
data['tau380'] = cls.to_wvl(380.0, data[AOD550], data['alpha'])
data['tau500'] = cls.to_wvl(500.0, data[AOD550], data['alpha'])
data['tau700'] = cls.to_wvl(700.0, data[AOD550], data['alpha'])
return data
