diff --git a/IMRPDT_waveform/.DS_Store b/IMRPDT_waveform/.DS_Store
new file mode 100644
index 0000000000000000000000000000000000000000..5008ddfcf53c02e82d7eee2e57c38e5672ef89f6
Binary files /dev/null and b/IMRPDT_waveform/.DS_Store differ
diff --git a/IMRPDT_waveform/IMRPhenomDNRTparity.py b/IMRPDT_waveform/IMRPhenomDNRTparity.py
new file mode 100644
index 0000000000000000000000000000000000000000..5878ee445663e24227c05b109dad1eae3582e65a
--- /dev/null
+++ b/IMRPDT_waveform/IMRPhenomDNRTparity.py
@@ -0,0 +1,22 @@
+def gen(**kwds):
+    from pycbc.waveform import get_fd_waveform
+    import lal
+
+    if 'approximant' in kwds:
+        kwds.pop("approximant")
+    hp, hc = get_fd_waveform(approximant="IMRPhenomD_NRTidal", **kwds)
+
+    temp =  kwds['parity_Amu'] / 1e9 / lal.QE_SI * lal.H_SI * lal.PI * lal.PI / lal.H0_SI
+
+    hp_parity = hp + hc * temp * hp.sample_frequencies **2
+    hc_parity = hc - hp * temp * hp.sample_frequencies **2
+
+    return hp_parity, hc_parity
+
+
+#pycbc.waveform.add_custom_waveform('IMRPhenomParity', IMRPhenomParity, 'frequency', force=True)
+
+def add_me(**kwds):
+    kwds['cpu_fd']['IMRPhenomD_NRTidal_Parity'] = gen
+    kwds['filter_time_lengths']['IMRPhenomD_NRTidal_Parity'] = kwds['filter_time_lengths']['IMRPhenomD']
+