pack_energy_diagnostics Subroutine

   subroutine pack_energy_diagnostics(energy, diag)
      type(energy_fields_t), intent(in) :: energy
      real(rk), intent(out) :: diag(12)

      diag = zero
      diag(1) = energy%last_conductive_boundary_flux_out
      diag(2) = real(energy%last_conductive_boundary_flux_available, rk)
      diag(3) = energy%cumulative_boundary_rho_h_advective_flux_out
      diag(4) = energy%cumulative_boundary_rho_h_conductive_flux_out
      diag(5) = energy%cumulative_qrad_integral
      diag(6) = energy%cumulative_rho_species_hdiff_integral
      diag(7) = real(energy%cumulative_energy_budget_available, rk)
      diag(8) = energy%last_energy_update_delta_rate_integral
      diag(9) = energy%last_energy_update_rhs_integral
      diag(10) = energy%last_energy_update_balance_defect
      diag(11) = energy%relative_last_energy_update_balance_defect
      diag(12) = energy%last_operator_consistent_rho_h_integral
   end subroutine pack_energy_diagnostics