From 7892f85679bf5e47d4ab8ad43ec58b039cf45882 Mon Sep 17 00:00:00 2001
From: richb2k <75949649+richb2k@users.noreply.github.com>
Date: Wed, 2 Oct 2024 16:59:16 -0400
Subject: [PATCH] set the value of MaxIndex based on the maximum value in the
 MeshPiece vector.

---
 fem/src/MeshUtils.F90 | 1 +
 1 file changed, 1 insertion(+)

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