foxes 1.3__tar.gz → 1.5__tar.gz

{foxes-1.3/foxes.egg-info → foxes-1.5}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: foxes
-Version: 1.3
+Version: 1.5
 Summary: Farm Optimization and eXtended yield Evaluation Software
 Author: Jonas Schulte
 Maintainer: Jonas Schulte
@@ -36,7 +36,6 @@ Classifier: Topic :: Scientific/Engineering
 Classifier: Intended Audience :: Developers
 Classifier: Intended Audience :: Science/Research
 Classifier: Programming Language :: Python :: 3
-Classifier: Programming Language :: Python :: 3.8
 Classifier: Programming Language :: Python :: 3.9
 Classifier: Programming Language :: Python :: 3.10
 Classifier: Programming Language :: Python :: 3.11
@@ -45,63 +44,42 @@ Classifier: Programming Language :: Python :: 3.13
 Classifier: License :: OSI Approved :: MIT License
 Classifier: Operating System :: OS Independent
 Classifier: Development Status :: 4 - Beta
-Requires-Python: >=3.8
+Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: matplotlib
-Requires-Dist: numpy
-Requires-Dist: pandas
-Requires-Dist: xarray
-Requires-Dist: scipy
-Requires-Dist: cycler
-Requires-Dist: tqdm
-Requires-Dist: pyyaml
+Requires-Dist: cycler>=0.10
+Requires-Dist: h5netcdf>=1.0
+Requires-Dist: matplotlib>=3.8
+Requires-Dist: numpy>=1.26
+Requires-Dist: pandas>=2.0
+Requires-Dist: pyyaml>=4.0
+Requires-Dist: scipy>=1.12
+Requires-Dist: tqdm>=2.0
+Requires-Dist: xarray>=2023
 Provides-Extra: opt
-Requires-Dist: foxes-opt; extra == "opt"
+Requires-Dist: foxes-opt>=0.5; extra == "opt"
 Provides-Extra: dask
-Requires-Dist: h5netcdf; extra == "dask"
-Requires-Dist: dask; extra == "dask"
-Requires-Dist: distributed; extra == "dask"
-Requires-Dist: dask-jobqueue; extra == "dask"
-Requires-Dist: setuptools; extra == "dask"
-Provides-Extra: eng
-Requires-Dist: h5netcdf; extra == "eng"
-Requires-Dist: multiprocess; extra == "eng"
-Requires-Dist: dask; extra == "eng"
-Requires-Dist: distributed; extra == "eng"
-Requires-Dist: dask-jobqueue; extra == "eng"
-Requires-Dist: setuptools; extra == "eng"
-Requires-Dist: mpi4py; extra == "eng"
-Requires-Dist: ray; extra == "eng"
-Provides-Extra: eng0
-Requires-Dist: h5netcdf; extra == "eng0"
-Requires-Dist: multiprocess; extra == "eng0"
-Requires-Dist: dask; extra == "eng0"
-Requires-Dist: distributed; extra == "eng0"
-Requires-Dist: dask-jobqueue; extra == "eng0"
-Requires-Dist: setuptools; extra == "eng0"
-Requires-Dist: ray; extra == "eng0"
-Provides-Extra: test
-Requires-Dist: flake8; extra == "test"
-Requires-Dist: pytest; extra == "test"
-Requires-Dist: h5netcdf; extra == "test"
+Requires-Dist: dask>=2022.0; extra == "dask"
+Requires-Dist: distributed>=2022.0; extra == "dask"
+Requires-Dist: dask-jobqueue>=0.8; extra == "dask"
+Requires-Dist: setuptools>=61.0; extra == "dask"
 Provides-Extra: doc
