fix td waveform epoch

bilby/gw/source.py

@@ -396,6 +396,13 @@ def _base_lal_cbc_fd_waveform(
     h_plus *= frequency_bounds
     h_cross *= frequency_bounds
 
+    if wf_func == lalsim_SimInspiralFD:
+        dt = 1. / delta_frequency - (hplus.epoch.gpsSeconds + hplus.epoch.gpsNanoSeconds * 1e-9)
+        h_plus *= np.exp(
+            -1j * 2 * np.pi * dt * frequency_array)
+        h_cross *= np.exp(
+            -1j * 2 * np.pi * dt * frequency_array)
+
     return dict(plus=h_plus, cross=h_cross)