@@ -500,16 +500,38 @@ function calculate_results(wa::WingAerodynamics,
500500 dot (ftotal_induced_va, [0.0 , 1.0 , 0.0 ]),
501501 dot (ftotal_induced_va, [0.0 , 0.0 , 1.0 ])
502502 ] * panel. width
503+
504+ # ### Converting forces to the global reference frame
505+ # fx_global_2D = jit_dot(ftotal_induced_va, np.array([1, 0, 0]))
506+ # fy_global_2D = jit_dot(ftotal_induced_va, np.array([0, 1, 0]))
507+ # fz_global_2D = jit_dot(ftotal_induced_va, np.array([0, 0, 1]))
508+
509+ # fx_global_3D = fx_global_2D * panel_width
510+ # fy_global_3D = fy_global_2D * panel_width
511+ # fz_global_3D = fz_global_2D * panel_width
503512
504513 # Update sums
505514 lift_wing_3D_sum += lift_prescribed_va * panel. width
506515 drag_wing_3D_sum += drag_prescribed_va * panel. width
507516 side_wing_3D_sum += side_prescribed_va * panel. width
517+
518+ # TODO make this work
519+ # fx_global_3D_sum += fx_global_3D
520+ # fy_global_3D_sum += fy_global_3D
521+ # fz_global_3D_sum += fz_global_3D
508522
509523 # Store coefficients
510524 push! (cl_prescribed_va, lift_prescribed_va / (q_inf * panel. chord))
511525 push! (cd_prescribed_va, drag_prescribed_va / (q_inf * panel. chord))
512526 push! (cs_prescribed_va, side_prescribed_va / (q_inf * panel. chord))
527+
528+ # TODO translate this
529+ # fx_global_3D_list.append(fx_global_3D)
530+ # fy_global_3D_list.append(fy_global_3D)
531+ # fz_global_3D_list.append(fz_global_3D)
532+ # f_global_3D_list.append(
533+ # np.array([fx_global_3D, fy_global_3D, fz_global_3D])
534+ # )
513535 end
514536
515537 if is_only_f_and_gamma_output
0 commit comments