-Requires-Dist: setuptools; extra == "doc"
-Requires-Dist: sphinx; extra == "doc"
-Requires-Dist: sphinx-immaterial; extra == "doc"
-Requires-Dist: nbsphinx; extra == "doc"
-Requires-Dist: ipykernel; extra == "doc"
-Requires-Dist: ipywidgets; extra == "doc"
-Requires-Dist: m2r2; extra == "doc"
-Requires-Dist: lxml_html_clean; extra == "doc"
-Requires-Dist: dask; extra == "doc"
-Requires-Dist: distributed; extra == "doc"
-Provides-Extra: dev
-Requires-Dist: flake8; extra == "dev"
-Requires-Dist: pytest; extra == "dev"
-Requires-Dist: jupyter; extra == "dev"
-Requires-Dist: objsize; extra == "dev"
-Requires-Dist: black[jupyter]; extra == "dev"
+Requires-Dist: setuptools>=61.0; extra == "doc"
+Requires-Dist: sphinx>=5.0; extra == "doc"
+Requires-Dist: sphinx-immaterial>=0.10; extra == "doc"
+Requires-Dist: nbsphinx>=0.5; extra == "doc"
+Requires-Dist: ipykernel>=5.0; extra == "doc"
+Requires-Dist: ipywidgets>=5.0; extra == "doc"
+Requires-Dist: m2r2>=0.2; extra == "doc"
+Requires-Dist: lxml_html_clean>=0.4; extra == "doc"
+Requires-Dist: dask>=2022.0; extra == "doc"
+Requires-Dist: distributed>=2022.0; extra == "doc"
+Provides-Extra: test
+Requires-Dist: flake8>=0.1; extra == "test"
+Requires-Dist: pytest>=7.0; extra == "test"
+Provides-Extra: utils
+Requires-Dist: objsize>=0.5; extra == "utils"
+Dynamic: license-file
 
 # Welcome to foxes
 
@@ -161,14 +139,7 @@ Evaluation Software"`
 
 ## Requirements
 
-The supported Python versions are: 
-
-- `Python 3.8`
-- `Python 3.9`
-- `Python 3.10`
-- `Python 3.11`
-- `Python 3.12`
-- `Python 3.13`
+The supported Python versions are `Python 3.9`...`3.13`.
 
 ## Installation
 

{foxes-1.3 → foxes-1.5}/README.md

@@ -56,14 +56,7 @@ Evaluation Software"`
 
 ## Requirements
 
-The supported Python versions are: 
-
-- `Python 3.8`
-- `Python 3.9`
-- `Python 3.10`
-- `Python 3.11`
-- `Python 3.12`
-- `Python 3.13`
+The supported Python versions are `Python 3.9`...`3.13`.
 
 ## Installation
 

{foxes-1.3 → foxes-1.5}/docs/source/conf.py

@@ -22,7 +22,7 @@ from foxes import __version__
 # -- Project information -----------------------------------------------------
 
 project = "foxes"
-copyright = "2024, Fraunhofer IWES"
+copyright = "2025, Fraunhofer IWES"
 author = "Fraunhofer IWES"
 
 # The short X.Y version
