diff --git a/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py b/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py
index 6eed16305..ed6e1b834 100644
--- a/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py
+++ b/CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py
@@ -116,7 +116,7 @@ def determine_heat_flux(V_g, u, k, flux):
 # Adapter definition and initialization
 precice = Adapter(adapter_config_filename="precice-adapter-config.json")
 
-precice_dt = precice.initialize(coupling_boundary, mesh, V, write_function=f_N_function)
+precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=f_N_function)
 
 # Create a FEniCS Expression to define and control the coupling boundary values
 coupling_expression = precice.create_coupling_expression()
diff --git a/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py b/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py
index 6157e9bd6..09893e4a5 100644
--- a/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py
+++ b/FSI/cylinderFlap/OpenFOAM-FEniCS/Solid/cyl-flap.py
@@ -88,7 +88,9 @@ def remaining_boundary(x, on_boundary):
 force_boundary = AutoSubDomain(remaining_boundary)
 
 # Initialize the coupling interface
-precice_dt = precice.initialize(coupling_boundary, mesh, V, dim, fixed_boundary=clamped_boundary_domain)
+# Function space V is passed twice as both read and write functions are defined using the same space
+precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=V,
+                                fixed_boundary=clamped_boundary_domain)
 
 fenics_dt = precice_dt  # if fenics_dt == precice_dt, no subcycling is applied
 # fenics_dt = 0.02  # if fenics_dt < precice_dt, subcycling is applied
diff --git a/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py b/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py
index 008ffb082..41dfeff4b 100644
--- a/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py
+++ b/FSI/flap_perp/OpenFOAM-FEniCS/Solid/perp-flap.py
@@ -72,7 +72,9 @@ def neumann_boundary(x, on_boundary):
 force_boundary = AutoSubDomain(neumann_boundary)
 
 # Initialize the coupling interface
-precice_dt = precice.initialize(coupling_boundary, mesh, V, dim, fixed_boundary=clamped_boundary_domain)
+# Function space V is passed twice as both read and write functions are defined using the same space
+precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=V,
+                                fixed_boundary=clamped_boundary_domain)
 
 fenics_dt = precice_dt  # if fenics_dt == precice_dt, no subcycling is applied
 # fenics_dt = 0.02  # if fenics_dt < precice_dt, subcycling is applied
diff --git a/HT/partitioned-heat/fenics-fenics/heat.py b/HT/partitioned-heat/fenics-fenics/heat.py
index 21c3742b0..8b98af863 100644
--- a/HT/partitioned-heat/fenics-fenics/heat.py
+++ b/HT/partitioned-heat/fenics-fenics/heat.py
@@ -124,10 +124,10 @@ def determine_gradient(V_g, u, flux):
 # Initialize the adapter according to the specific participant
 if problem is ProblemType.DIRICHLET:
     precice = Adapter(adapter_config_filename="precice-adapter-config-D.json")
-    precice_dt = precice.initialize(coupling_boundary, mesh, V, write_function=f_N_function)
+    precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=f_N_function)
 elif problem is ProblemType.NEUMANN:
     precice = Adapter(adapter_config_filename="precice-adapter-config-N.json")
-    precice_dt = precice.initialize(coupling_boundary, mesh, V_g, write_function=u_D_function)
+    precice_dt = precice.initialize(coupling_boundary, read_function_space=V_g, write_object=u_D_function)
 
 boundary_marker = False