diff --git a/IMRPDT_waveform/IMRPhenomDNRTparity.py b/IMRPDT_waveform/IMRPhenomDNRTparity.py
index 5878ee445663e24227c05b109dad1eae3582e65a..4e7e203be061605a10cae12d2f4be90f93ea3efc 100644
--- a/IMRPDT_waveform/IMRPhenomDNRTparity.py
+++ b/IMRPDT_waveform/IMRPhenomDNRTparity.py
@@ -6,7 +6,7 @@ def gen(**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
+    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
diff --git a/waveform/IMRPhenomParity.py b/waveform/IMRPhenomParity.py
index 5f962cfd9f3d78704d765fd3241e216c3e68c53a..fca51a2eca3cecb578bb8b274697bcb7865310ef 100644
--- a/waveform/IMRPhenomParity.py
+++ b/waveform/IMRPhenomParity.py
@@ -1,12 +1,13 @@
 def genIMRPhenomParity(**kwds):
     from pycbc.waveform import get_fd_waveform
     import lal
+    #print(kwds)
 
     if 'approximant' in kwds:
         kwds.pop("approximant")
     hp, hc = get_fd_waveform(approximant="IMRPhenomXPHM", **kwds)
 
-    temp =  kwds['parity_Amu'] / 1e9 / lal.QE_SI * lal.H_SI * lal.PI * lal.PI / lal.H0_SI
+    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