diff --git a/src/components/kalman_component/print_posterior_emissions.py b/src/components/kalman_component/print_posterior_emissions.py index 024cb017..a4e7e347 100644 --- a/src/components/kalman_component/print_posterior_emissions.py +++ b/src/components/kalman_component/print_posterior_emissions.py @@ -41,8 +41,6 @@ def print_posterior_emissions(config_path, period_number, base_directory): posterior_sf = xr.load_dataset(post_sf_path) posterior_emis_ds = get_posterior_emissions(hemco_emis, posterior_sf) posterior_emis = posterior_emis_ds["EmisCH4_Total"].isel(time=0, drop=True) - # Recall that posterior_sf from archive_sf multiplies sf over the full inversion record - posterior_emis = posterior_sf * hemco_emis total_emis = sum_total_emissions(posterior_emis, areas, mask) # Print