Get number of points along the blade

Exawind · Aug 3, 2023 · 92b7f68 · 92b7f68
1 parent c2f4324
commit 92b7f68
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/amr-wind/wind_energy/actuator/turbine/fast/turbine_fast_ops.H b/amr-wind/wind_energy/actuator/turbine/fast/turbine_fast_ops.H
@@ -146,7 +146,7 @@ struct InitDataOp<TurbineFast, SrcTrait>
         amrex::ParallelDescriptor::Comm_dup(
             amrex::ParallelDescriptor::Communicator(), tdata.tcomm);
 
-        int sz_info[4]{0, 0, 0, 0};
+        int sz_info[5]{0, 0, 0, 0, 0};
         if (info.is_root_proc) {
             tdata.fast->init_turbine(tdata.fast_data.tid_local);
 
@@ -155,6 +155,7 @@ struct InitDataOp<TurbineFast, SrcTrait>
             sz_info[1] = tf.to_cfd.fx_Len;
             sz_info[2] = tf.from_cfd.u_Len;
             sz_info[3] = tf.num_pts_tower;
+            sz_info[4] = tf.num_blade_elem;
         }
 
         // Broadcast data to everyone
@@ -165,7 +166,7 @@ struct InitDataOp<TurbineFast, SrcTrait>
             tdata.num_blades = sz_info[0];
             // back calculate what the value per blade for number of points in
             // the openfast data structure
-            tdata.num_vel_pts_blade = sz_info[2] / tdata.num_blades;
+            tdata.num_vel_pts_blade = sz_info[4];
             data.grid().resize(sz_info[1], sz_info[2]);
             tdata.chord.resize(sz_info[1]);
             tdata.num_pts_tower = sz_info[3];