foxes 1.1.1__py3-none-any.whl → 1.2__py3-none-any.whl

tests/1_verification/flappy_0_6/row_Jensen_quadratic_centre/test_row_Jensen_quadratic_centre.py

@@ -5,7 +5,7 @@ import inspect
 
 import foxes
 import foxes.variables as FV
-import foxes.constants as FC
+from foxes.config import config
 
 thisdir = Path(inspect.getfile(inspect.currentframe())).parent
 
@@ -39,54 +39,56 @@ def test():
         verbosity=0,
     )
 
-    algo = foxes.algorithms.Downwind(
-        farm,
-        states,
-        mbook=mbook,
-        rotor_model="centre",
-        wake_models=["Jensen_quadratic_k007"],
-        wake_frame="rotor_wd",
-        partial_wakes={"Jensen_linear_k007": "centre"},
-        verbosity=0,
-    )
-
-    data = algo.calc_farm()
-
-    df = data.to_dataframe()[[FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]]
-
-    print("\nReading file", cfile)
-    fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
-
-    print()
-    print("TRESULTS\n")
-    sel = (df[FV.P] > 0) & (fdata[FV.P] > 0)
-    df = df.loc[sel]
-    fdata = fdata.loc[sel]
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-
-    print("\nVERIFYING\n")
-    df[FV.WS] = df["REWS"]
-    df[FV.AMB_WS] = df["AMB_REWS"]
-
-    delta = df - fdata
-    print(delta)
-
-    chk = delta.abs()
-    print(chk.max())
-
-    var = FV.WS
-    print(f"\nCHECKING {var}")
-    sel = chk[var] >= 1e-7
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-    print(chk.loc[sel])
-    assert (chk[var] < 1e-7).all()
-
-    var = FV.P
-    sel = chk[var] >= 1e-5
-    print(f"\nCHECKING {var}\n", delta.loc[sel])
-    assert (chk[var] < 1e-5).all()
+    with foxes.Engine.new("threads", chunk_size_states=c):
+
+        algo = foxes.algorithms.Downwind(
+            farm,
+            states,
+            mbook=mbook,
+            rotor_model="centre",
+            wake_models=["Jensen_quadratic_k007"],
+            wake_frame="rotor_wd",
+            partial_wakes={"Jensen_linear_k007": "centre"},
+            verbosity=0,
+        )
+
+        data = algo.calc_farm()
+
+        df = data.to_dataframe()[[FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]]
+
+        print("\nReading file", cfile)
+        fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
+
+        print()
+        print("TRESULTS\n")
+        sel = (df[FV.P] > 0) & (fdata[FV.P] > 0)
+        df = df.loc[sel]
+        fdata = fdata.loc[sel]
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+
+        print("\nVERIFYING\n")
+        df[FV.WS] = df["REWS"]
+        df[FV.AMB_WS] = df["AMB_REWS"]
+
+        delta = df - fdata
+        print(delta)
+
+        chk = delta.abs()
+        print(chk.max())
+
+        var = FV.WS
+        print(f"\nCHECKING {var}")
+        sel = chk[var] >= 1e-7
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+        print(chk.loc[sel])
+        assert (chk[var] < 1e-7).all()
+
+        var = FV.P
+        sel = chk[var] >= 1e-5
+        print(f"\nCHECKING {var}\n", delta.loc[sel])
+        assert (chk[var] < 1e-5).all()
 
 
 if __name__ == "__main__":

tests/1_verification/flappy_0_6_2/grid_rotors/test_grid_rotors.py

@@ -4,7 +4,7 @@ import inspect
 
 import foxes
 import foxes.variables as FV
-import foxes.constants as FC
+from foxes.config import config
 
 thisdir = Path(inspect.getfile(inspect.currentframe())).parent
 
@@ -22,108 +22,106 @@ def test():
         (["Basta"], "grid64", "grid64"),
     ]
 
