diff --git a/ProcessLib/Process.h b/ProcessLib/Process.h
index 753b27d5db406c736378f6991cfba77226a72269..d3e767e3330e54970d6553ec462b25128f38f99e 100644
--- a/ProcessLib/Process.h
+++ b/ProcessLib/Process.h
@@ -10,6 +10,9 @@
 #ifndef PROCESS_LIB_PROCESS_H_
 #define PROCESS_LIB_PROCESS_H_
 
+#include "AssemblerLib/MeshComponentMap.h"
+#include "MeshLib/Elements/Element.h"
+#include "MeshLib/Node.h"
 #include "MeshLib/Mesh.h"
 
 namespace ProcessLib
@@ -26,6 +29,40 @@ public:
 
     virtual void initialize() = 0;
 
+    std::vector<std::vector<std::size_t>>
+    createDofMap(AssemblerLib::MeshComponentMap& mesh_component_map) const
+    {
+        // dof_map contains for each element vector of global indices to
+        // node/element process variables.
+        std::vector<std::vector<std::size_t>> dof_map;
+
+        // XXX Shouldn't the mesh_id (and the elements) be from the MeshSubsets?
+        // Otherwise we should create the MeshComponentMap also here using the
+        // same meshes and elements.
+        std::size_t const mesh_id = _mesh.getID();
+        std::vector<MeshLib::Element*> const& elements = _mesh.getElements();
+
+        dof_map.reserve(elements.size());
+
+        // For each element find the global indices for node/element components.
+        for (MeshLib::Element* e : elements)
+        {
+            std::size_t const nnodes = e->getNNodes();
+            std::vector<MeshLib::Location> vec_items;
+            vec_items.reserve(nnodes);
+            for (std::size_t j = 0; j < nnodes; j++)
+                vec_items.emplace_back(
+                    mesh_id,
+                    MeshLib::MeshItemType::Node,
+                    e->getNode(j)->getID());
+
+            dof_map.push_back(mesh_component_map.getGlobalIndices<
+                AssemblerLib::ComponentOrder::BY_COMPONENT>(vec_items));
+        }
+
+        return dof_map;
+    }
+
 protected:
     MeshLib::Mesh const& _mesh;
     std::size_t const _numberOfNodeComponents;