diff --git a/src/HydroContact/test_hydro_contact.cpp b/src/HydroContact/test_hydro_contact.cpp index 464df5774..2d26adbbd 100644 --- a/src/HydroContact/test_hydro_contact.cpp +++ b/src/HydroContact/test_hydro_contact.cpp @@ -117,8 +117,9 @@ void RadhydroSimulation::computeReferenceSolution(amrex::MultiFa stateExact(i, j, k, HydroSystem::x1Momentum_index) = rho * vx; stateExact(i, j, k, HydroSystem::x2Momentum_index) = 0.; stateExact(i, j, k, HydroSystem::x3Momentum_index) = 0.; - stateExact(i, j, k, HydroSystem::energy_index) = quokka::EOS::ComputeEintFromPres(rho, P) + 0.5 * rho * (vx * vx); - stateExact(i, j, k, HydroSystem::internalEnergy_index) = quokka::EOS::ComputeEintFromPres(rho, P); + stateExact(i, j, k, HydroSystem::energy_index) = + quokka::EOS::ComputeEintFromPres(rho, P) + 0.5 * rho * (vx * vx); + stateExact(i, j, k, HydroSystem::internalEnergy_index) = quokka::EOS::ComputeEintFromPres(rho, P); }); }