Great! I am glad you tested it.

Please review also my other changes. Note that the signature of flux_solution() has changed. Before:

flux_solution(PARTICLE_FLUX_e::T, STRESS_TOR_i::T, ENERGY_FLUX_e::T, ENERGY_FLUX_i::T)

and now

flux_solution(
    ENERGY_FLUX_e::T,
    ENERGY_FLUX_i::T,
    PARTICLE_FLUX_e::T,
    PARTICLE_FLUX_i::Vector{T},
    STRESS_TOR_i::T)

@orso82 I tested a few different cases with each of the models (:hirshmansigmar, :changhinton, :neo) and they give the same results as before for electron energy/particle and total ion energy fluxes, so this is good to go from my side.

I guess on the FUSE side you'll update things such that the neoclassical actor writes the ion_particle_flux solution to dd?