let the user pick low and high res parameters, but let the code use defaults if user does not provide them

Author: ptrbortolotti

openfast_toolbox/fastfarm/FASTFarmCaseCreation.py

@@ -325,7 +325,7 @@ def _checkInputs(self):
         self.fmax      = self.wts[0]['fmax']
         self.Cmeander  = self.wts[0]['Cmeander']
 
-        if self.inflowType == 'TS':
+        if self.inflowType == 'TS' and self.dt_high == None:
             self.dt_high = 1./(2.*self.fmax)
 
         # Check the platform heading and initialize as zero if needed
@@ -507,19 +507,21 @@ def _checkInputs(self):
 
         # Check the ds and dt for the high- and low-res boxes. If not given, call the
         # AMR-Wind auxiliary function with dummy domain limits. 
-        if None in (self.dt_high, self.ds_high, self.dt_low, self.ds_low):
+        if None in (self.dt_high, self.ds_high, self.dt_low, self.ds_low) and self.inflowType == 'LES':
             mod_wake_str = ['','polar', 'curled', 'cartesian']
             print(f'WARNING: One or more temporal or spatial resolution for low- and high-res domains were not given.')
             print(f'         Estimated values for {mod_wake_str[self.mod_wake]} wake model shown below.')
             self._determine_resolutions_from_dummy_amrwind_grid()
 
-        # Check the domain extents
-        if self.dt_low%(self.dt_high-1e-15) > 1e-12:
-            raise ValueError(f'The temporal resolution dT_Low should be a multiple of dT_High')
-        if self.dt_low < self.dt_high:
-            raise ValueError(f'The temporal resolution dT_High should not be greater than dT_Low on the LES side')
-        if self.ds_low < self.ds_high:
-            raise ValueError(f'The grid resolution dS_High should not be greater than dS_Low on the LES side')
+        # Check the domain extents when values are provided by the user:
+        if self.dt_low != None and self.dt_high!= None:
+            if self.dt_low%(self.dt_high-1e-15) > 1e-12:
+                raise ValueError(f'The temporal resolution dT_Low should be a multiple of dT_High')
+            if self.dt_low < self.dt_high:
+                raise ValueError(f'The temporal resolution dT_High should not be greater than dT_Low on the LES side')
+        if self.ds_low != None and self.ds_high!= None:
+            if self.ds_low < self.ds_high:
+                raise ValueError(f'The grid resolution dS_High should not be greater than dS_Low on the LES side')
 
 
         # Check the reference turbine for rotation
@@ -1899,7 +1901,7 @@ def TS_high_get_time_series(self):
                         wvel = np.roll(bts['u'][2, :, jTurb, kTurb], start_time_step)
 
                         # Map it to high-res time and dt (both) 
-                        time_hr = np.arange(0, self.tmax + self.dt_high, self.dt_high)
+                        time_hr = np.arange(bts.t[0], self.tmax + self.dt_high, self.dt_high)
                         uvel_hr = np.interp(time_hr, time, uvel)
                         vvel_hr = np.interp(time_hr, time, vvel)
                         wvel_hr = np.interp(time_hr, time, wvel)

openfast_toolbox/fastfarm/examples/Ex3_FFarmCompleteSetup.py

@@ -66,19 +66,6 @@ def main():
     # Set the yaw of each turbine for wind dir. One row for each wind direction.
     yaw_init = None
 
-    # ----------- Low- and high-res boxes parameters
-    # Should match LES if comparisons are to be made; otherwise, set desired values
-    # For an automatic computation of such parameters, omit them from the call to FFCaseCreation
-    # High-res boxes settings
-    dt_high_les = 0.6                # sampling frequency of high-res files
-    ds_high_les = 10.0               # dx, dy, dz that you want these high-res files at
-    extent_high = 1.2               # high-res box extent in y and x for each turbine, in D.
-    # Low-res boxes settings
-    dt_low_les  = 3                  # sampling frequency of low-res files
-    ds_low_les  = 20.0               # dx, dy, dz of low-res files
-    extent_low  = [3, 8,  3, 3, 2]   # extent in xmin, xmax, ymin, ymax, zmax, in D
-    
-    
     # ----------- Execution parameters
     ffbin = '/full/path/to/your/binary/.../bin/FAST.Farm'
 
@@ -132,8 +119,6 @@ def main():
 
     # Initial setup
     case = FFCaseCreation(path, wts, tmax, zbot, vhub, shear, TIvalue, inflow_deg,
-                          dt_high_les, ds_high_les, extent_high,
-                          dt_low_les, ds_low_les, extent_low,
                           ffbin=ffbin, mod_wake=mod_wake, yaw_init=yaw_init,
                           nSeeds=nSeeds,
                           inflowType=inflowType,