B.5 Calculation of

So far the only restriction over the radial coordinate ρ is that it is a monotonic flux function going from 0 at the center to 1 at the edge. In the code we take ρ = r ∕a_p where r is defined as a function of and a_p = r. Using

(144)

we find

(145)