fix the EOS call

AMReX-Astro · Dec 5, 2023 · bcf379b · bcf379b
1 parent aef3a2f
commit bcf379b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 5 deletions.
diff --git a/Exec/reacting_tests/nse_test/convergence_simplified_sdc_w_vel.sh b/Exec/reacting_tests/nse_test/convergence_simplified_sdc_w_vel.sh
@@ -8,6 +8,7 @@ RUNPARAMS="
 castro.time_integration_method=3
 problem.u0=1.e8
 problem.v0=1.e8
+integrator.nse_deriv_dt_factor=0.01
 "
 
 mpiexec -n 8 ${EXEC} inputs.64 ${RUNPARAMS} >& /dev/null

diff --git a/Source/hydro/trace_ppm_T.cpp b/Source/hydro/trace_ppm_T.cpp
@@ -72,7 +72,6 @@ Castro::trace_ppm_T(const Box& bx,
 
   GpuArray<int, NQSRC> do_source_trace;
 
-
   for (int n = 0; n < NQSRC; n++) {
     do_source_trace[n] = 0;
 
@@ -151,7 +150,7 @@ Castro::trace_ppm_T(const Box& bx,
     eos_state.rho = rho;
     eos_state.T = T;
     for (int n = 0; n < NumSpec; ++n) {
-        eos_state.xn[n] = q_arr(i,j,k,QFS+n);
+        eos_state.xn[n] = U_arr(i,j,k,UFS+n) * rho_inv_arr(i,j,k);
     }
 
     eos(eos_input_rt, eos_state);
@@ -172,7 +171,6 @@ Castro::trace_ppm_T(const Box& bx,
     Real sm;
     Real sp;
 
-
     // reconstruct density
 
     Real Ip_rho[3];
@@ -218,7 +216,6 @@ Castro::trace_ppm_T(const Box& bx,
     ppm_reconstruct(s, flat, sm, sp);
     ppm_int_profile_single(sm, sp, s[i0], un, dtdx, Ip_utt_1, Im_utt_1);
 
-
     // source terms -- we only trace if the terms in the stencil are non-zero
     int do_trace;
 
@@ -261,7 +258,6 @@ Castro::trace_ppm_T(const Box& bx,
 
     // T would go here, but we don't currently fill temperature sources
 
-
     // transverse velocities
 
     Real Ip_src_ut_1 = 0.0_rt;