diff --git a/ProcessLib/ProcessVariable.cpp b/ProcessLib/ProcessVariable.cpp
index 2f4bfaf0d0a6c1b850ebc4ba6708de52e45733aa..3f9d64a5eb8d8bb1a68cc511a0458af5983c2423 100644
--- a/ProcessLib/ProcessVariable.cpp
+++ b/ProcessLib/ProcessVariable.cpp
@@ -19,7 +19,7 @@
 namespace ProcessLib
 {
 ProcessVariable::ProcessVariable(BaseLib::ConfigTree const& config,
-                                 MeshLib::Mesh const& mesh,
+                                 MeshLib::Mesh& mesh,
                                  GeoLib::GEOObjects const& geometries)
     : _name(config.getConfParam<std::string>("name")),
       _mesh(mesh),
diff --git a/ProcessLib/ProcessVariable.h b/ProcessLib/ProcessVariable.h
index 8268035611cfef9a67c2f811c16e2ad9a65e66bd..3efe3671772638a6fbdef7e42cf63495cefd43c1 100644
--- a/ProcessLib/ProcessVariable.h
+++ b/ProcessLib/ProcessVariable.h
@@ -47,7 +47,7 @@ namespace ProcessLib
 class ProcessVariable
 {
 public:
-	ProcessVariable(BaseLib::ConfigTree const& config, MeshLib::Mesh const& mesh,
+	ProcessVariable(BaseLib::ConfigTree const& config, MeshLib::Mesh& mesh,
 	                GeoLib::GEOObjects const& geometries);
 
 	ProcessVariable(ProcessVariable&&);
@@ -98,7 +98,7 @@ public:
 
 private:
 	std::string const _name;
-	MeshLib::Mesh const& _mesh;
+	MeshLib::Mesh& _mesh;
 	int _tuple_size;
 	std::unique_ptr<InitialCondition> _initial_condition;
 	std::vector<std::unique_ptr<UniformDirichletBoundaryCondition>>