Fully Integrate into Ghost

posts/2021-04-01-air-propulsion-simulation/index.qmd

@@ -124,13 +124,22 @@ Below in figure 1, the result of the simulation is plotted. Notice the massive e
 
 thrust_values = df.Thrust .|> ustrip .|> value;
 thrust_uncertainties = df.Thrust .|> ustrip .|> uncertainty;
+time_values = df.Time .|> ustrip;
 
-air = DataFrame(Thrust=thrust_values, Uncertainty=thrust_uncertainties, Time=df.Time .|> u"s" .|> ustrip);
+# Downsample for plotting (keep ~1000 points for visual fidelity)
+n_points = length(thrust_values)
+step = max(1, n_points ÷ 1000)
+idx = 1:step:n_points
 
+air = DataFrame(
+    Thrust=thrust_values[idx],
+    Uncertainty=thrust_uncertainties[idx],
+    Time=time_values[idx]
+);
 
-plot(df.Time .|> ustrip, thrust_values,
+plot(air.Time, air.Thrust,
     title="Thrust Over Time",
-    ribbon=(thrust_uncertainties, thrust_uncertainties),
+    ribbon=(air.Uncertainty, air.Uncertainty),
     fillalpha=0.2, label="Thrust",
     xlabel="Time (s)",
     ylabel="Thrust (N)",