diff --git a/fem/src/MeshUtils.F90 b/fem/src/MeshUtils.F90
index f2f80ef40d..6f28894317 100644
--- a/fem/src/MeshUtils.F90
+++ b/fem/src/MeshUtils.F90
@@ -25552,6 +25552,7 @@ SUBROUTINE CalculateMeshPieces( Mesh, ElementMode, PieceIndex)
     CALL Info('CalculateMeshPieces','Mesh coloring loops: '//I2S(Loop),Level=6)
 
     ! Compute the true number of different pieces
+    MaxIndex = MAXVAL( MeshPiece )
     IF( MaxIndex == 1 ) THEN
       NoPieces = 1
       IF(PRESENT(PieceIndex)) PieceIndex = 1