format

docs/src/water_activity_plots/T_vs_wateractivity.jl

@@ -27,9 +27,8 @@ p_sat_ice = [TD.saturation_vapor_pressure(thermo_params, T, TD.Ice()) for T in T
 a_w = p_sol_1 ./ p_sat_liq                 # a_w current parameterization
 a_w_alternate = p_sol_1 ./ p_sol_0         # a_w if sat vapor pressure over pure liq water was calculated with p_sol eqn
 a_w_ice = p_sat_ice ./ p_sat_liq           # a_ice current parameterization
-a_w_ice_alternate = p_sat_ice ./ p_sol_0          # a_ice if sat vapor pressure over pure liq water was calculated with p_sol eqn
+a_w_ice_alternate = p_sat_ice ./ p_sol_0   # a_ice if sat vapor pressure over pure liq water was calculated with p_sol eqn
 a_w_ice_μ = [exp((210368 + 131.438*T - (3.32373e6 /T) - 41729.1*log(T))/(8.31441*T)) for T in T_range]  # a_ice using chemical potential parameterization
-#! format: on
 
 # Plotting. NOTE: all ".* -1" is only to flip x axis
 fig = MK.Figure(resolution = (800, 600))
@@ -42,3 +41,4 @@ MK.lines!(ax1, a_w_ice_μ .* -1, T_range, label = "a_w_ice using μ", color = :l
 MK.axislegend()
 
 MK.save("T_vs_wateractivity.svg", fig)
+#! format: on