Evolve ion (impurity) density in ActorFluxMatcher

[nanshi1177]

Using ActorFluxMatcher to evolve impurity density profile as follows:

act.ActorTGLF.sat_rule = :sat1
act.ActorTGLF.lump_ions = false
act.ActorFluxMatcher.optimizer_algorithm =:simple
act.ActorTGLF.electromagnetic = false
act.ActorFluxMatcher.max_iterations=120
act.ActorFluxMatcher.rho_transport = 0.2:0.1:0.8
act.ActorFluxMatcher.step_size = 0.1
#act.ActorFluxMatcher.evolve_densities = :flux_match
act.ActorFluxMatcher.evolve_densities = Dict(
    :DT => :quasi_neutrality,
    :DT_fast => :fixed,
    :He4 => :fixed,
    :He4_fast => :fixed,
    :Kr83 => :flux_match,
    :electrons => :flux_match)
act.ActorTGLF.model = :TGLF
#act.ActorTGLF.user_specified_model = "sat0_es_d3d"

FUSE.ActorHCD(dd,act);
act.ActorNeoclassical.model = :hirshmansigmar
FUSE.ActorFluxMatcher(dd,act;verbose=true,do_plot=true)

The output shows below: