diff --git a/Source/PeleLMeX_Eos.cpp b/Source/PeleLMeX_Eos.cpp index eb0af40e..3ce14f19 100644 --- a/Source/PeleLMeX_Eos.cpp +++ b/Source/PeleLMeX_Eos.cpp @@ -135,27 +135,27 @@ PeleLM::calcDivU( } else if (flagfab.getType(bx) != FabType::regular) { // EB containing // boxes amrex::ParallelFor( - bx, [rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, extRhoH, divu, - use_react, flag, + bx, [rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, + extRhoH, divu, use_react, flag, leosparm] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { if (flag(i, j, k).isCovered()) { divu(i, j, k) = 0.0; } else { compute_divu( - i, j, k, rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, extRhoH, - divu, use_react, leosparm); + i, j, k, rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, + extRhoY, extRhoH, divu, use_react, leosparm); } }); } else #endif { amrex::ParallelFor( - bx, - [rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, extRhoH, divu, - use_react, leosparm] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { + bx, [rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, + extRhoH, divu, use_react, + leosparm] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { compute_divu( - i, j, k, rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, extRhoY, extRhoH, - divu, use_react, leosparm); + i, j, k, rhoY, rhoH, T, SpecD, Fourier, DiffDiff, r, extRho, + extRhoY, extRhoH, divu, use_react, leosparm); }); } } diff --git a/Source/PeleLMeX_K.H b/Source/PeleLMeX_K.H index acca8f79..e02228ea 100644 --- a/Source/PeleLMeX_K.H +++ b/Source/PeleLMeX_K.H @@ -242,24 +242,23 @@ compute_divu( } // Additional terms related to external forcing on density - amrex::Real hExtRho = rhoinv * rhoH(i,j,k) * extRho(i, j, k); + amrex::Real hExtRho = rhoinv * rhoH(i, j, k) * extRho(i, j, k); amrex::Real yExtRho[NUM_SPECIES] = {0.0}; for (int n = 0; n < NUM_SPECIES; n++) { yExtRho[n] = y[n] * extRho(i, j, k); - } + } amrex::Real denominv = 1.0_rt / (rho * cpmix * T(i, j, k)); - divu(i, j, k) = - denominv * (specEnthDiff(i, j, k) + tempDiff(i, j, k) - + extRhoH(i, j, k) - hExtRho); + divu(i, j, k) = denominv * (specEnthDiff(i, j, k) + tempDiff(i, j, k) + + extRhoH(i, j, k) - hExtRho); for (int n = 0; n < NUM_SPECIES; n++) { - amrex::Real specTerm = + amrex::Real specTerm = specDiff(i, j, k, n) + extRhoY(i, j, k, n) - yExtRho[n]; if (do_react != 0) { specTerm += rhoYdot(i, j, k, n); } - divu(i, j, k) += specTerm * (mwtinv[n] * Wbar * rhoinv - hi[n] * denominv) - + rhoinv * extRho(i, j, k); + divu(i, j, k) += specTerm * (mwtinv[n] * Wbar * rhoinv - hi[n] * denominv) + + rhoinv * extRho(i, j, k); } }