detquantlib 3.15.0__tar.gz → 4.0.0__tar.gz

{detquantlib-3.15.0 → detquantlib-4.0.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: detquantlib
-Version: 3.15.0
+Version: 4.0.0
 Summary: An internal library containing functions and classes that can be used across Quant models.
 License-File: LICENSE.txt
 Author: Dynamic Energy Trading
@@ -80,6 +80,9 @@ poetry run python docs.py
 
 Classes:
 
+- `Battery`:
+  - `from detquantlib.assets import Battery`
+  - `from detquantlib.assets.battery import Battery`
 - `DetDatabase`:
   - `from detquantlib.data import DetDatabase`
   - `from detquantlib.data.databases.detdatabase import DetDatabase`
@@ -98,6 +101,9 @@ Classes:
 - `Sftp`:
   - `from detquantlib.data import Sftp`
   - `from detquantlib.data.sftp.sftp import Sftp`
+- `WindTurbine`:
+  - `from detquantlib.assets import WindTurbine`
+  - `from detquantlib.assets.wind_turbine import WindTurbine`
 
 Functions:
 

{detquantlib-3.15.0 → detquantlib-4.0.0}/README.md

@@ -54,6 +54,9 @@ poetry run python docs.py
 
 Classes:
 
+- `Battery`:
+  - `from detquantlib.assets import Battery`
+  - `from detquantlib.assets.battery import Battery`
 - `DetDatabase`:
   - `from detquantlib.data import DetDatabase`
   - `from detquantlib.data.databases.detdatabase import DetDatabase`
@@ -72,6 +75,9 @@ Classes:
 - `Sftp`:
   - `from detquantlib.data import Sftp`
   - `from detquantlib.data.sftp.sftp import Sftp`
+- `WindTurbine`:
+  - `from detquantlib.assets import WindTurbine`
+  - `from detquantlib.assets.wind_turbine import WindTurbine`
 
 Functions:
 

detquantlib-4.0.0/detquantlib/assets/__init__.py

@@ -0,0 +1,4 @@
+from .battery import Battery
+from .wind_turbine import WindTurbine, WindTurbineModels
+
+__all__ = ["Battery", "WindTurbine", "WindTurbineModels"]

detquantlib-4.0.0/detquantlib/assets/wind_turbine.py

@@ -0,0 +1,345 @@
+# Third-party packages
+import numpy as np
+import pandas as pd
+from scipy.interpolate import CubicSpline
+
+# Internal modules
+from detquantlib.converters import power_to_energy
+
+
+class WindTurbine:
+    """A class to store wind turbine information and perform wind turbine-related calculations."""
+
+    def __init__(
+        self,
+        rated_power_mw: float,
+        cut_in: float,
+        cut_out: float,
+        nr_turbines: int,
+        power_curve: CubicSpline = None,
+    ):
+        """
+        Constructor method.
+
+        Args:
+            rated_power_mw: Rated power of a single wind turbine (in MW)
+            cut_in: Cut-in wind speed
+            cut_out: Cut-out wind speed
+            nr_turbines: Number of wind turbines in the wind farm
+            power_curve: Power curve of a single wind turbine, stored as a cubic spline object.
+                If not provided, can be set via the 'fit_power_curve()' method below.
+        """
+        self.rated_power_mw = rated_power_mw
+        self.cut_in = cut_in
+        self.cut_out = cut_out
+        self.nr_turbines = nr_turbines
+        self.power_curve = power_curve
+
+    def fit_power_curve(
+        self, wind_speed: list | np.ndarray | pd.Series, power_mw: list | np.ndarray | pd.Series
+    ):
+        """
+        Fits the wind turbine's power curve with a cubic spline, and stores the resulting spline
+        object in attribute self.power_curve.
+
+        Args:
+            wind_speed: Wind speeds
+            power_mw: Corresponding power in MW (for a single wind turbine)
+        """
+        # Make sure input data is stored as numpy array
+        wind_speed, power_mw = np.array(wind_speed), np.array(power_mw)
+
+        # Trim data to cut-in - cut-out interval
+        idx = (wind_speed >= self.cut_in) & (wind_speed <= self.cut_out)
+        wind_speed_trim, power_trim = wind_speed[idx], power_mw[idx]
+
+        # Add data point with power=0.0 before cut-in to ensure smooth behaviour of spline
+        # interpolation
+        w0 = wind_speed_trim[0] - np.diff(wind_speed_trim[:2])
+        p0 = [0]
+        wind_speed_trim, power_trim = np.concat([w0, wind_speed_trim]), np.concat([p0, power_trim])
+
+        # Fit cubic spline
+        cs = CubicSpline(wind_speed_trim, power_trim)
+
+        # Store spline object
+        self.power_curve = cs
+
+    def wind_to_power(self, wind_speed: list | np.ndarray | pd.Series) -> np.ndarray:
+        """
+        Converts wind speeds to generated power, aggregated over the total number of turbines.
+
+        Args:
+            wind_speed: Wind speeds
+
+        Returns:
+            Corresponding power in MW (total over the number of turbines set in the object)
+
+        Raises:
+            ValueError: Raises an error if self.power_curve is None.
+        """
+        # Check that power curve attribute is defined
+        if self.power_curve is None:
+            raise ValueError(
+                "Cannot execute method self.wind_to_power() if attribute self.power_curve is "
+                "None. Use method self.fit_power_curve() to set self.power_curve based on power "
+                "curve data."
+            )
+
+        # Make sure input data is stored as numpy array
+        wind_speed = np.array(wind_speed)
+
+        # Cubic spline interpolation
+        power_mw = self.power_curve(wind_speed)
+
+        # Cap to rated power
+        power_mw = np.minimum(power_mw, self.rated_power_mw)
+
+        # Adjust values outside of cut-in - cut-out bounds
+        idx = (wind_speed < self.cut_in) | (wind_speed > self.cut_out)
+        power_mw[idx] = 0.0
+
+        # Rescale power to number of wind turbines
+        power_mw *= self.nr_turbines
+
+        return power_mw
+
+    def wind_to_energy(self, wind_speed: list | np.ndarray | pd.Series, tenor: str) -> np.ndarray:
+        """
+        Converts wind speeds to generated energy, aggregated over the total number of turbines.
+
+        Args:
+            wind_speed: Wind speeds
+            tenor: Product tenor
+
+        Returns:
+            Corresponding energy in MWh (total over the number of turbines set in the object)
+        """
+        power_mw = self.wind_to_power(wind_speed)
+        energy_mwh = power_to_energy(
+            power=power_mw, power_unit="MW", energy_unit="MWh", tenor=tenor
+        )
+        return energy_mwh
+
+    @staticmethod
+    def build_from_model(model: dict, nr_turbines: int):
+        """
+        Builds a WindTurbine object from a known wind turbine model.
+
+        Args:
+            model: Wind turbine model. The dictionary can be taken from the list of models in
+                the WindTurbineModels class. Otherwise, the dictionary should have the same
+                structure as those in the WindTurbineModels class.
+            nr_turbines: Number of wind turbines in the wind farm
+
+        Returns:
+            WindTurbine object
+        """
+        # Create object
+        obj = WindTurbine(
+            rated_power_mw=model["rated_power_mw"],
+            cut_in=model["cut_in"],
+            cut_out=model["cut_out"],
+            nr_turbines=nr_turbines,
+        )
+
+        # Set power_curve attribute
+        obj.fit_power_curve(
+            wind_speed=[d["wind_speed_(m/s)"] for d in model["power_curve"]],
+            power_mw=[d["power(mw)"] for d in model["power_curve"]],
+        )
+
+        return obj
+
+
+class WindTurbineModels:
+    """
+    An inventory of known wind turbine models. Each model is stored as an object attribute. All
+    attribute dictionaries must have the same structure.
+    """
+
+    enercon_e126_ep3_4mw = dict(
+        brand="Enercon",
+        model="E-126 EP3 4.0MW",
+        rated_power_mw=4.0,
+        cut_in=2.5,
+        cut_out=25.0,
+        power_curve_source="https://www.thewindpower.net/turbine_en_1576_enercon_e126-4000.php",
+        power_curve=[
+            {"wind_speed_(m/s)": 0.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 0.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 2.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 2.5, "power(mw)": 0.02},
+            {"wind_speed_(m/s)": 3.0, "power(mw)": 0.047},
+            {"wind_speed_(m/s)": 3.5, "power(mw)": 0.095},
+            {"wind_speed_(m/s)": 4.0, "power(mw)": 0.15},
+            {"wind_speed_(m/s)": 4.5, "power(mw)": 0.29},
+            {"wind_speed_(m/s)": 5.0, "power(mw)": 0.441},
+            {"wind_speed_(m/s)": 5.5, "power(mw)": 0.582},
+            {"wind_speed_(m/s)": 6.0, "power(mw)": 0.732},
+            {"wind_speed_(m/s)": 6.5, "power(mw)": 0.945},
+            {"wind_speed_(m/s)": 7.0, "power(mw)": 1.165},
+            {"wind_speed_(m/s)": 7.5, "power(mw)": 1.43},
+            {"wind_speed_(m/s)": 8.0, "power(mw)": 1.7},
+            {"wind_speed_(m/s)": 8.5, "power(mw)": 2.012},
+            {"wind_speed_(m/s)": 9.0, "power(mw)": 2.323},
+            {"wind_speed_(m/s)": 9.5, "power(mw)": 2.677},
+            {"wind_speed_(m/s)": 10.0, "power(mw)": 3.023},
+            {"wind_speed_(m/s)": 10.5, "power(mw)": 3.27},
+            {"wind_speed_(m/s)": 11.0, "power(mw)": 3.504},
+            {"wind_speed_(m/s)": 11.5, "power(mw)": 3.698},
+            {"wind_speed_(m/s)": 12.0, "power(mw)": 3.88},
+            {"wind_speed_(m/s)": 12.5, "power(mw)": 3.917},
+            {"wind_speed_(m/s)": 13.0, "power(mw)": 3.945},
+            {"wind_speed_(m/s)": 13.5, "power(mw)": 3.959},
+            {"wind_speed_(m/s)": 14.0, "power(mw)": 3.969},
+            {"wind_speed_(m/s)": 14.5, "power(mw)": 3.995},
+            {"wind_speed_(m/s)": 15.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 15.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 16.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 16.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 17.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 17.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 18.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 18.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 19.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 19.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 20.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 20.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 21.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 21.5, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 22.0, "power(mw)": 4.0},
+            {"wind_speed_(m/s)": 22.5, "power(mw)": 3.985},
+            {"wind_speed_(m/s)": 23.0, "power(mw)": 3.96},
+            {"wind_speed_(m/s)": 23.5, "power(mw)": 3.96},
+            {"wind_speed_(m/s)": 24.0, "power(mw)": 3.952},
+            {"wind_speed_(m/s)": 24.5, "power(mw)": 3.93},
+            {"wind_speed_(m/s)": 25.0, "power(mw)": 3.898},
+        ],
+    )
+
+    enercon_e53_800 = dict(
+        brand="Enercon",
+        model="E-53 800kW",
+        rated_power_mw=0.81,
+        cut_in=2.0,
+        cut_out=25.0,
+        power_curve_source="https://www.thewindpower.net/turbine_en_4_enercon_e53-800.php",
+        power_curve=[
+            {"wind_speed_(m/s)": 0.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 0.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 2.0, "power(mw)": 0.002},
+            {"wind_speed_(m/s)": 2.5, "power(mw)": 0.008},
+            {"wind_speed_(m/s)": 3.0, "power(mw)": 0.014},
+            {"wind_speed_(m/s)": 3.5, "power(mw)": 0.026},
+            {"wind_speed_(m/s)": 4.0, "power(mw)": 0.038},
+            {"wind_speed_(m/s)": 4.5, "power(mw)": 0.057},
+            {"wind_speed_(m/s)": 5.0, "power(mw)": 0.077},
+            {"wind_speed_(m/s)": 5.5, "power(mw)": 0.109},
+            {"wind_speed_(m/s)": 6.0, "power(mw)": 0.141},
+            {"wind_speed_(m/s)": 6.5, "power(mw)": 0.185},
+            {"wind_speed_(m/s)": 7.0, "power(mw)": 0.228},
+            {"wind_speed_(m/s)": 7.5, "power(mw)": 0.282},
+            {"wind_speed_(m/s)": 8.0, "power(mw)": 0.336},
+            {"wind_speed_(m/s)": 8.5, "power(mw)": 0.408},
+            {"wind_speed_(m/s)": 9.0, "power(mw)": 0.48},
+            {"wind_speed_(m/s)": 9.5, "power(mw)": 0.562},
+            {"wind_speed_(m/s)": 10.0, "power(mw)": 0.645},
+            {"wind_speed_(m/s)": 10.5, "power(mw)": 0.71},
+            {"wind_speed_(m/s)": 11.0, "power(mw)": 0.744},
+            {"wind_speed_(m/s)": 11.5, "power(mw)": 0.765},
+            {"wind_speed_(m/s)": 12.0, "power(mw)": 0.785},
+            {"wind_speed_(m/s)": 12.5, "power(mw)": 0.802},
+            {"wind_speed_(m/s)": 13.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 13.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 14.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 14.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 15.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 15.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 16.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 16.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 17.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 17.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 18.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 18.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 19.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 19.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 20.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 20.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 21.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 21.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 22.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 22.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 23.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 23.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 24.0, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 24.5, "power(mw)": 0.81},
+            {"wind_speed_(m/s)": 25.0, "power(mw)": 0.81},
+        ],
+    )
+
+    vestas_v52_850 = dict(
+        brand="Vestas",
+        model="V52 850kW",
+        rated_power_mw=0.85,
+        cut_in=3.0,
+        cut_out=25.0,
+        power_curve_source="https://www.thewindpower.net/turbine_en_27_vestas_v52-850.php",
+        power_curve=[
+            {"wind_speed_(m/s)": 0.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 0.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 1.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 2.0, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 2.5, "power(mw)": 0.0},
+            {"wind_speed_(m/s)": 3.0, "power(mw)": 0.005},
+            {"wind_speed_(m/s)": 3.5, "power(mw)": 0.015},
+            {"wind_speed_(m/s)": 4.0, "power(mw)": 0.03},
+            {"wind_speed_(m/s)": 4.5, "power(mw)": 0.049},
+            {"wind_speed_(m/s)": 5.0, "power(mw)": 0.067},
+            {"wind_speed_(m/s)": 5.5, "power(mw)": 0.097},
+            {"wind_speed_(m/s)": 6.0, "power(mw)": 0.127},
+            {"wind_speed_(m/s)": 6.5, "power(mw)": 0.162},
+            {"wind_speed_(m/s)": 7.0, "power(mw)": 0.197},
+            {"wind_speed_(m/s)": 7.5, "power(mw)": 0.244},
+            {"wind_speed_(m/s)": 8.0, "power(mw)": 0.29},
+            {"wind_speed_(m/s)": 8.5, "power(mw)": 0.35},
+            {"wind_speed_(m/s)": 9.0, "power(mw)": 0.41},
+            {"wind_speed_(m/s)": 9.5, "power(mw)": 0.475},
+            {"wind_speed_(m/s)": 10.0, "power(mw)": 0.539},
+            {"wind_speed_(m/s)": 10.5, "power(mw)": 0.6},
+            {"wind_speed_(m/s)": 11.0, "power(mw)": 0.66},
+            {"wind_speed_(m/s)": 11.5, "power(mw)": 0.71},
+            {"wind_speed_(m/s)": 12.0, "power(mw)": 0.751},
+            {"wind_speed_(m/s)": 12.5, "power(mw)": 0.784},
+            {"wind_speed_(m/s)": 13.0, "power(mw)": 0.816},
+            {"wind_speed_(m/s)": 13.5, "power(mw)": 0.84},
+            {"wind_speed_(m/s)": 14.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 14.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 15.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 15.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 16.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 16.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 17.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 17.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 18.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 18.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 19.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 19.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 20.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 20.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 21.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 21.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 22.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 22.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 23.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 23.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 24.0, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 24.5, "power(mw)": 0.85},
+            {"wind_speed_(m/s)": 25.0, "power(mw)": 0.85},
+        ],
+    )

{detquantlib-3.15.0 → detquantlib-4.0.0}/detquantlib/outputs/outputs_interface.py

@@ -184,26 +184,19 @@ class OutputSet:
 class PathDefinitions:
     """A class containing hard-coded path definitions."""
 
-    @staticmethod
-    def default_outputs_folder_name():
-        return "Outputs"
-
     @staticmethod
     def outputs_plotly_folder_name():
         return "PlotlyFigures"
 
     @staticmethod
-    def get_outputs_folder_name():
-        # Get name from environment variable if it exists, otherwise use default name
-        folder_name = os.getenv("OUTPUTS_FOLDER")
-        if folder_name is None:
-            folder_name = PathDefinitions.default_outputs_folder_name()
-        return folder_name
+    def default_outputs_folder_dir():
+        return Path.cwd().joinpath("Outputs")
 
     @staticmethod
     def get_outputs_folder_dir():
-        # Base directory of the folder containing all model outputs
-        return Path.cwd().joinpath(PathDefinitions.get_outputs_folder_name())
+        # Get path from environment variable if it exists, otherwise use default path
+        output_dir = os.getenv("OUTPUT_DIR")
+        return Path(output_dir) if output_dir else PathDefinitions.default_outputs_folder_dir()
 
     @staticmethod
     def get_outputs_plotly_folder_dir():

{detquantlib-3.15.0 → detquantlib-4.0.0}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "detquantlib"
-version = "3.15.0"
+version = "4.0.0"
 description = "An internal library containing functions and classes that can be used across Quant models."
 authors = ["Dynamic Energy Trading"]
 readme = "README.md"

detquantlib-3.15.0/detquantlib/assets/__init__.py

@@ -1,4 +0,0 @@
-from .battery import Battery
-from .wind_turbine import WindTurbine
-
-__all__ = ["Battery", "WindTurbine"]

detquantlib-3.15.0/detquantlib/assets/wind_turbine.py

@@ -1,92 +0,0 @@
-# Third-party packages
-import numpy as np
-import pandas as pd
-from scipy.interpolate import CubicSpline
-
-
-class WindTurbine:
-    """A class to store wind turbine information and perform wind turbine-related calculations."""
-
-    def __init__(
-        self, rated_power: float, cut_in: float, cut_out: float, power_curve: CubicSpline = None
-    ):
-        """
-        Constructor method.
-
-        Args:
-            rated_power: Rated power
-            cut_in: Cut-in wind speed
-            cut_out: Cut-out wind speed
-            power_curve: Power curve, stored as a cubic spline object. If not provided, can be
-                set via the 'fit_power_curve()' method below.
-        """
-        self.rated_power = rated_power
-        self.cut_in = cut_in
-        self.cut_out = cut_out
-        self.power_curve = power_curve
-
-    def fit_power_curve(
-        self, wind_speed: list | np.ndarray | pd.Series, power: list | np.ndarray | pd.Series
-    ):
-        """
-        Fits the wind turbine's power curve with a cubic spline, and stores the resulting spline
-        object in attribute self.power_curve.
-
-        Args:
-            wind_speed: Wind speeds
-            power: Corresponding power
-        """
-        # Make sure input data is stored as numpy array
-        wind_speed, power = np.array(wind_speed), np.array(power)
-
-        # Trim data to cut-in - cut-out interval
-        idx = (wind_speed >= self.cut_in) & (wind_speed <= self.cut_out)
-        wind_speed_trim, power_trim = wind_speed[idx], power[idx]
-
-        # Add data point with power=0.0 before cut-in to ensure smooth behaviour of spline
-        # interpolation
-        w0 = wind_speed_trim[0] - np.diff(wind_speed_trim[:2])
-        p0 = [0]
-        wind_speed_trim, power_trim = np.concat([w0, wind_speed_trim]), np.concat([p0, power_trim])
-
-        # Fit cubic spline
-        cs = CubicSpline(wind_speed_trim, power_trim)
-
-        # Store spline object
-        self.power_curve = cs
-
-    def wind_to_power(self, wind_speed: list | np.ndarray | pd.Series) -> np.ndarray:
-        """
-        Converts wind speeds to wind turbine power.
-
-        Args:
-            wind_speed: Wind speeds
-
-        Returns:
-            Corresponding power
-
-        Raises:
-            ValueError: Raises an error if self.power_curve is None.
-        """
-        # Check that power curve attribute is defined
-        if self.power_curve is None:
-            raise ValueError(
-                "Cannot execute method self.wind_to_power() if attribute self.power_curve is "
-                "None. Use method self.fit_power_curve() to set self.power_curve based on power "
-                "curve data."
-            )
-
-        # Make sure input data is stored as numpy array
-        wind_speed = np.array(wind_speed)
-
-        # Cubic spline interpolation
-        power = self.power_curve(wind_speed)
-
-        # Cap to rated power
-        power = np.minimum(power, self.rated_power)
-
-        # Adjust values outside of cut-in - cut-out bounds
-        idx = (wind_speed < self.cut_in) | (wind_speed > self.cut_out)
-        power[idx] = 0.0
-
-        return power