diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/Allrun b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/Allrun
index 527adc0a3eb6f79b039855562ac4fd09c89c8a04..0e3d605ef8bd4a3e95c008268fd86b0ea58e0025 100755
--- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/Allrun
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/Allrun
@@ -16,8 +16,14 @@ runApplication blockMesh
 runApplication surfaceFeatureExtract
 
 runApplication snappyHexMesh -overwrite
-runApplication topoSet -dict system/topoSetDict.f1
 
+#- Alternatively if snappyHexMesh has not set the faceZones correctly
+#  can do some topoSet magic:
+#   - pick up the outside of the cellZones
+#   - put into faceZones
+#runApplication topoSet -dict system/topoSetDict.f1
+
+# Use faceZones in combination with preserveFaceZones constraint
 runApplication $decompDict decomposePar
 
 # Restore initial fields
diff --git a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/system/decomposeParDict b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/system/decomposeParDict
index 7c20dbe0172698b1d312fb35284043e01f232c29..f980783689c1930b31b406f07a46ada3c504383e 100644
--- a/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/system/decomposeParDict
+++ b/tutorials/heatTransfer/chtMultiRegionFoam/snappyMultiRegionHeaterImplicit/system/decomposeParDict
@@ -27,7 +27,17 @@ constraints
     faces
     {
         type    preserveFaceZones;
-        zones   (f1 f2 f3 f4 f5);
+        zones
+        (
+            heater_to_bottomAir
+            heater_to_leftSolid
+            heater_to_rightSolid
+            heater_to_topAir
+            leftSolid_to_bottomAir
+            leftSolid_to_topAir
+            rightSolid_to_bottomAir
+            rightSolid_to_topAir
+        );
         enabled true;
     }
 }