initial_conditions

Creates the structs for the pdf and the velocity-space moments

Take the full ion distribution function, calculate the moments, then normalise and shift to the moment-kinetic grid.

Uses input value of f and modifies in place to the normalised distribution functions. Input density, upar, ppar, and vth are not used, the values are overwritten with the moments of f.

Inputs/outputs depend on z, vperp, and vpa (should be inside loops over species, r)

Take the full neutral-particle distribution function, calculate the moments, then normalise and shift to the moment-kinetic grid.

Uses input value of f and modifies in place to the normalised distribution functions. Input density, upar, ppar, and vth are not used, the values are overwritten with the moments of f.

Inputs/outputs depend on z, vzeta, vr and vz (should be inside loops over species, r)

Allocate arrays for distributions to be applied as boundary conditions to the pdf at various boundaries. Also initialise the Knudsen cosine distribution here so it can be used when initialising the neutral pdf.

Allocate arrays for pdfs

initialise the electron density

initelectronpdfoverdensityandboundaryphi initialises the normalised electron pdf = pdfe * vthe / dense and the boundary values of the electrostatic potential phi; care is taken to ensure that the parallel boundary condition is satisfied; NB: as the electron pdf is obtained via a time-independent equation, this 'initital' value for the electron will just be the first guess in an iterative solution

initialise the electron parallel flow density

initialise the electron thermal speed profile. for now the only initialisation option for the temperature is constant in z. returns vth0 = sqrt(2*Ts/Te)

creates the normalised pdfs and the velocity-space moments and populates them with a self-consistent initial condition

for now the only initialisation option for the temperature is constant in z returns vth0 = sqrt(2Ts/ms) / sqrt(2Te/ms) = sqrt(Ts/Te)