-    ck = {FC.STATE: c}
-
-    for i, (wakes, rotor, pwake) in enumerate(cases):
-        print(f"\nENTERING CASE {(wakes, rotor, pwake)}\n")
-
-        mbook = foxes.models.ModelBook()
-        ttype = foxes.models.turbine_types.PCtFile(
-            data_source=tfile, var_ws_ct=FV.REWS, var_ws_P=FV.REWS
-        )
-        mbook.turbine_types[ttype.name] = ttype
-
-        mbook.wake_models["Basta"] = foxes.models.wake_models.wind.Bastankhah2014(
-            sbeta_factor=0.25, superposition="ws_linear", induction="Betz"
-        )
-
-        states = foxes.input.states.StatesTable(
-            data_source=sfile,
-            output_vars=[FV.WS, FV.WD, FV.TI, FV.RHO],
-            var2col={FV.WS: "ws", FV.WD: "wd", FV.TI: "ti"},
-            fixed_vars={FV.RHO: 1.225, FV.Z0: 0.1, FV.H: 100.0},
-            profiles={FV.WS: "ABLLogNeutralWsProfile"},
-        )
-
-        farm = foxes.WindFarm()
-        foxes.input.farm_layout.add_from_file(
-            farm,
-            lfile,
-            turbine_models=["kTI_amb_02", ttype.name],
-            verbosity=1,
-        )
-
-        algo = foxes.algorithms.Downwind(
-            farm,
-            states,
-            mbook=mbook,
-            rotor_model=rotor,
-            wake_models=wakes,
-            wake_frame="rotor_wd",
-            partial_wakes=pwake,
-            chunks=ck,
-            verbosity=1,
-        )
-
-        data = algo.calc_farm()
-
-        df = data.to_dataframe()[
-            [FV.AMB_WD, FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]
-        ]
-
-        cfile = cpath / f"results_{i}.csv.gz"
-        print("\nReading file", cfile)
-        fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
-
-        print()
-        print("TRESULTS\n")
-        print(df)
-        print(fdata)
-
-        print("\nVERIFYING\n")
-        df[FV.WS] = df["REWS"]
-        df[FV.AMB_WS] = df["AMB_REWS"]
-
-        delta = df - fdata
-        print(delta)
-        chk = delta[[FV.AMB_WS, FV.AMB_P, FV.WS, FV.P]]
-        print(chk)
-        chk = chk.abs()
-        print(chk.max())
-
-        var = FV.AMB_WS
-        sel = chk[var] >= 1e-7
-        print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
-        print(df.loc[sel])
-        print(fdata.loc[sel])
-        print(delta.loc[sel])
-        assert (chk[var] < 1e-7).all()
-
-        var = FV.AMB_P
-        sel = chk[var] >= 1e-5
-        print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
-        print(df.loc[sel])
-        print(fdata.loc[sel])
-        print(delta.loc[sel])
-        assert (chk[var] < 1e-5).all()
-
-        var = FV.WS
-        sel = chk[var] >= 1.7e-3
-        print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
-        print(df.loc[sel])
-        print(fdata.loc[sel])
-        print(delta.loc[sel])
-        assert (chk[var] < 1.7e-3).all()
-
-        var = FV.P
-        sel = chk[var] >= 1.51
-        print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
-        print(df.loc[sel])
-        print(fdata.loc[sel])
-        print(delta.loc[sel])
-        assert (chk[var] < 1.51).all()
-
-        print()
+    with foxes.Engine.new("threads", chunk_size_states=c):
+        for i, (wakes, rotor, pwake) in enumerate(cases):
+            print(f"\nENTERING CASE {(wakes, rotor, pwake)}\n")
+
+            mbook = foxes.models.ModelBook()
+            ttype = foxes.models.turbine_types.PCtFile(
+                data_source=tfile, var_ws_ct=FV.REWS, var_ws_P=FV.REWS
+            )
+            mbook.turbine_types[ttype.name] = ttype
+
+            mbook.wake_models["Basta"] = foxes.models.wake_models.wind.Bastankhah2014(
+                sbeta_factor=0.25, superposition="ws_linear", induction="Betz"
+            )
+
+            states = foxes.input.states.StatesTable(
+                data_source=sfile,
+                output_vars=[FV.WS, FV.WD, FV.TI, FV.RHO],
+                var2col={FV.WS: "ws", FV.WD: "wd", FV.TI: "ti"},
+                fixed_vars={FV.RHO: 1.225, FV.Z0: 0.1, FV.H: 100.0},
+                profiles={FV.WS: "ABLLogNeutralWsProfile"},
+            )
+
+            farm = foxes.WindFarm()
+            foxes.input.farm_layout.add_from_file(
+                farm,
+                lfile,
+                turbine_models=["kTI_amb_02", ttype.name],
+                verbosity=1,
+            )
+
+            algo = foxes.algorithms.Downwind(
+                farm,
+                states,
+                mbook=mbook,
+                rotor_model=rotor,
+                wake_models=wakes,
+                wake_frame="rotor_wd",
+                partial_wakes=pwake,
+                verbosity=1,
+            )
+
+            data = algo.calc_farm()
+
+            df = data.to_dataframe()[
+                [FV.AMB_WD, FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]
+            ]
+
+            cfile = cpath / f"results_{i}.csv.gz"
+            print("\nReading file", cfile)
+            fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
+
+            print()
+            print("TRESULTS\n")
+            print(df)
+            print(fdata)
+
+            print("\nVERIFYING\n")
+            df[FV.WS] = df["REWS"]
+            df[FV.AMB_WS] = df["AMB_REWS"]
+
+            delta = df - fdata
+            print(delta)
+            chk = delta[[FV.AMB_WS, FV.AMB_P, FV.WS, FV.P]]
+            print(chk)
+            chk = chk.abs()
+            print(chk.max())
+
+            var = FV.AMB_WS
+            sel = chk[var] >= 1e-7
+            print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
+            print(df.loc[sel])
+            print(fdata.loc[sel])
+            print(delta.loc[sel])
+            assert (chk[var] < 1e-7).all()
+
+            var = FV.AMB_P
+            sel = chk[var] >= 1e-5
+            print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
+            print(df.loc[sel])
+            print(fdata.loc[sel])
+            print(delta.loc[sel])
+            assert (chk[var] < 1e-5).all()
+
+            var = FV.WS
+            sel = chk[var] >= 1.7e-3
+            print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
+            print(df.loc[sel])
+            print(fdata.loc[sel])
+            print(delta.loc[sel])
+            assert (chk[var] < 1.7e-3).all()
+
+            var = FV.P
+            sel = chk[var] >= 1.51
+            print(f"\nCHECKING {var}, {(wakes, rotor, pwake)}\n")
+            print(df.loc[sel])
+            print(fdata.loc[sel])
+            print(delta.loc[sel])
+            assert (chk[var] < 1.51).all()
+
+            print()
 
 
 if __name__ == "__main__":

