diff --git a/floris/core/solver.py b/floris/core/solver.py
index 47e096271..a785dbee9 100644
--- a/floris/core/solver.py
+++ b/floris/core/solver.py
@@ -529,7 +529,7 @@ def cc_solver(
             + (flow_field.u_initial_sorted - turb_u_wake) * mask2
         )
 
-        turb_avg_vels = average_velocity(turb_inflow_field)
+        turb_avg_vels = average_velocity(turb_inflow_field)[:, :, None, None]
         turb_Cts = thrust_coefficient(
             turb_avg_vels,
             flow_field.turbulence_intensity_field_sorted,