Ion flux estimation problem?

[avdeevag]

I noticed the following case which failed, that is a bit suspicious - it works well with other saturation rules and the total ion energy looks fine

chk = FUSE.Checkpoint();

ini, act=FUSE.case_parameters(:ITER, init_from=:scalars, boundary_from=:scalars);

dd = IMAS.dd();
FUSE.init(dd, ini, act);

chk[:init] = dd,ini,act;

R6="sat0_es_d3d_simple20";

dd_R1,ini_R1,act_R1 = chk[:init];

#----------------------

act_R1.ActorTGLF.user_specified_model = R6

act_R1.ActorCoreTransport.model = :FluxMatcher


FUSE.ActorPFdesign(dd_R1, act_R1);
FUSE.ActorStationaryPlasma(dd_R1, act_R1);

FUSE.plot_plasma_overview(dd_R1)