res-wind-up 0.1.5__py3-none-any.whl

res_wind_up-0.1.5.dist-info/LICENSE.txt

@@ -0,0 +1,26 @@
+BSD 3-Clause License
+
+Copyright (c) 2024, RENEWABLE ENERGY SYSTEMS LIMITED, all rights reserved.
+
+Redistribution and use in source and binary forms, with or without modification, are permitted
+provided that the following conditions are met:
+
+* Redistributions of source code must retain the above copyright notice, this list of conditions
+and the following disclaimer.
+
+* Redistributions in binary form must reproduce the above copyright notice, this list of
+conditions and the following disclaimer in the documentation and/or other materials provided
+with the distribution.
+
+* Neither the name of the copyright holder nor the names of its contributors may be used to
+endorse or promote products derived from this software without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
+IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
+AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER
+OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
+ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE
+OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+POSSIBILITY OF SUCH DAMAGE.

res_wind_up-0.1.5.dist-info/METADATA

@@ -0,0 +1,121 @@
+Metadata-Version: 2.1
+Name: res-wind-up
+Version: 0.1.5
+Summary: A tool to assess yield uplift of wind turbines
+Author-email: Alex Clerc <alex.clerc@res-group.com>
+License: BSD 3-Clause License
+        
+        Copyright (c) 2024, RENEWABLE ENERGY SYSTEMS LIMITED, all rights reserved.
+        
+        Redistribution and use in source and binary forms, with or without modification, are permitted
+        provided that the following conditions are met:
+        
+        * Redistributions of source code must retain the above copyright notice, this list of conditions
+        and the following disclaimer.
+        
+        * Redistributions in binary form must reproduce the above copyright notice, this list of
+        conditions and the following disclaimer in the documentation and/or other materials provided
+        with the distribution.
+        
+        * Neither the name of the copyright holder nor the names of its contributors may be used to
+        endorse or promote products derived from this software without specific prior written permission.
+        
+        THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
+        IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
+        AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER
+        OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+        CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+        SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
+        ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE
+        OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+        POSSIBILITY OF SUCH DAMAGE.
+        
+Project-URL: Homepage, https://github.com/resgroup/wind-up
+Classifier: Development Status :: 4 - Beta
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: BSD License
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE.txt
+Requires-Dist: geographiclib
+Requires-Dist: matplotlib
+Requires-Dist: pandas >=2.0.0
+Requires-Dist: pyarrow
+Requires-Dist: pydantic >=2.0.0
+Requires-Dist: python-dotenv
+Requires-Dist: pyyaml
+Requires-Dist: requests
+Requires-Dist: ruptures
+Requires-Dist: scipy
+Requires-Dist: seaborn
+Requires-Dist: tabulate
+Requires-Dist: toml
+Requires-Dist: utm
+Requires-Dist: tqdm
+Provides-Extra: dev
+Requires-Dist: pytest ; extra == 'dev'
+Requires-Dist: coverage ; extra == 'dev'
+Requires-Dist: poethepoet ; extra == 'dev'
+Requires-Dist: types-pyyaml ; extra == 'dev'
+Requires-Dist: types-tabulate ; extra == 'dev'
+Requires-Dist: types-toml ; extra == 'dev'
+Requires-Dist: types-requests ; extra == 'dev'
+Requires-Dist: types-tqdm ; extra == 'dev'
+Requires-Dist: ruff ; extra == 'dev'
+Requires-Dist: mypy ; extra == 'dev'
+Provides-Extra: jupyter
+Requires-Dist: jupyterlab ; extra == 'jupyter'
+Requires-Dist: notebook ; extra == 'jupyter'
+Requires-Dist: ipywidgets ; extra == 'jupyter'
+
+# wind-up
+A tool to assess yield uplift of wind turbines
+
+[![CI](https://github.com/resgroup/wind-up/actions/workflows/CI.yaml/badge.svg)](https://github.com/resgroup/wind-up/actions/workflows/CI.yaml)
+[![Python 3.10](https://img.shields.io/badge/python-≥3.10-blue.svg)](https://www.python.org/downloads/release/python-3100/)
+[![Lint & Format: Ruff](https://img.shields.io/endpoint?url=https://raw.githubusercontent.com/charliermarsh/ruff/main/assets/badge/v1.json)](https://github.com/charliermarsh/ruff)
+[![Typing: mypy](https://img.shields.io/badge/typing-mypy-yellow.svg)](https://github.com/python/mypy)
+[![TaskRunner: poethepoet](https://img.shields.io/badge/poethepoet-enabled-1abc9c.svg)](https://github.com/nat-n/poethepoet)
+
+## Getting Started
+See [`examples`](examples) folder for example analysis using the wind-up package. [`smarteole_example.ipynb`](examples%2Fsmarteole_example.ipynb) is a good place to start.
+
+The wind-up package can be installed in a virtual environment with the following commands:
+```shell
+# create and activate a virtual environment, if needed
+python -m venv .venv
+source .venv/Scripts/activate  # or .venv/bin/activate on linux or ".venv/Scripts/activate" in Windows command prompt
+# install the wind-up package in the virtual environment
+pip install res-wind-up # alternatively clone the repo, navigate to the wind-up folder and run "pip install ."
+```
+Note that the package is named `wind_up` (with an underscore) in Python code. For example to print the version of the installed package use the following code snippet:
+```python
+import wind_up
+print(wind_up.__version__)
+```
+
+## Contributing
+To start making changes fork the repository or make a new branch from `main`. Note `main` is protected; 
+if a commit fails to push and you want to undo it try `git reset origin/main --hard`
+
+After cloning the repository (and creating and activating the virtual environment), use the following commands to install the wind-up package in editable mode with the dev dependencies:
+```shell
+git clone https://github.com/resgroup/wind-up # or your fork of wind-up
+cd wind-up
+# create and activate a virtual environment
+python -m venv .venv
+source .venv/Scripts/activate  # or .venv/bin/activate on linux or ".venv/Scripts/activate" in Windows command prompt
+# install the package in editable mode with the dev dependencies
+pip install -e .[dev] # or .[jupyter,dev] if you want jupyter dependencies as well
+```
+Use `poe all` to run all required pre-push commands (make sure the virtual environment is activated)
+
+## Running tests
+Install dev dependencies and use `poe test` to run unit tests (make sure the virtual environment is activated)
+
+## License
+See [`LICENSE.txt`](LICENSE.txt)
+
+## Contact
+Alex.Clerc@res-group.com

res_wind_up-0.1.5.dist-info/RECORD

@@ -0,0 +1,46 @@
+wind_up/__init__.py,sha256=FHt7YmAQu2Eh93KclR8ifJld-7zUVZasJr4BQifdU6M,78
+wind_up/caching.py,sha256=4YEZwry4JObdK9vPrVVRzlxwsoXF90czC1BWgUACPq8,1186
+wind_up/combine_results.py,sha256=t_EY-6fOi1w58cqzSBS0OkpWrE0y_w2e0wGln4lPXJA,6035
+wind_up/constants.py,sha256=qneOsf2QYitrp2HiUSvzBrk1Vwy98Ue3S2sRoCqlSos,1443
+wind_up/conversions.py,sha256=RBCrU_8WqBUwGnIEHYaJTrefaaW_dkT7H7n5YnCF1iQ,230
+wind_up/detrend.py,sha256=zZESO9rCVW5YiGjPZvBoQ4Ne4gQ-sDbEhwkionlnqVk,13306
+wind_up/interface.py,sha256=-kkcsV4qhGfVuGM6g_6hRcREz_O30Qmzoj5L2AeDS8o,8632
+wind_up/long_term.py,sha256=C_oc7DoqpquoH_ILR0KeqBmTXdP0bodk4TjjEl8gcZI,7943
+wind_up/main_analysis.py,sha256=OmvPYzXhILrHxpbbJi1Xy79xlJG0N7RpgJXoie-gs7g,35985
+wind_up/math_funcs.py,sha256=AxMgW6IdfxGtNdg0d3VUJR8nNfIcbLoECoi716WrjTQ,616
+wind_up/models.py,sha256=5DSUtshgpI_25uF-7TR6aC8sqd3LkhIaRtEK5K_DRHw,20687
+wind_up/northing.py,sha256=eNHZnLivEVcGRK3KPtlpTsaFfOI3NSGVXcnihfn4vfc,7173
+wind_up/northing_utils.py,sha256=eAU6qbOoYMPJ7Qv88N25DGQeNT6ov-3uMzRlgUKmkQM,568
+wind_up/optimize_northing.py,sha256=hP9ncyGDcJ4bm1QNdS2VoVN6MaNUvTlpvZzSpcmeNz0,24852
+wind_up/pp_analysis.py,sha256=_82CfySIFPj2GPF_ZF1NQXbcsORxNYvhM_Aw-f5Vpv4,23149
+wind_up/reanalysis_data.py,sha256=LMWkfjN35y5LlxG3inat0gLtSowwyeGNTG7ozbBlXj4,6907
+wind_up/result_manager.py,sha256=S-pOYOOCgobg5roOGfzLtEy5ZYptZ-NOk-p39GvbH4Y,337
+wind_up/scada_funcs.py,sha256=BwBP_DjNjBVi38Hkkz0swDCZ4F75_KHkqt4vCtsSe6E,20543
+wind_up/scada_power_curve.py,sha256=0w6rsGGAUl-tih2l384Nlb6Wgzq2RXKWfK7a_TZB8G8,3497
+wind_up/smart_data.py,sha256=CoiYH1b333al-80kv6v2k_ZCfzgF_eZAmc6HipIgYzA,8435
+wind_up/waking_state.py,sha256=ebBXs0vdY8A1BZOTkMM0nxIXRmNBJBT-uNuIX3i1_MA,12243
+wind_up/wind_funcs.py,sha256=Tod6_HnTVuv6MC2cKJW7xUtOVaeeEu33NlBapFowtxo,361
+wind_up/windspeed_drift.py,sha256=zE7d8-LmcfvkTUxRgYqC1KhCWjLFcq3AA7xmi2UNOTU,3017
+wind_up/ws_est.py,sha256=TlQuSspdwj7acvmZhRI8jkbeiuJGSHR7heDnSNTNlc4,6726
+wind_up/yaml_loader.py,sha256=MbjAsruWVEZrUf8Ni7DMNN5_X58Dp3Z9MIoL83jD4VM,849
+wind_up/plots/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+wind_up/plots/combine_results_plots.py,sha256=Z4zXnNMZ-fCOg8NBT2wXmGu4MAMdHLBfbsVikdLrQL0,1715
+wind_up/plots/data_coverage_plots.py,sha256=U6JA5VJooYwng-fqzV2FFAcTo6B095Pr7Yr12ShtmH8,5602
+wind_up/plots/detrend_plots.py,sha256=e18j9WUJgjTi3tdHGuhPxNcb9lNvDla0zjXd7ZT79Hw,9845
+wind_up/plots/long_term_plots.py,sha256=z5fI5pqXogLcYmtLOhSxtTqPBTn18ST6Emnue5FPLBk,2999
+wind_up/plots/misc_plots.py,sha256=DNCb763IuIYrfX2O-d4Z8YtKLpFyoTRP-2MWUbpzrSw,6480
+wind_up/plots/northing_plots.py,sha256=Mex60h90hK3G7qah-z7fdqkXJveCoWFGoklTvERQyLA,4600
+wind_up/plots/optimize_northing_plots.py,sha256=_FLZmzuNUPoaUF2avlo649f3Tue21_VoO-TcT-OYs7c,2742
+wind_up/plots/pp_analysis_plots.py,sha256=LFKwPl9mIlq7HPyM243VNRZraBfmrEZULz5JVSfgC20,15872
+wind_up/plots/reanalysis_plots.py,sha256=9EbTNQMPyE9-Uagc-yAi_ukMTELGj-4Kh-foHsXfBhI,3685
+wind_up/plots/scada_funcs_plots.py,sha256=7PFDDd84rYzXnLqt0Fure9_FRtSyBfjknppuRuunctY,15759
+wind_up/plots/scada_power_curve_plots.py,sha256=1YF-jKPe4mUOH4o7_S3XIy-sSvfZVabfOO4hZ9t3L6I,3283
+wind_up/plots/waking_state_plots.py,sha256=xs-7q8lxNN41Djc5jAQ5BJuB7ifw7kl_g6UuCCEzFeM,1442
+wind_up/plots/windspeed_drift_plots.py,sha256=ryHwtmfs3rOYxjAyZ_V4do3_wIZbjefQJUm1g2Qks78,775
+wind_up/plots/ws_est_plots.py,sha256=yntSPCVSEMcgaC_bi0Wqdt-ZWFfa6KplqPGBDRowVaI,4063
+wind_up/plots/yaw_direction_plots.py,sha256=BFF_6-B71kLQQazAOvh1tTdgEpqx9ESxarvt-dD9M7M,7399
+res_wind_up-0.1.5.dist-info/LICENSE.txt,sha256=czey0_ea5L8l0GNuhpnHAox9CmWIWULUzBhaBUrQ0h8,1549
+res_wind_up-0.1.5.dist-info/METADATA,sha256=jB4MzDnPQvUZHbCYOLhRvTUUvQVS9iSHxup0BbfN4Go,6087
+res_wind_up-0.1.5.dist-info/WHEEL,sha256=y4mX-SOX4fYIkonsAGA5N0Oy-8_gI4FXw5HNI1xqvWg,91
+res_wind_up-0.1.5.dist-info/top_level.txt,sha256=WXJ6arh0IlkVqkCpnk2L20gM3VfjVxr4iogAAxzB_5o,8
+res_wind_up-0.1.5.dist-info/RECORD,,

res_wind_up-0.1.5.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (70.2.0)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

res_wind_up-0.1.5.dist-info/top_level.txt

@@ -0,0 +1 @@
+wind_up

wind_up/__init__.py

@@ -0,0 +1,3 @@
+from importlib.metadata import version
+
+__version__ = version(__package__)

wind_up/caching.py

@@ -0,0 +1,35 @@
+import pickle
+from collections.abc import Callable
+from pathlib import Path
+from typing import Any
+
+import pandas as pd
+
+from wind_up.result_manager import result_manager
+
+
+def with_parquet_cache(fp: Path, *, use_cache: bool = True) -> Callable:
+    def wrap(func: Callable[..., pd.DataFrame]) -> Callable[..., pd.DataFrame]:
+        def wrapped_f(*a: Any, **kw: Any) -> pd.DataFrame:  # noqa
+            if not Path(fp).is_file() or not use_cache:
+                func(*a, **kw).to_parquet(fp)
+            return pd.read_parquet(fp)
+
+        return wrapped_f
+
+    return wrap
+
+
+def with_pickle_cache(fp: Path, *, use_cache: bool = True) -> Callable:
+    def wrap(func: Callable[..., Any]) -> Callable[..., Any]:
+        def wrapped_f(*a: Any, **kw: Any) -> Any:  # noqa
+            if not Path(fp).is_file() or not use_cache or Path(fp).stat().st_size == 0:
+                with Path.open(fp, "wb") as f:
+                    pickle.dump(func(*a, **kw), f)
+            with Path.open(fp, "rb") as f:
+                result_manager.warning(f"loading cached pickle {fp}")
+                return pickle.load(f)
+
+        return wrapped_f
+
+    return wrap

wind_up/combine_results.py

@@ -0,0 +1,151 @@
+import itertools
+import logging
+import math
+
+import numpy as np
+import pandas as pd
+from scipy.stats import norm
+
+from wind_up.models import PlotConfig
+from wind_up.plots.combine_results_plots import plot_combine_results
+from wind_up.result_manager import result_manager
+
+logger = logging.getLogger(__name__)
+
+
+def calc_sigma_ref(rdf: pd.DataFrame, ref_list: list[str]) -> float:
+    # calculate the weighted absolute average reference uplift as a lower bound on uncertainty
+    ref_uplifts = rdf.loc[(rdf["test_wtg"].isin(ref_list)), "p50_uplift"]
+    ref_sigmas = rdf.loc[(rdf["test_wtg"].isin(ref_list)), "sigma_test"]
+    weights = 1 / (ref_sigmas**2)
+    return (abs(ref_uplifts) * weights).sum() / weights.sum() if weights.sum() > 0 else np.nan
+
+
+def calc_tdf(trdf: pd.DataFrame, ref_list: list[str], weight_col: str = "unc_weight") -> pd.DataFrame:
+    tdf = trdf.groupby("test_wtg").agg(
+        p50_uplift=pd.NamedAgg(
+            column="uplift_frc",
+            aggfunc=lambda x: (x * trdf.loc[x.index, weight_col]).sum() / trdf.loc[x.index, weight_col].sum(),
+        ),
+        sigma_uncorr=pd.NamedAgg(
+            column="unc_one_sigma_frc",
+            aggfunc=lambda x: np.sqrt(
+                ((x * trdf.loc[x.index, weight_col] / trdf.loc[x.index, weight_col].sum()) ** 2).sum(),
+            ),
+        ),
+        sigma_corr=pd.NamedAgg(
+            column="unc_one_sigma_frc",
+            aggfunc=lambda x: (x * trdf.loc[x.index, weight_col]).sum() / trdf.loc[x.index, weight_col].sum(),
+        ),
+        ref_count=pd.NamedAgg(column="uplift_frc", aggfunc=len),
+        is_ref=pd.NamedAgg(column="test_wtg", aggfunc=lambda x: x.isin(ref_list).any()),
+    )
+    tdf["sigma_test"] = (tdf["sigma_uncorr"] + tdf["sigma_corr"]) / 2
+    tdf = tdf.sort_values(by=["ref_count", "test_wtg"], ascending=[False, True])
+    tdf = tdf.reset_index()
+    sigma_ref = calc_sigma_ref(tdf, ref_list)
+    tdf["sigma"] = tdf["sigma_test"].clip(lower=sigma_ref)
+    tdf["p95_uplift"] = tdf["p50_uplift"] + norm.ppf(0.05) * tdf["sigma"]
+    tdf["p5_uplift"] = tdf["p50_uplift"] + norm.ppf(0.95) * tdf["sigma"]
+    return tdf
+
+
+def choose_best_refs(trdf: pd.DataFrame, ref_list: list[str], min_refs: int = 3) -> list[str]:
+    ref_count = len(ref_list)
+    if ref_count < min_refs:
+        msg = "ref_list must have at least min_refs elements"
+        raise ValueError(msg)
+
+    rdf = calc_tdf(trdf, ref_list)
+    best_sigma_ref = calc_sigma_ref(rdf, ref_list)
+    best_ref_list = ref_list.copy()
+    ref_count -= 1
+    while ref_count >= min_refs:
+        for c in itertools.combinations(ref_list, ref_count):
+            this_ref_list = list(c)
+            this_sigma_ref = calc_sigma_ref(rdf, this_ref_list)
+            if this_sigma_ref < best_sigma_ref:
+                best_sigma_ref = this_sigma_ref
+                best_ref_list = this_ref_list.copy()
+        ref_count -= 1
+    return best_ref_list
+
+
+def combine_results(
+    trdf: pd.DataFrame,
+    *,
+    auto_choose_refs: bool = True,
+    exclude_refs: list[str] | None = None,
+    plot_config: PlotConfig | None = None,
+) -> pd.DataFrame:
+    if exclude_refs is None:
+        exclude_refs = []
+
+    msg = "#" * 78 + "\n# combine results per test turbine\n" + "#" * 78
+    logger.info(msg)
+
+    trdf = trdf.copy()
+
+    if trdf.groupby(["test_wtg", "ref"]).size().max() > 1:
+        msg = "trdf must have no more than one row per test-ref pair"
+        raise ValueError(msg)
+
+    # remove reference predictions of themselves
+    trdf = trdf.loc[trdf["test_wtg"] != trdf["ref"], :]
+
+    if len(exclude_refs) > 0:
+        logger.info(f"excluding refs {exclude_refs}")
+        trdf = trdf.loc[~trdf["test_wtg"].isin(exclude_refs), :]
+        trdf = trdf.loc[~trdf["ref"].isin(exclude_refs), :]
+
+    if (trdf["unc_one_sigma_frc"] <= 0).any() or trdf["unc_one_sigma_frc"].isna().any():
+        msg = "unc_one_sigma_frc must be positive and non-NaN"
+        raise ValueError(msg)
+
+    weight_col = "unc_weight"
+    trdf[weight_col] = 1 / (trdf["unc_one_sigma_frc"] ** 2)
+
+    ref_list = sorted(trdf["ref"].unique())
+
+    min_refs = 3
+    if auto_choose_refs:
+        if len(ref_list) >= min_refs:
+            best_ref_list = choose_best_refs(trdf, ref_list, min_refs=min_refs)
+            refs_to_remove = [x for x in ref_list if x not in best_ref_list]
+            trdf = trdf.loc[~trdf["test_wtg"].isin(refs_to_remove), :]
+            trdf = trdf.loc[~trdf["ref"].isin(refs_to_remove), :]
+            ref_list = sorted(trdf["ref"].unique())
+        else:
+            result_manager.warning(f"len(ref_list) < {min_refs}, skipping auto_choose_refs")
+
+    logger.info(f"ref_list = {ref_list}")
+    tdf = calc_tdf(trdf, ref_list, weight_col)
+
+    # change column order for readability
+    cols = list(tdf.columns)
+    first_cols = ["test_wtg", "p50_uplift", "p95_uplift", "p5_uplift", "sigma"]
+    cols = first_cols + [x for x in cols if x not in first_cols]
+    tdf = tdf[cols]
+
+    if plot_config is not None:
+        plot_combine_results(trdf=trdf, tdf=tdf, plot_cfg=plot_config)
+
+    return tdf
+
+
+def calc_net_uplift(results_per_test_df: pd.DataFrame, *, confidence: float) -> tuple[float, float, float]:
+    if results_per_test_df.groupby("test_wtg").size().max() > 1:
+        msg = "results_per_test_df must have only one row per test turbine"
+        raise ValueError(msg)
+    net_p50 = (results_per_test_df["uplift_frc"] * results_per_test_df["mean_power_pre"]).sum() / results_per_test_df[
+        "mean_power_pre"
+    ].sum()
+    net_unc = (
+        math.sqrt(((results_per_test_df["unc_one_sigma_frc"] * results_per_test_df["mean_power_pre"]) ** 2).sum())
+        / results_per_test_df["mean_power_pre"].sum()
+    )
+    p_low = (1 - confidence) / 2
+    net_p_low = net_p50 + norm.ppf(p_low) * net_unc
+    p_high = 1 - p_low
+    net_p_high = net_p50 + norm.ppf(p_high) * net_unc
+    return net_p50, net_p_low, net_p_high

wind_up/constants.py

@@ -0,0 +1,45 @@
+from pathlib import Path
+
+PROJECTROOT_DIR = Path(__file__).parents[1]
+CONFIG_DIR = Path(__file__).parents[1] / "config"
+TURBINE_DATA_DIR = Path(__file__).parents[1] / "input_data/turbine_data"
+REANALYSIS_DIR = Path(__file__).parents[1] / "input_data/reanalysis"
+TOGGLE_DIR = Path(__file__).parents[1] / "input_data/toggle"
+OUTPUT_DIR = Path(__file__).parents[1] / "output"
+
+RANDOM_SEED = 0
+SCATTER_S = 1
+SCATTER_ALPHA = 0.2
+
+
+class DataColumns:
+    turbine_name = "TurbineName"
+    active_power_mean = "ActivePowerMean"
+    active_power_sd = "ActivePowerSD"
+    wind_speed_mean = "WindSpeedMean"
+    wind_speed_sd = "WindSpeedSD"
+    yaw_angle_mean = "YawAngleMean"
+    yaw_angle_min = "YawAngleMin"
+    yaw_angle_max = "YawAngleMax"
+    pitch_angle_mean = "PitchAngleMean"
+    gen_rpm_mean = "GenRpmMean"
+    ambient_temp = "AmbientTemp"
+    shutdown_duration = "ShutdownDuration"
+
+    @classmethod
+    def all(cls: type["DataColumns"]) -> list[str]:
+        return [v for k, v in vars(cls).items() if not k.startswith("_") and k != "all"]
+
+
+HOURS_PER_YEAR = 8766
+DEFAULT_AIR_DENSITY = 1.22
+
+TIMESTAMP_COL = "TimeStamp_StartFormat"
+RAW_WINDSPEED_COL = "raw_WindSpeedMean"
+RAW_POWER_COL = "raw_ActivePowerMean"
+RAW_DOWNTIME_S_COL = "raw_ShutdownDuration"
+RAW_YAWDIR_COL = "raw_YawAngleMean"
+
+REANALYSIS_WS_COL = "reanalysis_ws"
+REANALYSIS_WD_COL = "reanalysis_wd"
+WINDFARM_YAWDIR_COL = "wf_yawdir"

wind_up/conversions.py

@@ -0,0 +1,11 @@
+from typing import TypeVar
+
+import pandas as pd
+
+T = TypeVar("T", pd.Timestamp, pd.DatetimeIndex)
+
+
+def ensure_utc(t: T) -> T:
+    if t.tzinfo is None:
+        return t.tz_localize("UTC")
+    return t.tz_convert("UTC")