diff --git a/sbmy_control/parameters_card.py b/sbmy_control/parameters_card.py
index 5384adb..cfeda44 100644
--- a/sbmy_control/parameters_card.py
+++ b/sbmy_control/parameters_card.py
@@ -452,6 +452,35 @@ def create_parameter_card_dict(
             MeshDensity=MeshDensity,
             LoadPhiBCs=int(LoadPhiBCs),
             InputPhiBCsBase=InputPhiBCsBase,
+            
+            
+            # P3M default parameters
+            n_Tiles                      = 60,
+            # Number of tiles in each dimension for the P3M code
+            nPairsForceDiagnostic               =0,
+            # Number of particle-hole pairs to be used for the force diagnostic
+            nBinsForceDiagnostic               =0,
+            maxTrialsForceDiagnostic             =0,
+            OutputForceDiagnostic ='',
+            PrintOutputTimestepsLog             = 0,
+            OutputTimestepsLog               = "",
+            fac_dyn_custom                 = 0.03,
+            da_max_early_custom               = 0.0013,
+            fac_H_custom                 =  0.05,
+            fac_bend                    = 0.3125,
+            sub_bend1                    = 0.012,
+            sub_bend2                    = 0.014,
+            fac_p3m_fit                   = 0.15,
+            da_max_late_custom               = 0.02,
+            # Polynomial coefficients for the P3M time step limiter. Only used if use_p3m_fit=1
+            use_p3m_fit                   = 1,
+            p3m_fit_coeff0                 = -3.19544528,
+            p3m_fit_coeff1                 = -0.01948223,
+            p3m_fit_coeff2                 = 0.5180241,
+            p3m_fit_coeff3                 = -1.08372045,
+            p3m_fit_coeff4                 = -0.71808255,
+            p3m_fit_coeff5                 = -0.14370312,
+            RunForceDiagnostic               = 0,
 
             h=h,
             Omega_m=Omega_m,