premise 2.1.0.dev0__tar.gz → 2.1.0.dev3__tar.gz

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/MANIFEST.in

@@ -39,6 +39,10 @@ include premise/data/geomap/*.yml
 include premise/data/geomap/*.yaml
 include premise/data/renewables/*.csv
 include premise/data/direct_air_capture/*.yaml
+include premise/data/transport/car/*.csv
+include premise/data/transport/bus/*.csv
+include premise/data/transport/truck/*.csv
+include premise/data/transport/two-wheeler/*.csv
 include premise/data/utils/*.csv
 include premise/data/utils/*.yml
 include premise/data/utils/*.yaml

{premise-2.1.0.dev0/premise.egg-info → premise-2.1.0.dev3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: premise
-Version: 2.1.0.dev0
+Version: 2.1.0.dev3
 Summary: Coupling IAM output to ecoinvent LCA database ecoinvent for prospective LCA
 Author-email: Romain Sacchi <romain.sacchi@psi.ch>, Alois Dirnaichner <dirnaichner@pik-potsdam.de>,  Chris Mutel <cmutel@gmail.com>
 Maintainer-email: Romain Sacchi <romain.sacchi@psi.ch>

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/__init__.py

@@ -1,5 +1,5 @@
 __all__ = ("NewDatabase", "clear_cache", "get_regions_definition")
-__version__ = (2, 1, 0, "dev0")
+__version__ = (2, 1, 0, "dev3")
 
 
 from premise.new_database import NewDatabase

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/activity_maps.py

@@ -93,11 +93,12 @@ def act_fltr(
         else:
             filters.append(ws.contains(field, value))
 
-    for field, value in mask.items():
-        if isinstance(value, list):
-            filters.extend([ws.exclude(ws.contains(field, v)) for v in value])
-        else:
-            filters.append(ws.exclude(ws.contains(field, value)))
+    if mask:
+        for field, value in mask.items():
+            if isinstance(value, list):
+                filters.extend([ws.exclude(ws.contains(field, v)) for v in value])
+            else:
+                filters.append(ws.exclude(ws.contains(field, value)))
 
     return list(ws.get_many(database, *filters))
 
@@ -277,8 +278,27 @@ class InventorySet:
 
         database = database or self.database
 
+        names = []
+
+        for entry in filtr.values():
+            if "fltr" in entry:
+                if isinstance(entry["fltr"], dict):
+                    if "name" in entry["fltr"]:
+                        names.extend(entry["fltr"]["name"])
+                elif isinstance(entry["fltr"], list):
+                    names.extend(entry["fltr"])
+                else:
+                    names.append(entry["fltr"])
+
+        subset = list(
+            ws.get_many(
+                database,
+                ws.either(*[ws.contains("name", name) for name in names]),
+            )
+        )
+
         techs = {
-            tech: act_fltr(database, fltr.get("fltr"), fltr.get("mask"))
+            tech: act_fltr(subset, fltr.get("fltr"), fltr.get("mask"))
             for tech, fltr in filtr.items()
         }
 

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/biomass.py

@@ -164,31 +164,61 @@ class Biomass(BaseTransformation):
             ]
 
             for biomass_type, biomass_act in biomass_map.items():
-                total_prod_vol = np.clip(
-                    (
-                        self.iam_data.production_volumes.sel(
-                            variables=available_biomass_vars, region=region
-                        )
-                        .interp(year=self.year)
-                        .sum(dim="variables")
-                    ),
-                    1e-6,
-                    None,
-                )
-
-                if biomass_type in available_biomass_vars:
-                    share = np.clip(
+                if self.year in self.iam_data.production_volumes.coords["year"].values:
+                    total_prod_vol = np.clip(
+                        (
+                            self.iam_data.production_volumes.sel(
+                                variables=available_biomass_vars,
+                                region=region,
+                                year=self.year,
+                            ).sum(dim="variables")
+                        ),
+                        1e-6,
+                        None,
+                    )
+                else:
+                    total_prod_vol = np.clip(
                         (
                             self.iam_data.production_volumes.sel(
-                                variables=biomass_type, region=region
+                                variables=available_biomass_vars, region=region
                             )
                             .interp(year=self.year)
-                            .sum()
-                            / total_prod_vol
-                        ).values.item(0),
-                        0,
-                        1,
+                            .sum(dim="variables")
+                        ),
+                        1e-6,
+                        None,
                     )
+
+                if biomass_type in available_biomass_vars:
+                    if (
+                        self.year
+                        in self.iam_data.production_volumes.coords["year"].values
+                    ):
+                        share = np.clip(
+                            (
+                                self.iam_data.production_volumes.sel(
+                                    variables=biomass_type,
+                                    region=region,
+                                    year=self.year,
+                                ).sum()
+                                / total_prod_vol
+                            ).values.item(0),
+                            0,
+                            1,
+                        )
+                    else:
+                        share = np.clip(
+                            (
+                                self.iam_data.production_volumes.sel(
+                                    variables=biomass_type, region=region
+                                )
+                                .interp(year=self.year)
+                                .sum()
+                                / total_prod_vol
+                            ).values.item(0),
+                            0,
+                            1,
+                        )
                 elif (
                     self.system_model == "consequential"
                     and biomass_type == "biomass - residual"

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/cement.py

@@ -10,7 +10,16 @@ of the wurst database to the newly created cement markets.
 from collections import defaultdict
 
 from .logger import create_logger
-from .transformation import BaseTransformation, Dict, IAMDataCollection, List, np, ws
+from .transformation import (
+    BaseTransformation,
+    Dict,
+    IAMDataCollection,
+    InventorySet,
+    List,
+    calculate_input_energy,
+    np,
+    ws,
+)
 from .validation import CementValidation
 
 logger = create_logger("cement")
@@ -98,6 +107,18 @@ class Cement(BaseTransformation):
         )
         self.version = version
 
+        mapping = InventorySet(self.database)
+        self.cement_fuels_map: dict = mapping.generate_cement_fuels_map()
+
+        # reverse the fuel map to get a mapping from ecoinvent to premise
+        self.fuel_map_reverse: dict = {}
+
+        self.fuel_map: dict = mapping.generate_fuel_map()
+
+        for key, value in self.fuel_map.items():
+            for v in list(value):
+                self.fuel_map_reverse[v] = key
+
     def fetch_current_energy_details(self, dataset):
         # Fetches the current energy consumption per ton of clinker
 
@@ -115,7 +136,13 @@ class Cement(BaseTransformation):
 
                 # Calculate the energy once for the given exc
                 input_energy = (
-                    self.calculate_input_energy(exc["name"], exc["amount"], exc["unit"])
+                    calculate_input_energy(
+                        fuel_name=exc["name"],
+                        fuel_amount=exc["amount"],
+                        fuel_unit=exc["unit"],
+                        fuels_specs=self.fuels_specs,
+                        fuel_map_reverse=self.fuel_map_reverse,
+                    )
                     * 1000
                 )
 
@@ -469,23 +496,30 @@ class Cement(BaseTransformation):
             self.add_to_index(new_dataset)
 
         # cement markets
-        markets = ws.get_many(
-            self.database,
-            ws.contains("name", "market for cement"),
-            ws.contains("reference product", "cement"),
-            ws.doesnt_contain_any("name", ["factory", "tile", "sulphate", "plaster"]),
-            ws.doesnt_contain_any("location", self.regions),
+        markets = list(
+            ws.get_many(
+                self.database,
+                ws.contains("name", "market for cement"),
+                ws.contains("reference product", "cement"),
+                ws.doesnt_contain_any(
+                    "name", ["factory", "tile", "sulphate", "plaster"]
+                ),
+                ws.doesnt_contain_any("location", self.regions),
+            )
         )
 
-        markets = {(m["name"], m["reference product"]) for m in markets}
+        unique_markets = list(
+            set([(m["name"], m["reference product"]) for m in markets])
+        )
 
         new_datasets = []
 
-        for dataset in markets:
+        for dataset in unique_markets:
             new_cement_markets = self.fetch_proxies(
                 name=dataset[0],
                 ref_prod=dataset[1],
                 production_variable="cement, dry feed rotary kiln",
+                subset=markets,
             )
 
             # add to log
@@ -499,24 +533,31 @@ class Cement(BaseTransformation):
         self.database.extend(new_datasets)
 
         # cement production
-        production = ws.get_many(
-            self.database,
-            ws.contains("name", "cement production"),
-            ws.contains("reference product", "cement"),
-            ws.doesnt_contain_any("name", ["factory", "tile", "sulphate", "plaster"]),
+        production = list(
+            ws.get_many(
+                self.database,
+                ws.contains("name", "cement production"),
+                ws.contains("reference product", "cement"),
+                ws.doesnt_contain_any(
+                    "name", ["factory", "tile", "sulphate", "plaster"]
+                ),
+            )
         )
 
-        production = {(p["name"], p["reference product"]) for p in production}
+        reduced_production = list(
+            set([(p["name"], p["reference product"]) for p in production])
+        )
 
         new_datasets = []
 
-        for dataset in production:
+        for dataset in reduced_production:
             # Fetch proxy datasets (one per IAM region)
             # Delete old datasets
             new_cement_production = self.fetch_proxies(
                 name=dataset[0],
                 ref_prod=dataset[1],
                 production_variable="cement, dry feed rotary kiln",
+                subset=production,
             )
 
             # add to log

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/clean_datasets.py

@@ -10,6 +10,7 @@ from pathlib import Path
 from typing import Dict, List, Tuple, Union
 
 import bw2io
+import numpy as np
 import wurst
 import yaml
 from bw2data.database import DatabaseChooser
@@ -34,23 +35,17 @@ def remove_uncertainty(database):
     :param database:
     :return:
     """
-
-    keys_to_remove = [
-        "loc",
-        "scale",
-        "pedigree",
-        "minimum",
-        "maximum",
-    ]
+    uncertainty_keys = ["scale", "shape", "minimum", "maximum"]
+    nan_value = np.nan
 
     for dataset in database:
-        for exc in dataset["exchanges"]:
-            if "uncertainty type" in exc:
-                if exc["uncertainty type"] != 0:
-                    exc["uncertainty type"] = 0
-                    for key in keys_to_remove:
-                        if key in exc:
-                            del exc[key]
+        for exchange in dataset["exchanges"]:
+            if not exchange.get("preserve uncertainty", False):
+                exchange["uncertainty type"] = 0
+                exchange["loc"] = exchange["amount"]
+                for key in uncertainty_keys:
+                    exchange[key] = nan_value
+
     return database
 
 

{premise-2.1.0.dev0 → premise-2.1.0.dev3}/premise/data/fuels/fuel_markets.yml

@@ -12,7 +12,7 @@ natural gas:
   name: market for natural gas, high pressure
   reference product: natural gas, high pressure
   unit: cubic meter
-  lhv: 47.5
+  lhv: 36
 hydrogen:
   name: market for hydrogen, gaseous
   reference product: hydrogen, gaseous