tests/1_verification/flappy_0_6_2/row_Bastankhah_Crespo/test_row_Bastankhah_Crespo.py

@@ -5,7 +5,7 @@ import inspect
 
 import foxes
 import foxes.variables as FV
-import foxes.constants as FC
+from foxes.config import config
 
 thisdir = Path(inspect.getfile(inspect.currentframe())).parent
 
@@ -22,8 +22,6 @@ def test():
     cfile = thisdir / "flappy" / "results.csv.gz"
     tfile = thisdir / "NREL-5MW-D126-H90.csv"
 
-    ck = {FC.STATE: c}
-
     mbook = foxes.models.ModelBook()
     ttype = foxes.models.turbine_types.PCtFile(
         data_source=tfile, var_ws_ct=FV.REWS, var_ws_P=FV.REWS
@@ -51,65 +49,66 @@ def test():
         verbosity=1,
     )
 
-    algo = foxes.algorithms.Downwind(
-        farm,
-        states,
-        mbook=mbook,
-        rotor_model=rotor,
-        wake_models=["Basta", "Crespo"],
-        wake_frame="rotor_wd",
-        partial_wakes=["axiwake6", "top_hat"],
-        chunks=ck,
-        verbosity=1,
-    )
-
-    data = algo.calc_farm()
-
-    df = data.to_dataframe()[
-        [FV.X, FV.Y, FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_TI, FV.TI]
-    ]
-
-    print()
-    print("TRESULTS\n")
-    print(df)
-
-    print("\nReading file", cfile)
-    fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
-
-    print()
-    print("TRESULTS\n")
-    # sel   = (df[FV.P]>0) & (fdata[FV.P]>0)
-    # df    = df.loc[sel]
-    # fdata = fdata.loc[sel]
-    print(df)
-    print(fdata)
-
-    print("\nVERIFYING\n")
-    df[FV.WS] = df["REWS"]
-    df[FV.AMB_WS] = df["AMB_REWS"]
-
-    delta = df - fdata
-    print(df)
-    print(delta[[FV.WS, FV.TI]])
-
-    chk = delta.abs()
-    print(chk.max())
-
-    var = FV.WS
-    print(f"\nCHECKING {var}")
-    sel = chk[var] >= 3e-3
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-    print(chk.loc[sel])
-    assert (chk[var] < 3e-3).all()
-
-    var = FV.TI
-    print(f"\nCHECKING {var}")
-    sel = chk[var] >= 3e-4
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-    print(chk.loc[sel])
-    assert (chk[var] < 3e-4).all()
+    with foxes.Engine.new("threads", chunk_size_states=c):
+
+        algo = foxes.algorithms.Downwind(
+            farm,
+            states,
+            mbook=mbook,
+            rotor_model=rotor,
+            wake_models=["Basta", "Crespo"],
+            wake_frame="rotor_wd",
+            partial_wakes=["axiwake6", "top_hat"],
+            verbosity=1,
+        )
+
+        data = algo.calc_farm()
+
+        df = data.to_dataframe()[
+            [FV.X, FV.Y, FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_TI, FV.TI]
+        ]
+
+        print()
+        print("TRESULTS\n")
+        print(df)
+
+        print("\nReading file", cfile)
+        fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
+
+        print()
+        print("TRESULTS\n")
+        # sel   = (df[FV.P]>0) & (fdata[FV.P]>0)
+        # df    = df.loc[sel]
+        # fdata = fdata.loc[sel]
+        print(df)
+        print(fdata)
+
+        print("\nVERIFYING\n")
+        df[FV.WS] = df["REWS"]
+        df[FV.AMB_WS] = df["AMB_REWS"]
+
+        delta = df - fdata
+        print(df)
+        print(delta[[FV.WS, FV.TI]])
+
+        chk = delta.abs()
+        print(chk.max())
+
+        var = FV.WS
+        print(f"\nCHECKING {var}")
+        sel = chk[var] >= 3e-3
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+        print(chk.loc[sel])
+        assert (chk[var] < 3e-3).all()
+
+        var = FV.TI
+        print(f"\nCHECKING {var}")
+        sel = chk[var] >= 3e-4
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+        print(chk.loc[sel])
+        assert (chk[var] < 3e-4).all()
 
 
 if __name__ == "__main__":

