dist_fn_arrays Module

g_adjust transforms between representations of perturbed dist'n func'n. = g_wesson J0(Z) - q phi/T F_m where <> = gyroaverage g_gs2 = g_wesson - q phi/T J0(Z) F_m - m v_||^2/T B_||/B J1(Z)/Z F_m For numerical convenience the GS2 variable g uses the form g_gs2. g_wesson (see Wesson's book, Tokamaks) is often a more convenient form: e.g. for handling collisions, calculating v-space moments in real space.

To transform gnew from g_gs2 to g_wesson form: call g_adjust(gnew,phinew,bparnew,fphi,fbpar) or call g_adjust(gnew,phinew,bparnew, to_g_gs2 = .false.) or transform from gnew from g_wesson to g_gs2 form: call g_adjust(gnew,phinew,bparnew,-fphi,-fbpar) or call g_adjust(gnew,phinew,bparnew, to_g_gs2 = .true. ) CMR, 17/4/2012:

Components

Return the offset between g and h at a specific set of indices

Arguments

Return Value complex

A small helper function which constructs and returns the string ' : '. Intended for use in producing the identifier string passed to the check routines of this module. Might be worth finding a better location for this generic utility.

Arguments

Return Value character(len=:), allocatable

Transform between g_gs2 and g_wesson with direction indicated by passed logical [to_g_gs2].

Arguments

Transform between g_gs2 and g_wesson with direction indicated by passed floats [facphi] and [facbpar].

Arguments

This routine checks trapped particle bounce conditions: g(thetab,1:ik,it,il,ie,is)=g(thetab,2:ik,it,il,ie,is) and flags fractional errors exceeding a threshold tolerance, tol. CMR, 3/10/2013:

Arguments

Assign $g$ a value corresponding to $h$ (g_wesson of CMR's note in g_adjust) of $h = 0$ in order to be consistent with g_adjust. The correct function call to set h = 0 is call set_h_zero(g, phi, bpar, iglo)

Arguments

Check if the passed distribution function shaped array is non-zero in the forbidden region. If it is then trigger an abort.

Arguments

Check if the passed distribution function shaped array has the same value at bounce points for both signs of v||. This provides similar functionality to check_g_bouncepoints but is intended to offer a quieter approach which could be used more routinely (i.e. throughout a run).

Arguments