diff --git a/src/language/modules/SchemeModule.cpp b/src/language/modules/SchemeModule.cpp
index 1c888821f8d6c9c62e8ddc58ce23cd97c4f05b32..9fc6d68bb0ae5295925cb6bcefce8308cbbf5cf2 100644
--- a/src/language/modules/SchemeModule.cpp
+++ b/src/language/modules/SchemeModule.cpp
@@ -30,6 +30,7 @@
 #include <scheme/DiscreteFunctionVectorInterpoler.hpp>
 #include <scheme/ExternalBoundaryConditionDescriptor.hpp>
 #include <scheme/FixedBoundaryConditionDescriptor.hpp>
+#include <scheme/FluxingAdvectionSolver.hpp>
 #include <scheme/FourierBoundaryConditionDescriptor.hpp>
 #include <scheme/FreeBoundaryConditionDescriptor.hpp>
 #include <scheme/HyperelasticSolver.hpp>
@@ -589,6 +590,13 @@ SchemeModule::SchemeModule()
                                                  }
 
                                                  ));
+  this->_addBuiltinFunction("fluxing_advection", std::function(
+
+                                                   [](const std::shared_ptr<const IMesh> new_mesh,
+                                                      const std::shared_ptr<const DiscreteFunctionVariant>& old_q)
+                                                     -> std::shared_ptr<const DiscreteFunctionVariant> {
+                                                     return FluxingAdvectionSolverHandler(new_mesh, old_q);
+                                                   }));
 
   MathFunctionRegisterForVh{*this};
 }
diff --git a/src/scheme/FluxingAdvectionSolver.cpp b/src/scheme/FluxingAdvectionSolver.cpp
index cb490acd59eefb6faee56f25a3d15f27152460dc..d3fcd4bb165e4c2b02c88157f97234e605eaf4a6 100644
--- a/src/scheme/FluxingAdvectionSolver.cpp
+++ b/src/scheme/FluxingAdvectionSolver.cpp
@@ -85,7 +85,7 @@ class FluxingAdvectionSolver
   FluxingAdvectionSolver(const std::shared_ptr<const MeshType> old_mesh,
                          const std::shared_ptr<const MeshType> new_mesh,
                          const DiscreteFunctionP0<Dimension, const double>& old_q)
-    : m_old_mesh{old_mesh}, m_new_mesh{new_mesh}, m_old_q{old_q}
+    : m_old_mesh{old_mesh}, m_new_mesh{new_mesh}
   {}
 
   ~FluxingAdvectionSolver() = default;