diff --git a/dune/tectonic/data-structures/body/body.cc b/dune/tectonic/data-structures/body/body.cc
index 83fb4fa86971ec3e903f92e61bb88c7f4cde4abe..191cf4d60f213918db2e7e10be669cb90e73f814 100644
--- a/dune/tectonic/data-structures/body/body.cc
+++ b/dune/tectonic/data-structures/body/body.cc
@@ -94,6 +94,12 @@ auto LeafBody<GridTEMPLATE, VectorType>::grid() const
     return grid_;
 }
 
+template <class GridTEMPLATE, class VectorType>
+auto LeafBody<GridTEMPLATE, VectorType>::hostGrid() const
+-> std::shared_ptr<GridTEMPLATE> {
+
+    return hostGrid_;
+}
 
 template <class GridTEMPLATE, class VectorType>
 auto LeafBody<GridTEMPLATE, VectorType>::gridView() const
diff --git a/dune/tectonic/data-structures/body/body.hh b/dune/tectonic/data-structures/body/body.hh
index 74809e5048f1b488baf4226da9d9473fb9f1721a..7f6984351f65b30be876d678848b83b3534f11dc 100644
--- a/dune/tectonic/data-structures/body/body.hh
+++ b/dune/tectonic/data-structures/body/body.hh
@@ -90,6 +90,7 @@ class LeafBody {
     void setDeformation(const VectorType& def);
     auto deformation() const -> DeformationFunction&;
     auto grid() const -> std::shared_ptr<GridType>;
+    auto hostGrid() const -> std::shared_ptr<HostGridType>;
 
     auto gridView() const -> const GridView&;
     auto nVertices() const -> size_t;