diff --git a/ProcessLib/ProcessVariable.cpp b/ProcessLib/ProcessVariable.cpp
index 28b79663ed25facabf688164d75ee0bff11e5808..1810f6569c3801f6e1f002535d49b4b61004ee04 100644
--- a/ProcessLib/ProcessVariable.cpp
+++ b/ProcessLib/ProcessVariable.cpp
@@ -176,12 +176,6 @@ MeshLib::Mesh const& ProcessVariable::getMesh() const
     return _mesh;
 }
 
-MeshLib::PropertyVector<double>& ProcessVariable::getOrCreateMeshProperty()
-{
-    return *MeshLib::getOrCreateMeshProperty<double>(
-        _mesh, _name, MeshLib::MeshItemType::Node, _n_components);
-}
-
 std::vector<std::unique_ptr<BoundaryCondition>>
 ProcessVariable::createBoundaryConditions(
     const NumLib::LocalToGlobalIndexMap& dof_table,
diff --git a/ProcessLib/ProcessVariable.h b/ProcessLib/ProcessVariable.h
index 383e9761c8fab8a3e1ff86d8be6fbc330c8e0be3..a343b165a01b34f88f663e99100c1e622860587b 100644
--- a/ProcessLib/ProcessVariable.h
+++ b/ProcessLib/ProcessVariable.h
@@ -69,11 +69,6 @@ public:
         return _initial_condition;
     }
 
-    // Get or create a property vector for results.
-    // The returned mesh property size is number of mesh nodes times number of
-    // components.
-    MeshLib::PropertyVector<double>& getOrCreateMeshProperty();
-
     unsigned getShapeFunctionOrder() const { return _shapefunction_order; }
 
 private: