diff --git a/dumux/python/assembly/fvassembler.hh b/dumux/python/assembly/fvassembler.hh
index 1ede28dc7574a7a7f20bfa754c90a5489fa86afe..23642efb0e1da4b03ca85356fe7ae29339d020cc 100644
--- a/dumux/python/assembly/fvassembler.hh
+++ b/dumux/python/assembly/fvassembler.hh
@@ -59,10 +59,6 @@ void registerFVAssembler(pybind11::handle scope, pybind11::class_<FVAssembler, o
     cls.def_property_readonly("isStationaryProblem", &FVAssembler::isStationaryProblem);
     cls.def_property_readonly("gridVariables", [](FVAssembler& self) { return self.gridVariables(); });
 
-    cls.def("updateGridVariables", [](FVAssembler& self, const SolutionVector& curSol){
-        self.updateGridVariables(curSol);
-    });
-
     cls.def("assembleResidual", [](FVAssembler& self, const SolutionVector& curSol){
         self.assembleResidual(curSol);
     });
diff --git a/test/python/test_1p.py b/test/python/test_1p.py
index 80ff3088c16cd59e2ac085371dcc5716d01dd257..029be21f61d2271d1b8ce83d883fa6b955c294aa 100755
--- a/test/python/test_1p.py
+++ b/test/python/test_1p.py
@@ -133,7 +133,6 @@ gridVars = GridVariables(problem=problem, model=model)
 assembler = FVAssembler(problem=problem, gridVariables=gridVars, model=model, diffMethod=diffMethod)
 sol = blockVector(assembler.numDofs)
 gridVars.init(sol)
-assembler.updateGridVariables(sol)
 print("numdofs", assembler.numDofs)
 
 # Assemble the Jacobian and the residual