From c3c73d4f7752ec41888438b3294a423b12776a37 Mon Sep 17 00:00:00 2001 From: Gernot Wagner Date: Fri, 30 Mar 2018 10:47:20 -0400 Subject: [PATCH] Update bau.py Changed "(self.emission_by_decisions[n-1:n].mean())" to "(self.emission_by_decisions[n-1] + self.emission_by_decisions[n]) / 2". The former did not calculate the mean but only took the first element self.emission_by_decisions[n-1]. --- ezclimate/bau.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ezclimate/bau.py b/ezclimate/bau.py index c7f2fa2..61109de 100755 --- a/ezclimate/bau.py +++ b/ezclimate/bau.py @@ -127,7 +127,7 @@ def bau_emissions_setup(self, tree): for n in range(1, num_periods): self.emission_by_decisions[n] = self.emission_by_time(tree.decision_times[n]) - self.emission_per_period[n] = period_len[n-1] * (self.emission_by_decisions[n-1:n].mean()) + self.emission_per_period[n] = period_len[n-1] * (self.emission_by_decisions[n-1] + self.emission_by_decisions[n]) / 2 #the total increase in ghg level of 600 (from 400 to 1000) in the bau path is allocated over time self.emission_to_ghg = (self.ghg_end - self.ghg_start) * self.emission_per_period / self.emission_per_period.sum()