do not use the norm of the pitch rate

ufechner7 · Nov 23, 2024 · 7e4205f · 7e4205f
1 parent 79316ec
commit 7e4205f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/KPS4.jl b/src/KPS4.jl
@@ -370,9 +370,9 @@ Updates the vector s.forces of the first parameter.
         s.forces[s.set.segments + 3] .+= (L2 + D2)
     end
     # if the impact of cmq is just drag
-    f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * norm(s.pitch_rate) * s.set.cord_length) * va_xz1
+    # f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * norm(s.pitch_rate) * s.set.cord_length) * va_xz1
     # otherwise
-    # f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * s.pitch_rate * s.set.cord_length) * va_xz1
+    f_d = -0.5 * rho * s.set.area * norm(va_xz1) * (s.set.cmq * s.pitch_rate * s.set.cord_length) * va_xz1
     s.forces[s.set.segments + 2] .+= f_d
     s.forces[s.set.segments + 4] .+= (L3 + D3 -0.5*f_d)
     s.forces[s.set.segments + 5] .+= (L4 + D4 -0.5*f_d)