Merge branch 'f/ffarmscript' of github.com:rthedin/OpenFAST-python-toolbox into f/ffarmscript

Author: rthedin

pyFAST/fastfarm/AMRWindSimulation.py

@@ -36,15 +36,18 @@ def __init__(self, wts:dict,
         self.n_cell             = n_cell
         self.max_level          = max_level
         self.incflo_velocity_hh = incflo_velocity_hh
-        self.postproc_name      = postproc_name
+        self.postproc_name_lr   = f"{postproc_name}_lr"
+        self.postproc_name_hr   = f"{postproc_name}_hr"
         self.buffer_lr          = buffer_lr
         self.buffer_hr          = buffer_hr
         self.ds_hr              = ds_hr
         self.ds_lr              = ds_lr
 
         # Placeholder variables, to be calculated by FFCaseCreation
-        self.output_frequency = None
-        self.sampling_labels = None
+        self.output_frequency_lr = None
+        self.output_frequency_hr = None
+        self.sampling_labels_lr = None
+        self.sampling_labels_hr = None
         self.nx_lr = None
         self.ny_lr = None
         self.nz_lr = None
@@ -105,14 +108,18 @@ def _calc_sampling_params(self):
         '''
 
         ### ~~~~~~~~~ Calculate high-level info for AMR-Wind sampling ~~~~~~~~~
-        sampling_labels = ["Low"]
+        sampling_labels_lr = ["Low"]
+        self.sampling_labels_lr = sampling_labels_lr
+
+        sampling_labels_hr = []
         for turbkey in self.wts:
             if 'name' in self.wts[turbkey].keys():
                 wt_name = self.wts[turbkey]['name']
             else:
                 wt_name = f'T{turbkey}'
-            sampling_labels.append(f"High{wt_name}_inflow0deg")
-        self.sampling_labels = sampling_labels
+            sampling_labels_hr.append(f"High{wt_name}_inflow0deg")
+        
+        self.sampling_labels_hr = sampling_labels_hr
 
         ### ~~~~~~~~~ Calculate timestep values and AMR-Wind plane sampling frequency ~~~~~~~~~
         ## Low resolution domain, dt_low_les
@@ -133,7 +140,8 @@ def _calc_sampling_params(self):
         self.dt_high_les = self.dt * np.floor(dt_hr_max/self.dt)  # Ensure that dt_hr is a multiple of the AMR-Wind timestep
 
         ## Sampling frequency
-        self.output_frequency = int(self.dt_high_les/self.dt)
+        self.output_frequency_lr = int(np.floor(self.dt_low_les/self.dt))
+        self.output_frequency_hr = int(np.floor(self.dt_high_les/self.dt))
 
         ### ~~~~~~~~~ Calculate grid resolutions ~~~~~~~~~
         ## Low resolution domain, ds_lr (s = x/y/z)
@@ -373,25 +381,31 @@ def write_sampling_params(self, outdir=None):
         print(f"Writing to {outfile} ...")
         with open(outfile,"w") as out:
             # Write high-level info for sampling
-            sampling_labels_str = " ".join(str(item) for item in self.sampling_labels)
+            sampling_labels_lr_str = " ".join(str(item) for item in self.sampling_labels_lr)
+            sampling_labels_hr_str = " ".join(str(item) for item in self.sampling_labels_hr)
             out.write(f"# Sampling info generated by AMRWindSamplingCreation.py\n")
-            out.write(f"incflo.post_processing    = {self.postproc_name} # averaging\n")
-            out.write(f"{self.postproc_name}.output_format    = netcdf\n")
-            out.write(f"{self.postproc_name}.output_frequency = {self.output_frequency}\n")
-            out.write(f"{self.postproc_name}.fields           = velocity # temperature tke\n")
-            out.write(f"{self.postproc_name}.labels           = {sampling_labels_str}\n")
+            out.write(f"incflo.post_processing                = {self.postproc_name_lr} {self.postproc_name_hr} # averaging\n\n")
+            out.write(f"{self.postproc_name_lr}.output_format    = netcdf\n")
+            out.write(f"{self.postproc_name_lr}.output_frequency = {self.output_frequency_lr}\n")
+            out.write(f"{self.postproc_name_lr}.fields           = velocity # temperature tke\n")
+            out.write(f"{self.postproc_name_lr}.labels           = {sampling_labels_lr_str}\n\n")
+
+            out.write(f"{self.postproc_name_hr}.output_format    = netcdf\n")
+            out.write(f"{self.postproc_name_hr}.output_frequency = {self.output_frequency_hr}\n")
+            out.write(f"{self.postproc_name_hr}.fields           = velocity # temperature tke\n")
+            out.write(f"{self.postproc_name_hr}.labels           = {sampling_labels_hr_str}\n")
 
             # Write out low resolution sampling plane info
             zoffsets_lr_str = " ".join(str(int(item)) for item in self.zoffsets_lr)
 
             out.write(f"\n# Low sampling grid spacing = {self.ds_lr} m\n")
-            out.write(f"{self.postproc_name}.Low.type         = PlaneSampler\n")
-            out.write(f"{self.postproc_name}.Low.num_points   = {self.nx_lr} {self.ny_lr}\n")
-            out.write(f"{self.postproc_name}.Low.origin       = {self.xlow_lr:.1f} {self.ylow_lr:.1f} {self.zlow_lr:.1f}\n")  # Round the float output
-            out.write(f"{self.postproc_name}.Low.axis1        = {self.xdist_lr:.1f} 0.0 0.0\n")  # Assume: axis1 oriented parallel to AMR-Wind x-axis
-            out.write(f"{self.postproc_name}.Low.axis2        = 0.0 {self.ydist_lr:.1f} 0.0\n")  # Assume: axis2 oriented parallel to AMR-Wind y-axis
-            out.write(f"{self.postproc_name}.Low.normal       = 0.0 0.0 1.0\n")
-            out.write(f"{self.postproc_name}.Low.offsets      = {zoffsets_lr_str}\n")
+            out.write(f"{self.postproc_name_lr}.Low.type         = PlaneSampler\n")
+            out.write(f"{self.postproc_name_lr}.Low.num_points   = {self.nx_lr} {self.ny_lr}\n")
+            out.write(f"{self.postproc_name_lr}.Low.origin       = {self.xlow_lr:.1f} {self.ylow_lr:.1f} {self.zlow_lr:.1f}\n")  # Round the float output
+            out.write(f"{self.postproc_name_lr}.Low.axis1        = {self.xdist_lr:.1f} 0.0 0.0\n")  # Assume: axis1 oriented parallel to AMR-Wind x-axis
+            out.write(f"{self.postproc_name_lr}.Low.axis2        = 0.0 {self.ydist_lr:.1f} 0.0\n")  # Assume: axis2 oriented parallel to AMR-Wind y-axis
+            out.write(f"{self.postproc_name_lr}.Low.normal       = 0.0 0.0 1.0\n")
+            out.write(f"{self.postproc_name_lr}.Low.offsets      = {zoffsets_lr_str}\n")
 
             # Write out high resolution sampling plane info
             for turbkey in self.hr_domains:
@@ -414,10 +428,10 @@ def write_sampling_params(self, outdir=None):
                 zoffsets_hr_str = " ".join(str(int(item)) for item in zoffsets_hr)
 
                 out.write(f"\n# Turbine {wt_name} at (x,y) = ({wt_x}, {wt_y}), with D = {wt_D}, grid spacing = {self.ds_hr} m\n")
-                out.write(f"{self.postproc_name}.{sampling_name}.type         = PlaneSampler\n")
-                out.write(f"{self.postproc_name}.{sampling_name}.num_points   = {nx_hr} {ny_hr}\n")
-                out.write(f"{self.postproc_name}.{sampling_name}.origin       = {xlow_hr:.1f} {ylow_hr:.1f} {zlow_hr:.1f}\n")  # Round the float output
-                out.write(f"{self.postproc_name}.{sampling_name}.axis1        = {xdist_hr:.1f} 0.0 0.0\n")  # Assume: axis1 oriented parallel to AMR-Wind x-axis
-                out.write(f"{self.postproc_name}.{sampling_name}.axis2        = 0.0 {ydist_hr:.1f} 0.0\n")  # Assume: axis2 oriented parallel to AMR-Wind y-axis
-                out.write(f"{self.postproc_name}.{sampling_name}.normal       = 0.0 0.0 1.0\n")
-                out.write(f"{self.postproc_name}.{sampling_name}.offsets      = {zoffsets_hr_str}\n")
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.type         = PlaneSampler\n")
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.num_points   = {nx_hr} {ny_hr}\n")
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.origin       = {xlow_hr:.1f} {ylow_hr:.1f} {zlow_hr:.1f}\n")  # Round the float output
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.axis1        = {xdist_hr:.1f} 0.0 0.0\n")  # Assume: axis1 oriented parallel to AMR-Wind x-axis
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.axis2        = 0.0 {ydist_hr:.1f} 0.0\n")  # Assume: axis2 oriented parallel to AMR-Wind y-axis
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.normal       = 0.0 0.0 1.0\n")
+                out.write(f"{self.postproc_name_hr}.{sampling_name}.offsets      = {zoffsets_hr_str}\n")

pyFAST/fastfarm/examples/Ex4_AMRWindSamplingSetup.py

@@ -41,50 +41,5 @@ def main():
     # Write out sampling parameters
     amr.write_sampling_params(outdir)
 
-    '''
-    Below are the generated contents of sampling_config.i:
-        # Sampling info generated by AMRWindSamplingCreation.py
-        incflo.post_processing = sampling # averaging
-        sampling.output_frequency = 1
-        sampling.fields = velocity # temperature tke
-        sampling.labels = Low HighT0_inflow0deg HighT1_inflow0deg HighT2_inflow0deg
-
-        # Low sampling grid spacing = 10.0 m
-        sampling.Low.type         = PlaneSampler
-        sampling.Low.num_points   = 78 206
-        sampling.Low.origin       = 885.0 2165.0 5.0
-        sampling.Low.axis1        = 770.0 0.0 0.0
-        sampling.Low.axis2        = 0.0 2050.0 0.0
-        sampling.Low.normal       = 0.0 0.0 1.0
-        sampling.Low.offsets      = 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 240 250 260 270 280
-
-        # Turbine T0 at (x,y) = (1280.0, 2560), with D = 126.9, grid spacing = 5.0 m
-        sampling.HighT0_inflow0deg.type         = PlaneSampler
-        sampling.HighT0_inflow0deg.num_points   = 32 32
-        sampling.HighT0_inflow0deg.origin       = 1197.5 2477.5 2.5
-        sampling.HighT0_inflow0deg.axis1        = 155.0 0.0 0.0
-        sampling.HighT0_inflow0deg.axis2        = 0.0 155.0 0.0
-        sampling.HighT0_inflow0deg.normal       = 0.0 0.0 1.0
-        sampling.HighT0_inflow0deg.offsets      = 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 205 210 215 220 225 230
-
-        # Turbine T1 at (x,y) = (1280.0, 3200), with D = 126.9, grid spacing = 5.0 m
-        sampling.HighT1_inflow0deg.type         = PlaneSampler
-        sampling.HighT1_inflow0deg.num_points   = 32 32
-        sampling.HighT1_inflow0deg.origin       = 1197.5 3117.5 2.5
-        sampling.HighT1_inflow0deg.axis1        = 155.0 0.0 0.0
-        sampling.HighT1_inflow0deg.axis2        = 0.0 155.0 0.0
-        sampling.HighT1_inflow0deg.normal       = 0.0 0.0 1.0
-        sampling.HighT1_inflow0deg.offsets      = 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 205 210 215 220 225 230
-
-        # Turbine T2 at (x,y) = (1280.0, 3840), with D = 126.9, grid spacing = 5.0 m
-        sampling.HighT2_inflow0deg.type         = PlaneSampler
-        sampling.HighT2_inflow0deg.num_points   = 32 32
-        sampling.HighT2_inflow0deg.origin       = 1197.5 3757.5 2.5
-        sampling.HighT2_inflow0deg.axis1        = 155.0 0.0 0.0
-        sampling.HighT2_inflow0deg.axis2        = 0.0 155.0 0.0
-        sampling.HighT2_inflow0deg.normal       = 0.0 0.0 1.0
-        sampling.HighT2_inflow0deg.offsets      = 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 205 210 215 220 225 230
-    '''
-
 if __name__ == '__main__':
     main()