Small change in zeta_dot calculation
@@ -476,8 +476,8 @@ def generate_strip(node_info, airfoil_db, aligned_grid, orientation_in=np.array(
for i_M in range(node_info['M'] + 1):
# zeta_dot_a_frame[:, i_M] += (
# np.dot(Cab, np.cross(omega_b, strip_coordinates_b_frame[:, i_M])))
- zeta_dot_a_frame[:, i_M] -= (
- np.dot(Cab, np.cross(omega_b, strip_coordinates_b_frame[:, i_M])))
+ zeta_dot_a_frame[:, i_M] += (
+ np.cross(np.dot(Cab, np.dot(omega_b, Cab.T)), strip_coordinates_a_frame[:, i_M]))
# np.cross(omega_b, strip_coordinates_b_frame[:, i_M]))
# np.cross(omega_b, strip_coordinates_a_frame[:, i_M]))
