Adding internal energies to the output of the conservative variables (6-equation model)

src/post_process/m_start_up.fpp

@@ -419,6 +419,21 @@ contains
 
         end if
 
+        ! Adding the individual energies to the formatted database file --------
+        if (model_eqns == 3 .and. cons_vars_wrt) then
+            do i = 1, num_fluids
+
+                q_sf = q_cons_vf(i + intxb - 1)%sf(-offset_x%beg:m + offset_x%end, &
+                                                   -offset_y%beg:n + offset_y%end, &
+                                                   -offset_z%beg:p + offset_z%end)
+
+                write (varname, '(A,I0)') 'alpha_rho_e', i
+                call s_write_variable_to_formatted_database_file(varname, t_step)
+
+                varname(:) = ' '
+            end do
+        end if
+
         !Adding Energy cascade FFT
         if (fft_wrt) then