Update grid.py

Author: MYMahfouz

floris/simulation/grid.py

@@ -355,16 +355,28 @@ def set_grid(self) -> None:
         Then, create the grid based on this wind-from-left orientation
         """
         # These are the rotated coordinates of the wind turbines based on the wind direction
-        x, y, z = rotate_coordinates_rel_west(self.wind_directions, self.turbine_coordinates_array)
+        fix_orientation = True #dh. TODO add variable in class
+        if fix_orientation : wd = np.ones_like(self.wind_directions)*270 #dh. do not rotate
+        else : wd = self.wind_directions #dh. do rotate
+            
+        x, y, z = rotate_coordinates_rel_west(wd, self.turbine_coordinates_array)
 
         max_diameter = np.max(self.reference_turbine_diameter)
 
         if self.normal_vector == "z":  # Rules of thumb for horizontal plane
             if self.x1_bounds is None:
-                self.x1_bounds = (np.min(x) - 2 * max_diameter, np.max(x) + 10 * max_diameter)
+                # dh. broaden the flowfiled_planar
+                if fix_orientation : 
+                    self.x1_bounds = (np.min(x) - 10 * max_diameter, np.max(x) + 10 * max_diameter)
+                else : 
+                    self.x1_bounds = (np.min(x) - 2 * max_diameter, np.max(x) + 10 * max_diameter)
 
             if self.x2_bounds is None:
-                self.x2_bounds = (np.min(y) - 2 * max_diameter, np.max(y) + 2 * max_diameter)
+                # dh 
+                if fix_orientation :
+                   self.x2_bounds = (np.min(y) - 10 * max_diameter, np.max(y) + 10 * max_diameter)
+                else : 
+                   self.x2_bounds = (np.min(y) - 2 * max_diameter, np.max(y) + 2 * max_diameter)
 
             # TODO figure out proper z spacing for GCH, currently set to +/- 10.0
             x_points, y_points, z_points = np.meshgrid(
@@ -377,6 +389,9 @@ def set_grid(self) -> None:
                 ]),
                 indexing="ij"
             )
+            
+            if fix_orientation : 
+                x_points, y_points, z_points = rotate_coordinates_rel_west(self.wind_directions, (x_points, y_points, z_points), inv_rot=False )
 
             self.x_sorted = x_points[None, None, :, :, :]
             self.y_sorted = y_points[None, None, :, :, :]