chalmersplasmatheory/DREAM

Bug with setting an initial current density when a runaway distribution function is used

Opened this issue · 0 comments

When setting an initial RE avalanche distribution function with f_re.setInitialAvalancheDistribution() and an initial total current density j_tot = j_ohm + j_re set with j_ohm.setInitialProfile(), the initial j_tot is instead exactly equal to j_ohm. This is despite having an RE current dominating the initial total current.

This is as if the integral to determine the RE current from the distribution function is posterior to the j_tot sum.