diff --git a/Python/pywarpx/picmi.py b/Python/pywarpx/picmi.py
index 42c0313a6e9..d06fe90d2d6 100644
--- a/Python/pywarpx/picmi.py
+++ b/Python/pywarpx/picmi.py
@@ -1351,6 +1351,9 @@ def init(self, kw):
 
     def solver_initialize_inputs(self):
 
+        # Open BC means FieldBoundaryType::Open for electrostatic sims, rather than perfectly-matched layer
+        BC_map['open'] = 'open'
+
         self.grid.grid_initialize_inputs()
 
         if self.relativistic:
@@ -1371,6 +1374,8 @@ def solver_initialize_inputs(self):
             pywarpx.boundary.potential_hi_y = self.grid.potential_ymax
             pywarpx.boundary.potential_hi_z = self.grid.potential_zmax
 
+        pywarpx.warpx.poisson_solver = self.method
+
 
 class GaussianLaser(picmistandard.PICMI_GaussianLaser):
     def laser_initialize_inputs(self):