@@ -67,9 +67,8 @@ intersphinx_mapping = {
 # autosummary_generate = False
 
 # The suffix(es) of source filenames.
-# You can specify multiple suffix as a list of string:
 #
-source_suffix = [".rst", ".md"]
+source_suffix = {".rst": "restructuredtext", ".md": "restructuredtext"}
 
 # The master toctree document.
 master_doc = "index"
@@ -337,6 +336,7 @@ python_apigen_modules = {
     "foxes.models.vertical_profiles": "_foxes/models/vertical_profiles/",
     "foxes.models.axial_induction": "_foxes/models/axial_induction",
     "foxes.models.ground_models": "_foxes/models/ground_models",
+    "foxes.models.wake_deflections": "_foxes/models/wake_deflections",
     "foxes.models.wake_frames": "_foxes/models/wake_frames/",
     "foxes.models.wake_models": "_foxes/models/wake_models/",
     "foxes.models.wake_models.induction": "_foxes/models/wake_models/induction/",

{foxes-1.3 → foxes-1.5}/examples/abl_states/run.py

@@ -154,7 +154,7 @@ if __name__ == "__main__":
     # power results
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
     print(f"Farm efficiency   : {o.calc_farm_efficiency():.2f}")
     print(f"Annual farm yield : {turbine_results[FV.YLD].sum():.2f} GWh.")

{foxes-1.3 → foxes-1.5}/examples/compare_rotors_pwakes/run.py

@@ -203,7 +203,7 @@ if __name__ == "__main__":
 
         else:
             raise ValueError(
-                f"Please either give one rotor, or one pwake, or same number of both"
+                "Please either give one rotor, or one pwake, or same number of both"
             )
 
     if args.title is not None:

{foxes-1.3 → foxes-1.5}/examples/compare_wakes/run.py

@@ -7,7 +7,6 @@ import foxes.variables as FV
 
 
 def calc(mbook, farm, states, wakes, points, args):
-
     algo = foxes.algorithms.Downwind(
         farm,
         states,
@@ -232,7 +231,7 @@ if __name__ == "__main__":
                 ax.plot(xlist / D, results[0], label=wake)
                 ax.set_ylabel(args.var)
 
-            ax.set_title(f"y = 0")
+            ax.set_title("y = 0")
             ax.set_xlabel("x/D")
             ax.legend(loc="best")
             ax.grid()

{foxes-1.3 → foxes-1.5}/examples/dyn_wakes/run.py

@@ -73,7 +73,6 @@ if __name__ == "__main__":
     parser.add_argument(
         "-p", "--pwakes", help="The partial wakes models", default=None, nargs="+"
     )
-    parser.add_argument("-sc", "--scheduler", help="The scheduler choice", default=None)
     parser.add_argument(
         "-w",
         "--wakes",
@@ -82,6 +81,16 @@ if __name__ == "__main__":
         nargs="+",
     )
     parser.add_argument("-f", "--frame", help="The wake frame", default="dyn_wakes")
+    parser.add_argument(
+        "-d", "--deflection", help="The wake deflection", default="no_deflection"
+    )
+    parser.add_argument(
+        "-y",
+        "--yawm",
+        help="The uniform yaw misalignment value",
+        type=float,
+        default=None,
+    )
     parser.add_argument(
         "-m", "--tmodels", help="The turbine models", default=[], nargs="+"
     )
@@ -121,6 +130,19 @@ if __name__ == "__main__":
         max_length_km=8
     )
 
+    # optionally set turbines in yaw:
+    N = int(args.n_turbines**0.5)
+    if args.yawm is None:
+        ymodels = []
+    else:
+        yawm = np.zeros((1, N * N), dtype=np.float64)
+        yawm[:, :N] = args.yawm
+        mbook.turbine_models["set_yawm"] = foxes.models.turbine_models.SetFarmVars(
+            pre_rotor=True
+        )
+        mbook.turbine_models["set_yawm"].add_var(FV.YAWM, yawm)
+        ymodels = ["set_yawm"]
+
     if args.background0:
         States = foxes.input.states.Timeseries
         kwargs = {}
@@ -147,13 +169,12 @@ if __name__ == "__main__":
     )
 
     farm = foxes.WindFarm()
-    N = int(args.n_turbines**0.5)
     foxes.input.farm_layout.add_grid(
         farm,
         xy_base=np.array([0.0, 0.0]),
         step_vectors=np.array([[1000.0, 0], [0, 800.0]]),
         steps=(N, N),
-        turbine_models=args.tmodels + [ttype.name],
+        turbine_models=ymodels + args.tmodels + [ttype.name],
     )
 
     if not args.nofig and args.show_layout:
@@ -175,6 +196,7 @@ if __name__ == "__main__":
         rotor_model=args.rotor,
         wake_models=args.wakes,
         wake_frame=args.frame,
+        wake_deflection=args.deflection,
         partial_wakes=args.pwakes,
         mbook=mbook,
         max_it=args.max_it,
@@ -221,9 +243,9 @@ if __name__ == "__main__":
     # power results
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
-    print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+    print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
 
     engine.finalize()
 
@@ -262,6 +284,7 @@ if __name__ == "__main__":
             title=None,
             animated=True,
             precalc=True,
+            rotor_color="red",
         )
         next(ofg)
 

{foxes-1.3 → foxes-1.5}/examples/field_data_nc/run.py

@@ -70,7 +70,7 @@ if __name__ == "__main__":
     )
     args = parser.parse_args()
 
