diff --git a/simulations/airshower_beacon_simulation/bc_baseline_phase_deltas.py b/simulations/airshower_beacon_simulation/bc_baseline_phase_deltas.py
index ce1783b..e8b5fec 100755
--- a/simulations/airshower_beacon_simulation/bc_baseline_phase_deltas.py
+++ b/simulations/airshower_beacon_simulation/bc_baseline_phase_deltas.py
@@ -81,7 +81,7 @@ if __name__ == "__main__":
     ##############################
     # Compare actual time shifts #
     ##############################
-    actual_antenna_true_phases = { a.name: -2*np.pi*a.attrs['clock_offset']*f_beacon for a in antennas }
+    actual_antenna_true_phases = { a.name: -2*np.pi*a.attrs['clock_offset']*f_beacon for a in sorted(antennas, key=lambda a: int(a.name)) }
 
     # Compare actual time shifts
     my_phase_diffs = []
@@ -112,10 +112,10 @@ if __name__ == "__main__":
                 phase_residuals = lib.phase_mod(phase_diffs[:,1] - my_phase_diffs)
 
                 fig.suptitle("Difference between Measured and Actual phase difference\n for Baselines (i,j" + (')' if ref_ant_id is None else '='+str([ int(a.name) for a in ref_ants])+')'))
-                axs[-1].set_xlabel("Baseline Phase Residual $\\varphi_{ij_{meas}} - \\varphi_{ij_{true}}$ [rad]")
+                axs[-1].set_xlabel("Baseline Phase Residual $\\Delta\\varphi_{ij_{meas}} - \\Delta\\varphi_{ij_{true}}$ [rad]")
             else:
                 fig.suptitle("Comparison Measured and Actual phase difference\n for Baselines (i,j"  + (')' if ref_ant_id is None else '='+str([ int(a.name) for a in ref_ants])+')'))
-                axs[-1].set_xlabel("Baseline Phase $\\varphi_{ij}$ [rad]")
+                axs[-1].set_xlabel("Baseline Phase $\\Delta\\varphi_{ij}$ [rad]")
 
 
             i=0