diff --git a/scri/asymptotic_bondi_data/from_initial_values.py b/scri/asymptotic_bondi_data/from_initial_values.py
index 2188e37c..1d1041b4 100644
--- a/scri/asymptotic_bondi_data/from_initial_values.py
+++ b/scri/asymptotic_bondi_data/from_initial_values.py
@@ -103,7 +103,7 @@ def asany_atleast2d_complex(a):
         # ψ₂ = ∫ (-ðψ₃ + σψ₄) du
         ψ2_0 = (
             ModesTimeSeries(psi2, u, spin_weight=0, multiplication_truncator=max).real
-            - 1j * (σ_0.bar.eth.eth + σ_0 * σ_1.bar).imag
+            - 1j * (σ_0.bar.eth_GHP.eth_GHP + σ_0 * σ_1.bar).imag
         )
         ψ2_1 = σ_0 * ψ4_0 + ð(ψ3_0)
         ψ2_2 = (σ_1 * ψ4_0 + ð(ψ3_1)) / 2
@@ -134,7 +134,7 @@ def adjust_psi2_imaginary_part(psi2, abd):
             sigma_bar_dot_initial = abd.sigma.bar.dot[..., 0, :]
             return (
                 ModesTimeSeries(psi2, abd.time, spin_weight=0).real
-                - 1j * (sigma_initial.bar.eth.eth + sigma_initial * sigma_bar_dot_initial).imag
+                - 1j * (sigma_initial.bar.eth_GHP.eth_GHP + sigma_initial * sigma_bar_dot_initial).imag
             )
 
         abd.sigma = sigma0