-    states = foxes.input.states.FieldDataNC(
+    states = foxes.input.states.FieldData(
         args.file_pattern,
         output_vars=[FV.WS, FV.WD, FV.TI, FV.RHO],
         load_mode=args.load_mode,

{foxes-1.3 → foxes-1.5}/examples/induction/run.py

@@ -147,9 +147,9 @@ if __name__ == "__main__":
     # power results
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
-    print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+    print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
     print(f"Annual farm yield : {turbine_results[FV.YLD].sum():.2f} GWh.")
 
     # horizontal flow plot

{foxes-1.3 → foxes-1.5}/examples/multi_height/run.py

@@ -110,4 +110,4 @@ if __name__ == "__main__":
     o = foxes.output.FarmResultsEval(farm_results)
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power: {P/1000:.1f} MW, Efficiency = {P/P0*100:.2f} %")
+    print(f"\nFarm power: {P / 1000:.1f} MW, Efficiency = {P / P0 * 100:.2f} %")

{foxes-1.3 → foxes-1.5}/examples/power_mask/run.py

@@ -162,7 +162,7 @@ if __name__ == "__main__":
         o = foxes.output.FarmResultsEval(farm_results)
         P0 = o.calc_mean_farm_power(ambient=True)
         P = o.calc_mean_farm_power()
-        print(f"\nFarm power: {P/1000:.1f} MW, Efficiency = {P/P0*100:.2f} %")
+        print(f"\nFarm power: {P / 1000:.1f} MW, Efficiency = {P / P0 * 100:.2f} %")
 
         o1 = foxes.output.StateTurbineMap(farm_results)
 
@@ -180,7 +180,7 @@ if __name__ == "__main__":
         o = foxes.output.FarmResultsEval(farm_results)
         P0 = o.calc_mean_farm_power(ambient=True)
         P = o.calc_mean_farm_power()
-        print(f"\nFarm power: {P/1000:.1f} MW, Efficiency = {P/P0*100:.2f} %")
+        print(f"\nFarm power: {P / 1000:.1f} MW, Efficiency = {P / P0 * 100:.2f} %")
 
         if not args.nofig:
             o0 = foxes.output.StateTurbineMap(farm_results)

{foxes-1.3 → foxes-1.5}/examples/quickstart/run.py

@@ -1,7 +1,6 @@
 import foxes
 
 if __name__ == "__main__":
-
     states = foxes.input.states.Timeseries(
         "timeseries_3000.csv.gz", ["WS", "WD", "TI", "RHO"]
     )

{foxes-1.3 → foxes-1.5}/examples/random_timeseries/run.py

@@ -99,7 +99,6 @@ if __name__ == "__main__":
         chunk_size_states=args.chunksize_states,
         chunk_size_points=args.chunksize_points,
     ):
-
         if not args.nofig:
             fig = plt.figure(figsize=(14.5, 7))
             ax1 = fig.add_subplot(121)
@@ -204,7 +203,7 @@ if __name__ == "__main__":
         # power results
         P0 = o.calc_mean_farm_power(ambient=True)
         P = o.calc_mean_farm_power()
-        print(f"\nFarm power        : {P/1000:.1f} MW")
-        print(f"Farm ambient power: {P0/1000:.1f} MW")
-        print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+        print(f"\nFarm power        : {P / 1000:.1f} MW")
+        print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+        print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
         print(f"Annual farm yield : {turbine_results[FV.YLD].sum():.2f} GWh")

{foxes-1.3 → foxes-1.5}/examples/scan_row/run.py

@@ -161,9 +161,9 @@ if __name__ == "__main__":
     # power results
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
-    print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+    print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
     print(f"Annual farm yield : {turbine_results[FV.YLD].sum():.2f} GWh.")
     print()
 

{foxes-1.3 → foxes-1.5}/examples/sequential/run.py

@@ -24,9 +24,7 @@ if __name__ == "__main__":
         type=int,
         default=4,
     )
-    parser.add_argument(
-        "-d", "--debug", help="Switch on wake debugging", action="store_true"
-    )
+    parser.add_argument("--debug", help="Switch on wake debugging", action="store_true")
     parser.add_argument(
         "-S",
         "--n_states",
@@ -83,6 +81,16 @@ if __name__ == "__main__":
         type=float,
     )
     parser.add_argument("-e", "--engine", help="The engine", default="NumpyEngine")
