diff --git a/libsrc/occ/occ_solid.hpp b/libsrc/occ/occ_solid.hpp
index 1ce2d50c..505a1617 100644
--- a/libsrc/occ/occ_solid.hpp
+++ b/libsrc/occ/occ_solid.hpp
@@ -16,6 +16,7 @@ namespace netgen
         OCCSolid(TopoDS_Shape dshape)
             : solid(TopoDS::Solid(dshape))
         { }
+        TopoDS_Solid& GetShape() { return solid; }
     };
 }
 
diff --git a/libsrc/occ/occgeom.cpp b/libsrc/occ/occgeom.cpp
index 9d71fa77..821cbf8b 100644
--- a/libsrc/occ/occgeom.cpp
+++ b/libsrc/occ/occgeom.cpp
@@ -1261,6 +1261,24 @@ namespace netgen
           }
       }
 
+      // Propagate maxh to children
+      for(auto& solid : solids)
+      {
+        auto& shape = static_cast<OCCSolid&>(*solid).GetShape();
+        if(!OCCGeometry::HaveProperties(shape))
+          continue;
+        for(auto& f : GetFaces(shape))
+        {
+          auto& face = GetFace(f);
+          face.properties.maxh = min2(face.properties.maxh,
+                                      GetProperties(shape).maxh);
+        }
+      }
+      for(auto& face : faces)
+        for(auto& edge : face->edges)
+          edge->properties.maxh = min2(edge->properties.maxh,
+                                       face->properties.maxh);
+
       // Add identifications
       auto add_identifications = [&](auto & shapes, auto & shape_map)
       {
diff --git a/libsrc/occ/python_occ_shapes.cpp b/libsrc/occ/python_occ_shapes.cpp
index fbeb8e5a..f2992d73 100644
--- a/libsrc/occ/python_occ_shapes.cpp
+++ b/libsrc/occ/python_occ_shapes.cpp
@@ -836,12 +836,7 @@ DLL_HEADER void ExportNgOCCShapes(py::module &m)
                   },
                   [](TopoDS_Shape& self, double val)
                   {
-                    for (auto typ : { TopAbs_SOLID, TopAbs_FACE,  TopAbs_EDGE, TopAbs_VERTEX })
-                      for (TopExp_Explorer e(self, typ); e.More(); e.Next())
-                      {
-                        auto & maxh = OCCGeometry::GetProperties(e.Current()).maxh;
-                        maxh = min2(val, maxh);
-                      }
+                    OCCGeometry::GetProperties(self).maxh = val;
                   }, "maximal mesh-size for shape")
     
     .def_property("hpref",