interpolation_mode parameter for technodata

Author: tsmbland

docs/inputs/toml.rst

@@ -292,6 +292,13 @@ A sector accepts these attributes:
 
    Defaults to "last".
 
+*interpolation*
+   Interpolation method used to fill missing years in the *technodata* (defaults to "linear").
+   Available interpolation methods depend on the underlying `scipy method's kind attribute`_.
+   Years outside the data range will always be back/forward filled with the closest available data.
+
+   .. _scipy method's kind attribute: https://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.interp1d.html
+
 *dispatch_production*
    The method used to calculate supply of commodities after investments have been made.
 

src/muse/readers/toml.py

@@ -646,7 +646,7 @@ def read_technodata(
     time_framework: Sequence[int] | None = None,
     commodities: str | Path | None = None,
     regions: Sequence[str] | None = None,
-    **kwargs,
+    interpolation_mode: str = "linear",
 ) -> xr.Dataset:
     """Helper function to create technodata for a given sector."""
     from muse.readers.csv import read_technologies, read_trade
@@ -741,5 +741,5 @@ def read_technodata(
         technologies["year"] = "year", years
 
     year = sorted(set(time_framework).union(technologies.year.data.tolist()))
-    technologies = technologies.interp(year=year, **kwargs)
+    technologies = technologies.interp(year=year, method=interpolation_mode)
     return technologies

src/muse/sectors/sector.py

@@ -37,9 +37,15 @@ def factory(cls, name: str, settings: Any) -> Sector:
             raise RuntimeError(f"Missing 'subsectors' section in sector {name}")
         if len(sector_settings["subsectors"]._asdict()) == 0:
             raise RuntimeError(f"Empty 'subsectors' section in sector {name}")
+        interpolation_mode = sector_settings.pop("interpolation_mode", "linear")
 
         # Read technologies
-        technologies = read_technodata(settings, name, settings.time_framework)
+        technologies = read_technodata(
+            settings,
+            name,
+            settings.time_framework,
+            interpolation_mode=interpolation_mode,
+        )
 
         # Create subsectors
         subsectors = [