tests/1_verification/flappy_0_6_2/row_Bastankhah_linear_centre/test_row_Bastankhah_linear_centre.py

@@ -5,7 +5,7 @@ import inspect
 
 import foxes
 import foxes.variables as FV
-import foxes.constants as FC
+from foxes.config import config
 
 thisdir = Path(inspect.getfile(inspect.currentframe())).parent
 
@@ -22,8 +22,6 @@ def test():
     cfile = thisdir / "flappy" / "results.csv.gz"
     tfile = thisdir / "NREL-5MW-D126-H90.csv"
 
-    ck = {FC.STATE: c}
-
     mbook = foxes.models.ModelBook()
     ttype = foxes.models.turbine_types.PCtFile(
         data_source=tfile, var_ws_ct=FV.REWS, var_ws_P=FV.REWS
@@ -48,52 +46,53 @@ def test():
         verbosity=1,
     )
 
-    algo = foxes.algorithms.Downwind(
-        farm,
-        states,
-        mbook=mbook,
-        rotor_model=rotor,
-        wake_models=["Basta"],
-        wake_frame="rotor_wd",
-        partial_wakes="rotor_points",
-        chunks=ck,
-        verbosity=1,
-    )
-
-    data = algo.calc_farm()
-
-    df = data.to_dataframe()[[FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]]
-
-    print("\nReading file", cfile)
-    fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
-
-    print()
-    print("TRESULTS\n")
-    sel = (df[FV.P] > 0) & (fdata[FV.P] > 0)
-    df = df.loc[sel]
-    fdata = fdata.loc[sel]
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-
-    print("\nVERIFYING\n")
-    df[FV.WS] = df["REWS"]
-    df[FV.AMB_WS] = df["AMB_REWS"]
-
-    delta = df - fdata
-    print(delta)
-
-    chk = delta.abs()
-    print(chk.max())
-
-    var = FV.WS
-    print(f"\nCHECKING {var}")
-    sel = chk[var] >= 1e-7
-    print(df.loc[sel])
-    print(fdata.loc[sel])
-    print(chk.loc[sel])
-    assert (chk[var] < 1e-7).all()
-
-    var = FV.P
-    sel = chk[var] >= 1e-5
-    print(f"\nCHECKING {var}\n", delta.loc[sel])
-    assert (chk[var] < 1e-5).all()
+    with foxes.Engine.new("threads", chunk_size_states=c):
+
+        algo = foxes.algorithms.Downwind(
+            farm,
+            states,
+            mbook=mbook,
+            rotor_model=rotor,
+            wake_models=["Basta"],
+            wake_frame="rotor_wd",
+            partial_wakes="rotor_points",
+            verbosity=1,
+        )
+
+        data = algo.calc_farm()
+
+        df = data.to_dataframe()[[FV.WD, FV.AMB_REWS, FV.REWS, FV.AMB_P, FV.P]]
+
+        print("\nReading file", cfile)
+        fdata = pd.read_csv(cfile).set_index(["state", "turbine"])
+
+        print()
+        print("TRESULTS\n")
+        sel = (df[FV.P] > 0) & (fdata[FV.P] > 0)
+        df = df.loc[sel]
+        fdata = fdata.loc[sel]
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+
+        print("\nVERIFYING\n")
+        df[FV.WS] = df["REWS"]
+        df[FV.AMB_WS] = df["AMB_REWS"]
+
+        delta = df - fdata
+        print(delta)
+
+        chk = delta.abs()
+        print(chk.max())
+
+        var = FV.WS
+        print(f"\nCHECKING {var}")
+        sel = chk[var] >= 1e-7
+        print(df.loc[sel])
+        print(fdata.loc[sel])
+        print(chk.loc[sel])
+        assert (chk[var] < 1e-7).all()
+
+        var = FV.P
+        sel = chk[var] >= 1e-5
+        print(f"\nCHECKING {var}\n", delta.loc[sel])
+        assert (chk[var] < 1e-5).all()

