diff --git a/lisainstrument/instrument.py b/lisainstrument/instrument.py
index 42db45ae8e411073b2d3787ce760f3831e1698a6..1bc6cfc26c8a5e925857b1716d1e6e24226a017b 100755
--- a/lisainstrument/instrument.py
+++ b/lisainstrument/instrument.py
@@ -905,14 +905,14 @@ class Instrument:
         logger.debug("Computing local timer deviations")
         t = self.physics_et_withinitial
 
-        print(self.tps_wrt_tcb['1'].shape)
-        print(t.shape)
         self.the_wrt_tcb_withinitial = \
             self.tps_wrt_tcb \
-            + self.clock_offsets \
-            + self.clock_freqoffsets * (t + self.tps_wrt_tcb) \
-            + self.clock_freqlindrifts * (t + self.tps_wrt_tcb)**2 / 2 \
-            + self.clock_freqquaddrifts * (t + self.tps_wrt_tcb)**3 / 3 \
+            + ForEachSC(lambda sc:
+                self.clock_offsets[sc] \
+                + self.clock_freqoffsets[sc] * (t + self.tps_wrt_tcb[sc]) \
+                + self.clock_freqlindrifts[sc] / 2 * (t + self.tps_wrt_tcb[sc])**2 \
+                + self.clock_freqquaddrifts[sc] / 3 * (t + self.tps_wrt_tcb[sc])**3 \
+            ) \
             + ForEachSC(lambda sc:
                 cumulative_trapezoid(np.broadcast_to(
                     self.clock_noise_fluctuations_withinitial[sc],
@@ -920,7 +920,7 @@ class Instrument:
                 dx=self.physics_dt, initial=0)
             ) \
             + ForEachSC(lambda sc:
-                self.tps_wrt_tcb * self.clock_noise_fluctuations_withinitial[sc]
+                self.tps_wrt_tcb[sc] * self.clock_noise_fluctuations_withinitial[sc]
             )
 
         self.the_wrt_tcb_local = self.the_wrt_tcb_withinitial.transformed(