+    parser.add_argument(
+        "-d", "--deflection", help="The wake deflection", default="no_deflection"
+    )
+    parser.add_argument(
+        "-y",
+        "--yawm",
+        help="The uniform yaw misalignment value",
+        type=float,
+        default=None,
+    )
     parser.add_argument(
         "-n", "--n_cpus", help="The number of cpus", default=None, type=int
     )
@@ -109,21 +117,34 @@ if __name__ == "__main__":
     ttype = foxes.models.turbine_types.PCtFile(args.turbine_file)
     mbook.turbine_types[ttype.name] = ttype
 
+    # optionally set turbines in yaw:
+    N = 3
+    if args.yawm is None:
+        ymodels = []
+    else:
+        yawm = np.zeros((args.n_states, N * N), dtype=np.float64)
+        yawm[:, :N] = args.yawm
+        mbook.turbine_models["set_yawm"] = foxes.models.turbine_models.SetFarmVars(
+            pre_rotor=True
+        )
+        mbook.turbine_models["set_yawm"].add_var(FV.YAWM, yawm)
+        ymodels = ["set_yawm"]
+
     states = foxes.input.states.Timeseries(
         data_source="timeseries_3000.csv.gz",
         output_vars=[FV.WS, FV.WD, FV.TI, FV.RHO],
-        var2col={FV.WS: "WS", FV.WD: "WD", FV.TI: "TI", FV.RHO: "RHO"},
+        # fixed_vars={FV.WD: 270},
+        var2col={FV.WD: "WD", FV.WS: "WS", FV.TI: "TI", FV.RHO: "RHO"},
         states_sel=range(240, 240 + args.n_states),
     )
 
     farm = foxes.WindFarm()
-    N = 3
     foxes.input.farm_layout.add_grid(
         farm,
         xy_base=np.array([0.0, 0.0]),
         step_vectors=np.array([[1000.0, 0], [0, 800.0]]),
         steps=(N, N),
-        turbine_models=args.tmodels + [ttype.name],
+        turbine_models=ymodels + args.tmodels + [ttype.name],
     )
 
     if not args.nofig and args.show_layout:
@@ -143,6 +164,7 @@ if __name__ == "__main__":
             rotor_model=args.rotor,
             wake_models=args.wakes,
             wake_frame=args.frame,
+            wake_deflection=args.deflection,
             partial_wakes=args.pwakes,
             mbook=mbook,
             verbosity=0,
@@ -157,14 +179,15 @@ if __name__ == "__main__":
                 levels=None,
                 quiver_pars=dict(scale=0.01),
                 quiver_n=307,
-                xmin=-5000,
-                ymin=-5000,
-                xmax=7000,
-                ymax=7000,
+                xmin=-1000,
+                ymin=-1000,
+                xmax=4000,
+                ymax=3000,
                 vmin=0,
                 vmax=args.max_variable,
                 title=None,
                 animated=True,
+                rotor_color="red",
             )
             algo.plugins.append(anigen)
 

{foxes-1.3 → foxes-1.5}/examples/single_state/run.py

@@ -172,13 +172,12 @@ if __name__ == "__main__":
     # power results
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
-    print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+    print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
     print(f"Annual farm yield : {turbine_results[FV.YLD].sum():.2f} GWh.")
 
     if not args.nofig:
