diff --git a/bilby/gw/source.py b/bilby/gw/source.py
index 69109364a149304926318ebd8b73fc1adf96425e..99c52b6b758d3bbe0696b2630e19c8423b172fb8 100644
--- a/bilby/gw/source.py
+++ b/bilby/gw/source.py
@@ -342,6 +342,19 @@ def _base_lal_cbc_fd_waveform(
         waveform_dictionary, lambda_1)
     lalsim_SimInspiralWaveformParamsInsertTidalLambda2(
         waveform_dictionary, lambda_2)
+        
+    for key, value in waveform_kwargs.items():
+        func = getattr(lalsim, "SimInspiralWaveformParamsInsert" + key, 0)
+        if func != 0:
+            func(waveform_dictionary, value)
+    
+    print("\nHOOOOOLAAAAA\n")        
+    print("Mband = %f" % lalsim.SimInspiralWaveformParamsLookupPhenomXHMThresholdMband(waveform_dictionary))
+    print("PV = %f" % lalsim.SimInspiralWaveformParamsLookupPhenomXPrecVersion(waveform_dictionary))
+        
+    if ('numerical_relativity_file' in waveform_kwargs) and waveform_kwargs['numerical_relativity_file'] is not None:
+        lalsim.SimInspiralWaveformParamsInsertNumRelData(
+            waveform_dictionary, waveform_kwargs['numerical_relativity_file'])
 
     if ('mode_array' in waveform_kwargs) and waveform_kwargs['mode_array'] is not None:
         mode_array = waveform_kwargs['mode_array']