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