fix error for divu calc for manifold when its not being used with Nul…

…l mech
AMReX-Combustion · Dec 19, 2024 · 9372a81 · 9372a81
1 parent 652f0ea
commit 9372a81
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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) {