something

py_wake/examples/hornsrev1.py → py_wake/examples/hornsrev1_example.py

 import numpy as np
 from py_wake.aep._aep import AEP
-from py_wake.examples.data.hornsrev_setup import wt_x, wt_y, HornsrevV80,\
+from py_wake.examples.data.hornsrev1 import wt_x, wt_y, HornsrevV80,\
     HornsrevSite
 from py_wake.wake_models.noj import NOJ
 

py_wake/wake_models/fuga.py

@@ -55,7 +55,7 @@ class FugaWakeModel(WakeModel, LinearSum):
         mdu = [np.fromfile(path + prefix + '%04dUL.dat' % j, np.dtype('<f'), -1)
                for j in zlevels]
 
-        self.du = np.array(mdu, dtype=np.float32).reshape((len(mdu), nyW // 2, nxW)) * factor
+        self.du = -np.array(mdu, dtype=np.float32).reshape((len(mdu), nyW // 2, nxW)) * factor
         self.z0 = z0
         self.x0 = nxW // 4
         self.dx = dx
@@ -73,8 +73,9 @@ class FugaWakeModel(WakeModel, LinearSum):
 
     def calc_deficit(self, WS_lk, WS_eff_lk, dw_jl, cw_jl, ct_lk):
         mdu_jl = self.interpolate(dw_jl, cw_jl, 70)
-        deficit_jlk = mdu_jl[:, :, na] * ct_lk[na] * WS_eff_lk**2 / WS_lk
-        return -deficit_jlk
+        deficit_jlk = mdu_jl[:, :, na] * ct_lk[na] * (WS_eff_lk**2 / WS_lk)
+
+        return deficit_jlk
 
 
 def main():

py_wake/wind_turbines/iea37_wind_turbine.py

-from py_wake.examples.data.iea37.iea37_reader import read_iea37_windturbine
-from py_wake.wind_turbines import OneTypeWindTurbines
-from py_wake.examples.data.iea37 import iea37_path
-
-
-class IEA37_WindTurbines(OneTypeWindTurbines):
-    def __init__(self, yaml_filename):
-        name, hub_height, diameter, ct_func, power_func = read_iea37_windturbine(yaml_filename)
-        super().__init__(name, diameter, hub_height, ct_func, power_func)
-
-
-def main():
-    if __name__ == '__main__':
-        wt = IEA37_WindTurbines(iea37_path + 'iea37-335mw.yaml')
-        print(wt.diameter(0))
-        print(wt.hub_height(0))
-
-
-main()