fix td waveform epoch

bilby/gw/source.py

@@ -364,7 +364,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)