From feca41237becfe1686167224b8144ea298149c54 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Wed, 8 Jan 2025 09:54:21 -0500 Subject: [PATCH] fix the perturbation amplitude --- .../simple_convection/problem_initialize_state_data.H | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/Exec/hydro_tests/simple_convection/problem_initialize_state_data.H b/Exec/hydro_tests/simple_convection/problem_initialize_state_data.H index c83e36ba30..5f56630bd9 100644 --- a/Exec/hydro_tests/simple_convection/problem_initialize_state_data.H +++ b/Exec/hydro_tests/simple_convection/problem_initialize_state_data.H @@ -66,8 +66,8 @@ void problem_initialize_state_data (int i, int j, int k, amrex::Real cs = std::sqrt(eos_rp::eos_gamma * eos_state.p / state(i,j,k,URHO)); - state(i,j,k,UMX) = state(i,j,k,URHO) * (2.0 * amrex::Random(engine) - 1.0_rt); - state(i,j,k,UMY) = state(i,j,k,URHO) * (2.0 * amrex::Random(engine) - 1.0_rt); + state(i,j,k,UMX) = state(i,j,k,URHO) * 0.05_rt * cs * (2.0 * amrex::Random(engine) - 1.0_rt); + state(i,j,k,UMY) = state(i,j,k,URHO) * 0.05_rt * cs * (2.0 * amrex::Random(engine) - 1.0_rt); state(i,j,k,UMZ) = 0.0_rt; state(i,j,k,UEDEN) += 0.5_rt * (state(i,j,k,UMX) * state(i,j,k,UMX) +