diff --git a/examples/structured_2d_dgsem/elixir_euler_peng_robinson_transcritical_mixing.jl b/examples/structured_2d_dgsem/elixir_euler_peng_robinson_transcritical_mixing.jl
index de37773e27a..79135be6f9b 100644
--- a/examples/structured_2d_dgsem/elixir_euler_peng_robinson_transcritical_mixing.jl
+++ b/examples/structured_2d_dgsem/elixir_euler_peng_robinson_transcritical_mixing.jl
@@ -97,9 +97,12 @@ alive_callback = AliveCallback(analysis_interval = analysis_interval)
 
 stepsize_callback = StepsizeCallback(cfl = 0.25)
 
+save_solution = SaveSolutionCallback(interval = 500,
+                                     solution_variables = cons2thermo)
+
 callbacks = CallbackSet(summary_callback,
                         analysis_callback, alive_callback,
-                        stepsize_callback)
+                        stepsize_callback, save_solution)
 
 ###############################################################################
 # run the simulation
diff --git a/examples/tree_1d_dgsem/elixir_euler_nonideal_density_wave_FVO2.jl b/examples/tree_1d_dgsem/elixir_euler_nonideal_density_wave_FVO2.jl
index b401ea19021..6f19340a5b3 100644
--- a/examples/tree_1d_dgsem/elixir_euler_nonideal_density_wave_FVO2.jl
+++ b/examples/tree_1d_dgsem/elixir_euler_nonideal_density_wave_FVO2.jl
@@ -80,9 +80,12 @@ alive_callback = AliveCallback(analysis_interval = analysis_interval)
 
 stepsize_callback = StepsizeCallback(cfl = 0.5)
 
+save_solution = SaveSolutionCallback(interval = 500,
+                                     solution_variables = cons2thermo)
+
 callbacks = CallbackSet(summary_callback,
                         analysis_callback, alive_callback,
-                        stepsize_callback)
+                        stepsize_callback, save_solution)
 
 ###############################################################################
 # run the simulation
diff --git a/src/equations/nonideal_compressible_euler_1d.jl b/src/equations/nonideal_compressible_euler_1d.jl
index b8c7990c415..1fd8470a12c 100644
--- a/src/equations/nonideal_compressible_euler_1d.jl
+++ b/src/equations/nonideal_compressible_euler_1d.jl
@@ -91,6 +91,10 @@ equation of state routines are assumed to be evaluated in terms of `V` and `T`.
     return SVector(V, v1, T)
 end
 
+varnames(::typeof(cons2thermo), ::NonIdealCompressibleEulerEquations1D) = ("V",
+                                                                           "v1",
+                                                                           "T")
+
 # Calculate 1D flux for a single point
 @inline function flux(u, orientation::Integer,
                       equations::NonIdealCompressibleEulerEquations1D)
diff --git a/src/equations/nonideal_compressible_euler_2d.jl b/src/equations/nonideal_compressible_euler_2d.jl
index 6bcbf368b83..9210ee35461 100644
--- a/src/equations/nonideal_compressible_euler_2d.jl
+++ b/src/equations/nonideal_compressible_euler_2d.jl
@@ -72,6 +72,11 @@ varnames(::typeof(cons2prim), ::NonIdealCompressibleEulerEquations2D) = ("rho",
                                                                          "v2",
                                                                          "p")
 
+varnames(::typeof(cons2thermo), ::NonIdealCompressibleEulerEquations2D) = ("V",
+                                                                           "v1",
+                                                                           "v2",
+                                                                           "T")
+
 # Calculate flux for a single point
 @inline function flux(u, orientation::Integer,
                       equations::NonIdealCompressibleEulerEquations2D)