critical_energy exponent of avg_cmr should be 2/3 and not 3/2
|
Ec = 14.8 .* mf .* Te .* avg_cmr .^ 1.5 |
the function right above is correct:
|
Ec = 14.8 * mf * Te * avg_cmr^(2.0 / 3.0) |
See: T.H. Stix, "Heating of Toroidal Plasmas by Neutral Injection", Plasma Physics 14 (1972), equation 15

critical_energy exponent of avg_cmr should be 2/3 and not 3/2
IMAS.jl/src/physics/fast.jl
Line 183 in c42fbdc
the function right above is correct:
IMAS.jl/src/physics/fast.jl
Line 171 in c42fbdc
See: T.H. Stix, "Heating of Toroidal Plasmas by Neutral Injection", Plasma Physics 14 (1972), equation 15