Initialise arrays that are needed during the main time advance loop in the field solve for the flux tube simulations only These are initialised once and then used throughout the rest of the simulation gamtot = 1 - gamma0 = Z^2*n/T * int e^(-v^2) * (1 - J0^2) dv If using adiabatic electrons then this factor is modified and we use gamtot3 which includes the Boltmann response
Layouts Arrays Parameters Grids Calculations Geometry
Initialise gamtot Gamtot does not depend on flux tube index, so only compute for one flux tube index Calculate (1- j0^2) for each kx,ky,z TODO-GA: remove maxwellian_normalisation flag Integrate over vpa mu with weight = Z^2*n/T Avoid divide by zero when kx=ky=0 We do not evolve this mode, so the value is irrelevant Initialise gamtot3 Avoid divide by zero for kx=ky=0 mode, which we do not need anyway