Zhi Chen

I agree that the addition should be removed. #2468 makes `mom_flux_has_p` to return true instead of false for 1D cartesian, which is certainly the correct behavior. So now pressure is...

what was the value you're getting after moving `k_B` to the end?

I think the problem is with nse state mass fraction in the previous timestep. If it overflows we get e500, but we get e297, which means a different NSE composition,...

hmm, never mind. I think it never really returns those answers as the solution.

In `nse_derivs` we need the derivatives to calculate energies. And if we were to do the same, we will need Ydots if we calculate energy using `e_nuc = -N_A *...

just to reiterate, so we want to do `e_nuc = -N_A * c^2 * dY_k/dt * m_k`. But using only weak rate contributions in `dY_k/dt` because in NSE, strong reactions...

Detonation with this pr:  but it was a lot slower. Took ~400s runtime. For comparison, I could run the same problem about ~100s with NSE before, and its about...

Detonation with this pr also gets statements like: `burn entered NSE during integration (after 117 steps), zone = (227, 0, 0) recovering burn failure in NSE, zone = (227, 0,...

I realized that the current version on development is likely missing a factor of `1/dt` in the definition of `enuc_dot`, but it seems to produce consistent results when I don't...

Update with some progress: With current version, I have: ``` density 3.575417e+19 1.942 9.303788e+18 1.980 2.35768e+18 1.601 7.773908e+17 xmom 1.804038e+28 2.065 4.311366e+27 1.983 1.090598e+27 1.902 2.917623e+26 ymom 1.803993e+28 2.065 4.311497e+27...