artp_dict = dict()
iartp_dict = dict()
artp_dict['min'] = []
artp_dict['avg'] = []
artp_dict['max'] = []
iartp_dict['min'] = []
iartp_dict['avg'] = []
iartp_dict['max'] = []
for t in np.linspace(time_step, time_horizon, int(time_horizon * n_steps)):
iartp, artp = slp.compute_atp(
pollutant=pollutant,
rad_eff=rr_rad_eff,
th=t,
lifetime_range=True
)
artp_dict['min'].append(artp[0])
artp_dict['avg'].append(artp[1])
artp_dict['max'].append(artp[2])
iartp_dict['min'].append(iartp[0])
iartp_dict['avg'].append(iartp[1])
iartp_dict['max'].append(iartp[2])
print("\rt = {:5.2f}".format(t), flush=True, end=" ")
# Plot figure
sensitivity.plot_lifetime_range(artp_dict, iartp_dict, time_horizon)
# Compute temperature and precipitation potentials
if pol == 'CO2':
iartp, artp = co2.compute_atp(rad_eff=rr_rad_eff, th=th)
iarpp, slow_iarpp, fast_iarpp, arpp, slow_arpp, fast_arpp = co2.compute_app(
rad_eff=rad_eff,
rad_eff_a=rad_eff_a,
th=th,
rr_precip_avg=rr_precip_avg,
precip_avg=precip_avg)
else:
iartp, artp = slp.compute_atp(pollutant=pol,
rad_eff=rr_rad_eff,
th=th)
iarpp, slow_iarpp, fast_iarpp, arpp, slow_arpp, fast_arpp = slp.compute_app(
pollutant=pol,
rad_eff=rad_eff,
rad_eff_a=rad_eff_a,
th=th,
rr_precip_avg=rr_precip_avg,
precip_avg=precip_avg)
# Compute uncertainties
artp_std, arpp_std = propagation.get_potential_uncertainties(
pollutant=pol,
emission_region=emission_region,
response_regions=response_regions,
artp_dict = dict()
iartp_dict = dict()
artp_dict['artp'] = []
iartp_dict['iartp'] = []
artp_dict['no_c_artp'] = []
iartp_dict['no_c_iartp'] = []
artp_dict['no_erf_artp'] = []
iartp_dict['no_erf_iartp'] = []
artp_dict['no_c_no_erf_artp'] = []
iartp_dict['no_c_no_erf_iartp'] = []
for t in np.linspace(time_step, time_horizon, int(time_horizon * n_steps)):
iartp, artp = slp.compute_atp(pollutant=pollutant,
rad_eff=rr_rad_eff,
th=t)
no_c_iartp, no_c_artp = slp.compute_atp(pollutant=pollutant,
rad_eff=rr_rad_eff,
th=t,
c_scaling=False)
no_erf_iartp, no_erf_artp = slp.compute_atp(pollutant=pollutant,
rad_eff=rr_rad_eff,
th=t,
erf_scaling=False)
no_c_no_erf_iartp, no_c_no_erf_artp = slp.compute_atp(pollutant=pollutant,
rad_eff=rr_rad_eff,
th=t,