examples/windio/run.py

@@ -1,29 +0,0 @@
-import argparse
-
-import foxes
-
-if __name__ == "__main__":
-    # define arguments and options:
-    parser = argparse.ArgumentParser()
-    parser.add_argument(
-        "-wio",
-        "--windio_yaml",
-        help="The windio wind energy systems yaml file",
-        default="windio_5turbines_timeseries.yaml",
-    )
-    parser.add_argument(
-        "-nf", "--nofig", help="Do not show figures", action="store_true"
-    )
-    parser.add_argument(
-        "-V", "--verbosity", help="The verbosity level, 0=silent", type=int, default=1
-    )
-    args = parser.parse_args()
-
-    try:
-        wio_runner = foxes.input.windio.read_windio(
-            args.windio_yaml, verbosity=args.verbosity
-        )
-        with wio_runner as runner:
-            runner.run()
-    except ModuleNotFoundError as e:
-        print(e.msg)

foxes/data/windio/DTU_10MW_turbine.yaml

@@ -1,10 +0,0 @@
-name: DTU 10MW Offshore Reference Turbine
-performance:
-  power_curve:
-    power_values: [263388., 751154., 1440738., 2355734., 3506858., 4993092., 6849310., 9116402., 10000754., 10009590., 10000942., 10042678., 10003480., 10001600., 10001506., 10013632., 10007428., 10005360., 10002728., 10001130., 10004984., 9997558.]
-    power_wind_speeds: [4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.]
-  Ct_curve:
-    Ct_values: [0.923,0.919,0.904,0.858,0.814,0.814,0.814,0.814,0.577,0.419,0.323,0.259,0.211,0.175,0.148,0.126,0.109,0.095,0.084,0.074,0.066,0.059]
-    Ct_wind_speeds: [4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.]
-hub_height: 119.0
-rotor_diameter: 178.3