When the phase tracer doesn't function perfectly, the code sometimes returns an incorrect Jouguet velocity and/or LTE wall velocity in hydrodynamics.
Even if the solution of EOM with out-of-equilibrium contributions isn't affected by this, its error bar will be affected, since vwLTE enters there.
Maybe we can replace this vwLTE in the truncation error estimate by the LTE velocity of the template model whenever vwLTE = 1?
When the phase tracer doesn't function perfectly, the code sometimes returns an incorrect Jouguet velocity and/or LTE wall velocity in hydrodynamics.
Even if the solution of EOM with out-of-equilibrium contributions isn't affected by this, its error bar will be affected, since vwLTE enters there.
Maybe we can replace this vwLTE in the truncation error estimate by the LTE velocity of the template model whenever vwLTE = 1?