From 9372a8164cdb924e3d9be8cc5d8852a629fb6d86 Mon Sep 17 00:00:00 2001
From: Bruce Perry <bruce.perry@nrel.gov>
Date: Thu, 19 Dec 2024 10:33:20 -0700
Subject: [PATCH] fix error for divu calc for manifold when its not being used
 with Null mech

---
 Source/PeleLMeX_K.H | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Source/PeleLMeX_K.H b/Source/PeleLMeX_K.H
index 3ffbfa0a..d0d3eeda 100644
--- a/Source/PeleLMeX_K.H
+++ b/Source/PeleLMeX_K.H
@@ -345,13 +345,13 @@ compute_divu<pele::physics::eos::Manifold>(
 
   auto eos = pele::physics::PhysicsType::eos(eosparm);
   eos.RY2RRinvY(massdens, rho, rhoinv, y);
-  amrex::Real derivs[NUM_SPECIES - 1] = {0.0};
+  amrex::Real derivs[MANIFOLD_DIM] = {0.0};
   rho /= 1.0e3; // Unit conversion MKS -> CGS
   eos.RY2dRdY(rho, y, derivs);
 
   divu(i, j, k) = 0.0;
 
-  for (int n = 0; n < NUM_SPECIES - 1; n++) {
+  for (int n = 0; n < MANIFOLD_DIM; n++) {
     derivs[n] *= 1.0e3; // Unit conversion CGS -> MKS
     divu(i, j, k) += derivs[n] * specDiff(i, j, k, n);
     if (do_react != 0) {