-
         # horizontal flow plot
         o = foxes.output.FlowPlots2D(algo, farm_results)
         g = o.gen_states_fig_xy(

{foxes-1.3 → foxes-1.5}/examples/states_lookup_table/run.py

@@ -131,7 +131,7 @@ if __name__ == "__main__":
     o = foxes.output.FarmResultsEval(farm_results)
     P0 = o.calc_mean_farm_power(ambient=True)
     P = o.calc_mean_farm_power()
-    print(f"\nFarm power        : {P/1000:.1f} MW")
-    print(f"Farm ambient power: {P0/1000:.1f} MW")
-    print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+    print(f"\nFarm power        : {P / 1000:.1f} MW")
+    print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+    print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
     print(f"Annual farm yield : {o.calc_farm_yield(algo=algo):.2f} GWh")

{foxes-1.3 → foxes-1.5}/examples/streamline_wakes/run.py

@@ -25,7 +25,16 @@ if __name__ == "__main__":
     parser.add_argument(
         "-p", "--pwakes", help="The partial wakes models", default="centre", nargs="+"
     )
-    parser.add_argument("-sc", "--scheduler", help="The scheduler choice", default=None)
+    parser.add_argument(
+        "-d", "--deflection", help="The wake deflection", default="no_deflection"
+    )
+    parser.add_argument(
+        "-y",
+        "--yawm",
+        help="The uniform yaw misalignment value",
+        type=float,
+        default=None,
+    )
     parser.add_argument(
         "-w",
         "--wakes",
@@ -45,7 +54,7 @@ if __name__ == "__main__":
     parser.add_argument(
         "-lm",
         "--load_mode",
-        help="Do load mode",
+        help="Dataset load mode",
         default="preload",
     )
     parser.add_argument("-e", "--engine", help="The engine", default="process")
@@ -75,7 +84,20 @@ if __name__ == "__main__":
     ttype = foxes.models.turbine_types.PCtFile(args.turbine_file)
     mbook.turbine_types[ttype.name] = ttype
 
-    states = foxes.input.states.FieldDataNC(
+    # optionally set turbines in yaw:
+    N = int(args.n_turbines**0.5)
+    if args.yawm is None:
+        ymodels = []
+    else:
+        yawm = np.zeros((1, N * N), dtype=np.float64)
+        yawm[:, :N] = args.yawm
+        mbook.turbine_models["set_yawm"] = foxes.models.turbine_models.SetFarmVars(
+            pre_rotor=True
+        )
+        mbook.turbine_models["set_yawm"].add_var(FV.YAWM, yawm)
+        ymodels = ["set_yawm"]
+
+    states = foxes.input.states.FieldData(
         args.file_pattern,
         states_coord="state",
         x_coord="x",
@@ -86,17 +108,16 @@ if __name__ == "__main__":
         var2ncvar={FV.WS: "ws", FV.WD: "wd"},
         fixed_vars={FV.RHO: 1.225, FV.TI: 0.1},
         load_mode=args.load_mode,
-        bounds_error=False,
+        interpn_pars=dict(bounds_error=False),
     )
 
     farm = foxes.WindFarm()
-    N = int(args.n_turbines**0.5)
     foxes.input.farm_layout.add_grid(
         farm,
         xy_base=np.array([500.0, 500.0]),
         step_vectors=np.array([[500.0, 0], [0, 500.0]]),
         steps=(N, N),
-        turbine_models=args.tmodels + [ttype.name],
+        turbine_models=ymodels + args.tmodels + [ttype.name],
     )
 
     algo = foxes.algorithms.Downwind(
@@ -105,6 +126,7 @@ if __name__ == "__main__":
         rotor_model=args.rotor,
         wake_models=args.wakes,
         wake_frame=args.wake_frame,
+        wake_deflection=args.deflection,
         partial_wakes=args.pwakes,
         mbook=mbook,
         engine=args.engine,
@@ -136,6 +158,7 @@ if __name__ == "__main__":
             xmax=2500,
             ymin=0,
             ymax=2500,
+            rotor_color="red",
         ):
             plt.show()
             plt.close(fig)

{foxes-1.3 → foxes-1.5}/examples/tab_file/run.py

@@ -131,9 +131,9 @@ if __name__ == "__main__":
         o = foxes.output.FarmResultsEval(farm_results)
         P0 = o.calc_mean_farm_power(ambient=True)
         P = o.calc_mean_farm_power()
-        print(f"\nFarm power        : {P/1000:.1f} MW")
-        print(f"Farm ambient power: {P0/1000:.1f} MW")
-        print(f"Farm efficiency   : {o.calc_farm_efficiency()*100:.2f} %")
+        print(f"\nFarm power        : {P / 1000:.1f} MW")
+        print(f"Farm ambient power: {P0 / 1000:.1f} MW")
+        print(f"Farm efficiency   : {o.calc_farm_efficiency() * 100:.2f} %")
         print(f"Annual farm yield : {o.calc_farm_yield(algo=algo):.2f} GWh")
 
         if not args.nofig and args.calc_mean: